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ABSTRACT 

(Distribution  Limitati-m  Statement  A) 

Experimental  data  from  two  shock  tube/shock  tunnel  facilities 
capable  of  simulating  the  operational  conditions  of  continuous  flow 
chemical  laser  devices  are  presented.  One  facility  is  specifically 
designed  for  reproducing  the  reaction  mixing  region  developed  when 
hydrogen  is  injected  through  various  sonic  ports  into  a  supersonic 
fluorine  carrying  stream;  the  other  facility  simulates  the  parallel 
shear  layer  mixing  developed  between  a  supersonic  fluorine  carry¬ 
ing  stream  and  an  adjacent  supersonic  hydrogen  carrying  stream. 
Established  experimental  techniques,  including  the  Mach-Zehnder 
interferometer,  rapid  scanning  infrared  spectrometer,  pitot  pres¬ 
sure  probes,  and  thermocouple  temperature  probes,  have  been  em¬ 
ployed  to  map  flow  conditions  within  the  test  section  mixing  layer. 

A  theoretical  model  of  the  parallel  shear  flow  turbulent  mixing  re¬ 
gion  is  described  aid  results  from  this  model  are  compared  with 
measured  monochromatic  fringe-shift  profiles.  Reasonably  good 
agreement  has  been  demonstrated  between  theoretical  predictions  and 
the  experimental  profiles  for  two  widely  different  experimental  sit¬ 
uations  . 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 


The  use  of  gas-phase  chemical  reactions  in  high-speed  flow 
systems,  both  as  energy  sources  and  as  mechanisms  for  achieving 
population  inversion,  has  introduced  complex  fluid  dynamic  phenom¬ 
ena  into  the  technology  of  high-energy  laser  systems.  The  complex¬ 
ity  arises  from  the  coupling  of  gas  dynamics  and  nonequilibrium  gas- 
phase  chemical  reactions.  It  is  now  apparent  that  the  design  of 
satisfactory  reaction/mixing  regions  for  an  HF  chemical  laser  will 
require  optimization  of  the  plenum  and  nozzle  block  designs  in  terms 
of  the  gas  dynamics  and  chemical  kinetics. 

The  results  of  a  combined  theoretical  and  experimental  inves¬ 
tigation  of  the  supersonic  mixing  of  hydrogen,  fluorine,  and  helium 
diluent  are  presented  in  this  report.  The  principle  objective  of 
this  effort  was  to  demonstrate  the  feasibility  of  developing  an 
experimentally  verified  analytical  procedure  for  evaluating  the 
mixing  legion,  such  as  would  exist  in  the  optical  cavity  of  a  laser 
device,  in  terms  of  its  potential  for  producing  laser  power.  The 
experimental  program  used  two  separate  facilities  to  investigate, 
first,  two  different  methods  for  injecting  hydrogen  into  a  super¬ 
sonic  fluorine  carrying  stream,  and  second,  mixing  of  two  parallel 
supersonic  fluorine  and  hydrogen  carrying  streams.  Diagnostic 
equipment  used  interchangably  with  each  of  the  two  basic  facilities 
included  a  Mach-Zehnder  interferometer,  rapid  scanning  infrared 
spectrometer,  pitot  and  static  pressure  probes,  and  a  thermocouple 
temperature  probe. 

The  configurations  that  were  tested  in  the  experimental  phase 

are: 

1)  Multiple-orifice  injection  of  hydrogen  between  adjacent 
supersonic  streams  of  partially  dissociated  (>  70  percent) 
fluorine  plus  helium  diluent; 

2)  Slot  (or  slit)  injection  of  hydrogen  between  adjacent 
supersonic  streams  of  partially  dissociated  (>  70  percent) 
fluorine  plus  helium  diluent; 

3)  Mixing  of  adjacent  parallel  supersonic  streams,  one  con¬ 
taining  partially  dissociated  fluorine  plus  diluent  and  the 
other  containing  hydrogen  plus  argon. 
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Using  the  multiple-orifice  hydrogen  injection  facility,  mono¬ 
chromatic  snapshot  interferograms  of  the  mixing  region  were  recorded 
for  15  different  test  conditions.  The  test  conditions  were  composed 
of  five  separate  injector  conditions  for  each  of  three  different 
exit  static  pressures  in  the  fluorine  nozzle.  Infrared  emission 
intensity  was  not  recorded  for  these  tests  due  to  poor  transmission 
of  the  viewing  windows  available  during  the  test  period.  Complete 
static  and  total  pressure  surveys  were  conducted  for  one  nozzle  exit 
pressure  and  three  injector  conditions,  using  pure  fluorine  in  the 
nozzle  reservoir.  The  results  of  these  surveys  verified  the  uniform 
quality  of  the  flow  field  at  the  nozzle  exit  plane. 

Monochromatic  snapshot  interferograms  were  also  obtained  for  the 
slot  hydrogen  injection  scheme  for  15  different  test  conditions. 

These  tests  comprised  five  different  injection  conditions  for  each 
of  three  fluorine  nozzle  exit  static  pressures.  Rapid-scan  infrared 
emission  intensities  were  recorded  through  a  sapphire  window  at 
25  stations  along  the  nozzle  axis  downstream  of  the  injector.  Two 
different  injector  supply  pressures  were  considered  for  a  constant 
value  of  the  nozzle  exit  static  pressure.  These  data  are  reported 
in  terms  of  average  radiating  power  density  at  2.79  pm  in  the  mix¬ 
ing  region.  Static  and  total  pressure  surveys  were  recorded  for  two 
different  injector  pressure  conditions  and  a  single  nominal  value 
of  the  nozzle  exit  static  pressure.  Using  a  tungsten/tungsten-26% 
rehenium  thermocouple,  a  temperature  survey  of  the  mixing  region  was 
attempted.  A  2-mil-dia  thermocouple  junction  did  record  an  average 
value  of  the  temperature  across  the  mixing  region,  but  the  slow 
response  4  msec)  of  this  probe  produced  questionable  results. 

A  1-mil  thermocouple  junction  was  obtained  late  in  the  program  and 
gave  a  better  time  response  (v  1  msec)  than  did  the  2-mil  wire. 
However,  it  was  very  fragile  and  failed  after  four  shots.  A 
shielded,  vented  version  of  this  probe  should  work  quite  well  and 
will  be  used  in  future  experiments  with  this  facility. 

The  data  obtained  from  the  parallel  stream  mixing  configuration 
were  somewhat  limited  in  extent  due  to  time  constraints  and  the 
stronger  emphasis  placed  on  testing  two  hydrogen  injection  concepts. 

A  total  of  nine  monochromatic  snapshot  interferograms  were  recorded 
of  the  shear  layer  mixing  zone.  For  one  of  these  tests,  helium  was 
substituted  for  hydrogen  to  illustrate  a  nonreacting  mixing  condition. 
Rapid-scan  infrared  emission  spectra  were  obtained  at  five  different 
stations  downstream  of  the  initial  point  of  mixing  for  each  of  three 
different  parallel  mixing  conditions.  Static  pressure  measurements 
corresponding  to  each  of  the  nozzle  exit  conditions  were  recorded. 
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The  major  theoretical  efforts  in  this  study  were  directed 
toward  predicting  the  profiles  of  those  flow  parameters  that  could 
be  observed  within  the  mixing  layer  in  configuration  3  of  the  exper¬ 
imental  phase.  A  previously  developed  analysis  of  the  turbulent 
mixing  of  hydrogen/oxygen  supersonic  shear  layers  was  extensively 
revised  and  improved  during  this  effort  in  an  attempt  to  accurately 
account  for  the  chemical  kinetic  and  mixing  parameters  unique  to 
the  hydrogen/fluorine  system.  This  analysis  and  the  associated 
computer  code  employ  an  explicit  finite  difference  technique  for 
solving  the  boundary  layer  type  equations  that  describe  chemically 
reacting  turbulent  mixing  layers.  Program  options  include  consid¬ 
eration  of  either  equilibrium  or  finite-rate  chemical  kinetics  and 
uniform  or  arbitrary  initial  flow-property  profiles.  This  theo¬ 
retical  mixing  analysis  technique  was  used  to  predict  test  section 
parameters  corresponding  to  two  experimental  configurations  for 
which  snapshot  monochromatic  interferograms  were  recorded.  Fringe- 
shift  profiles  were  computed  through  the  mixing  region  as  the  most 
direct  means  of  comparison  with  expeiimental  data.  A  nonequilibrium 
calculation  (finite-rate  chemistry)  was  performed  for  one  set  of 
experimental  conditions  and  both  nonequilibrium  and  frozen-flow 
(nonreacting  chemistry)  calculations  were  performed  for  a  second 
set  of  experimental  conditions.  In  both  of  these  cases,  the  experi¬ 
mental  fringe-shift  data  are  compared  with  the  theoretical  predic¬ 
tions. 

The  mixing  analysis  program  contains  a  similarity  parameter,  o 
(constant  related  to  the  mixing  layer  growth  rate) ,  whose  value 
can  be  estimated  from  previous  experimental  evidence.  To  date, 
only  one  such  estimated  value  of  a  has  been  used  in  these  calcu¬ 
lations.  As  a  result,  comparisons  reported  here  must  be  consid¬ 
ered  to  represent  the  zeroth  order  in  the  theoretical/data  corre¬ 
lation  procedure,  which  can  be  refined  with  additional  efforts. 

As  a  result  of  the  studies  reported  here,  some  insight  into 
coupling  between  gas  dynamic  processes  and  nonequilibrium  chemistry 
associated  with  high-energy  chemical  lasers  has  been  derived.  Many 
characteristics  of  operational  devices  have  been  observed  in  these 
studies,  under  conditions  that  lend  themselves  to  quantitative 
analysis.  The  capability  of  the  two-dimensional  shock  tube/shock 
tunnel  facilities  to  produce  two-dimensional  mixing  layers  over  a 
range  of  conditions  for  both  injection  and  parallel  stream  mixing 
has  been  demonstrated.  Monochromatic  snapshot  interferograms 
graphically  depict  the  growth  of  the  mixing  region  in  each  of  these 
facilities  and,  when  reduced,  result  in  detailed  profiles  of  average 
density  through  the  mixing  layer.  The  rapid  scanning  infrared  spec¬ 
trometer  has  produced  detailed  emission  spectra  for  selected  portions 
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of  the  mixing  region.  From  these  spectra,  the  total  number  density 
of  excited  hydrogen-fluoride  molecules  within  the  mixing  region  can 
be  deduced.  Pitot  probe  surveys  have  verified  the  uniformity  of  the 
nozzle  exit  flow  fields. 

A  significant  feature  of  our  experimental  approach  is  that  it 
does  produce  mixing  layers  that  lend  themselves  to  quantitative 
analysis.  This  has  been  demonstrated  by  applying  a  theoretical 
mixing  analysis  computer  program  to  the  test  conditions  correspond¬ 
ing  to  two  different  parallel  mixing  experiments.  Even  though  a 
complete  understanding  of  the  many  complex  problems  has  not  yet 
been  achieved,  the  ability  of  the  analysis  to  correctly  predict 
important  features  of  the  experimental  profiles  within  the  mixing 
layer  is  most  encouraging.  Past  experience  with  hydrogen/ oxygen 
mixing  layers  indicates  that,  with  some  adjustment  of  the  empirical 
parameters  within  the  analytical  program,  a  close  correlation  be¬ 
tween  measured  profiles  and  computed  profiles  can  be  obtained. 

Once  this  goal  has  been  reached,  the  program  can  be  used  to  extend 
the  results  over  a  much  wider  range  of  flow  parameters.  This 
achievement  will  be  a  valuable  asset  for  the  design  of  mixing 
regions  of  large-scale  continuous- flow  devices. 
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SECTION  II 


EXPERIMENTAL  FACILITY  DESCRIPTION 


The  two  experimental  facilities  used  in  this  program  are  both 
adaptations  of  test  equipment  developed  at  Martin  Marietta  for  ex¬ 
amining  the  turbulent  mixing  of  parallel  supersonic  flows  in  the 
presence  of  hydrogen/oxygen  reactions  as  encountered  in  air-aug¬ 
mented  propulsion  systems.  Each  of  these  two  facilities  uses  the 
shock  heated  gases  produced  at  the  end-wall  region  of  a  shock  tube 
as  supply  reservoirs  for  two-dimensional  nozzles.  In  the  injection 
mixing  facility,  the  nozzle  is  separated  into  two  half  nozzles  by 
a  splitter  plate  from  which  hydrogen  can  be  injected  into  the 
supersonic  fluorine  stream  by  different  means  at  the  downstream 
end  of  the  splitter  plate.  The  supersonic  parallel  mixing  shock 
tunnel  consists  of  two  adjacent  two-dimensional  nozzles  supplied 
by  gases  produced  within  the  shock  tube  end-wall  region.  The 
most  advantageous  characteristic  of  these  configurations  is  that 
they  produce  well-defined  two-dimensional  flow  fields  with  con¬ 
venient  provisions  for  varying  flow  parameters  such  as  test  section 
pressures,  temperatures,  degree  of  dissociation,  and  fuel-to- 
oxidizer  ratios.  The  mixing  region  for  each  facility  can  be  viewed 
through  optical  windows  mounted  in  the  test  section  region,  and 
the  two-dimensionality  of  the  flow  fields  provides  a  convenient 
means  fcr  mapping  the  mixing  regions  by  optical  techniques  such 
as  the  Mach-Zehnder  interferometer  and  a  rapid  scanning  spectro¬ 
meter.  Complete  details  for  each  of  these  basic  facilities,  the 
associated  instrumentation,  and  theoretical  analysis  is  given  in 
the  following  four  subsections. 

1.  INJECTION  SHOCK  TUNNEL 

The  principal  dimensions  of  the  injection  mixing  facility  are 
given  in  Figure  1.  The  driver  consists  of  an  8-ft-long  section 
of  circular  tubing  with  an  8-in.  internal  diameter  connected  by  a 
1-ft-long  transition  section  to  provide  a  smooth  interface  between 
the  circular  driver  and  the  square  driven  tube.  The  main  diaphragm 
station  is  located  at  the  downstream  ena  of  the  transition  station, 
and  a  mechanical  harpoon  is  used  to  initiate  the  rupture  of  the 
aluminum  diaphragm  at  the  appropriate  time.  The  driven  tube  con¬ 
sists  of  five  sections  of  5-in.‘  aluminum  tube  with  a  total  length 
of  37  ft,  and  terminates  at  the  two-dimensional  nozzle  block  as¬ 
sembly.  A  thin  (0.010  in.)  aluminum  diaphragm  seals  off  the  driven 
section  from  the  test  section  and  dump  tank  assembly.  A  mild  steel 
dump  tank  (20-in.  outside  diameter  by  15-ft-long)  is  connected  to 
the  downstream  end  of  the  aluminum  test  section  assembly. 


5 


6 


Because  of  the  particularly  noxious  and  corrosive  nature  of 
the  fluorine  test  gas,  special  handling  and  cleaning  procedures 
were  required.  Both  the  driven  shock  tube  sections  and  the  nozzle 
assembly  were  throughly  degreased  and  then  exposed  to  pure  fluorine 
gas  at  30  psia  for  at  least  1  hr.  This  well-established  procedure 
forms  an  aluminum  fluoride  coating  on  the  exposed  aluminum  surfaces 
by  a  process  where  the  fluorine  replaces  oxygen  in  the  initial 
aluminum  oxide  coating.  This  procedure  allows  the  driven  section 
to  be  filled  with  up  to  50  torr  of  pure  fluorine  without  a  mea¬ 
surable  decay  in  pressure  for  up  to  30  minutes  for  a  leak-tight 
assembly. 

The  fluorine  handling  procedure  used  with  each  of  these  two 
facilities  is  shown  schematically  in  Figure  2.  Bottles  containing 
fluorine  are  always  maintained  in  the  laboratory  at  subatmospheric 
pressures  so  that,  in  the  event  of  a  leak,  air  will  pass  into  the 
plumbing  system  rather  than  fluorine  escaping.  The  shock  tube  dump 
tank  assembly  has  a  volume  that  is  sufficiently  large  to  ensure 
that  the  final  resulting  pressure  at  the  conclusion  of  a  run  is 
again  subatmospheric.  As  indicated  in  Figure  2,  the  mechanical 
pumps  all  exhaust  throi_ah  a  common  line  that  passes  the  exhaust 
gas  through  a  charcoal  filter  before  it  is  released  to  the  atmo¬ 
sphere.  The  postshot  operational  procedure  requires  an  immediate 
evacuation  to  at  least  20  torr  and  a  subsequent  refilling  to  at¬ 
mospheric  pressure  with  dry  nitrogen.  This  mixture  is  evacuated 
a  second  time  to  10  torr  and  then  refilled  with  dry  nitrogen 
before  opening  the  assembly  for  diaphragm  replacement.  This  general 
cleaning  and  gas  handling  procedure  has  been  very  successful,  as 
indicated  by  almost  daily  incident-free  operation  for  more  than 
12  months. 

Details  of  the  Mach  5.0  nozzle  blocks  and  splitter  plate  as¬ 
sembly  are  shown  in  Figures  3  and  4,  and  Figure  5  is  a  photograph 
of  the  actual  assembly  with  one  side  plate  removed.  The  field 
of  view  covered  by  the  optical  windows  is  shown  outlined  in  the 
drawing  of  Figure  3,  and  the  actual  nozzle  contour  dimensions  are 
given  in  the  insert  tabulation  of  Figure  4. 

Figure  6,  a  photograph  of  the  thin  aluminum  test  section 
diaphragm  taken  after  a  test  shot,  shows  that  the  aluminum  is 
sheared  very  nicely  by  the  sharp  leading  edge  of  the  injector  split¬ 
ter  plate  and  rests  smoothly  against  the  leading  edges  of  the 
nozzle  blocks.  These  diaphragms  are  very  consistent  in  opening 
and  no  measurable  material  loss  has  ever  been  detected,  with  the 
following  exceptions.  On  two  different  test  occasions,  the  test 
section  diaphragm  was  torn  very  badly;  this  could  always  be  traced 
to  material  loss  of  the  driver  diaphragm. 


7 


Figure  2.  Gas-Handling  Schematic 


Figure  3.  Injection  Test  Section 


THROAT  STATION  «  0.000 
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Figure  4.  Mech  5.0  Nozzle  Block  Detail* 
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Figure  5.  Injection  Test  Section 


Opened  Test 


A  typical  test  shot  sequence  of  events  consists  of  evacuating 
the  three  separate  facility  chambers  as  indicated  in  Figure  1 
using  mechanical  vacuum  pumps.  The  shock  tube  driven  section  and 
the  dump  tank  plus  test  section  are  evacuated  to  a  pressure  level 
of  at  least  10  microns,  while  the  driver  section  is  evacuated  to 
a  pressure  of  at  least  10  torr.  The  desired  test  gas  (pure  fluorine 
or  a  50%  by  volume  mixture  of  fluorine  and  helium)  is  loaded  into 
the  driven  tube  and  the  initial  load  pressure  is  measured  with  one 
of  three  available  types  of  stainless  steel  diaphragm  pressure 
transducers.  The  most  sensitive  transducers  (Statham  PA731TC's, 

0  to  2  psia  range)  are  used  to  measure  the  lower  initial  load 
pressure  between  3  and  25  torr,  and  two  less  sensitive  trans¬ 
ducers  (Viatran  105' s  with  No.  603  signal  conditioner,  0  to  5 
and  0  to  15  psia  ranges)  are  used  for  higher  initial  load  pres¬ 
sures.  These  gages  are  fluorine  compatible,  at  least  at  these 
low  pressure  levels,  and  they  are  checked  frequently  against 
a  Wallace  Tiernan  absolute  pressure  dial  gage.  This  procedure 
assures  an  accuracy  in  the  initial  fill  pressures  of  at  least  ±1%. 

Immediately  after  filling  the  test  section,  the  dump  tank 
pump  valve  is  shut  off,  and  the  high-pressure  driver  gas  is  loaded 
into  the  driver  section  to  the  desired  level.  The  driver  pressure 
is  measured  with  a  Heise  pressure  gage  of  appropriate  range. 

Cnee  the  required  pressure  level  is  obtained,  the  mechanical  har¬ 
poon  is  released,  breaking  the  driver  diaphragm.  Immediately  after 
the  shot,  the  driven  section  mechanical  pump  valve  is  opened,  and 
the  residual  gas  is  pumped  through  the  charcoal  filter.  As  de¬ 
scribed  earlier,  the  final  pressure  in  the  entire  facility  is 
always  maintained  at  or  below  atmospheric  pressure  to  minimize  the 
fluorine  hazard.  This  entire  procedure  requires  less  than  2  min¬ 
utes  from  start  of  test  gas  fill  to  actual  test  shot.  The  leak 
rate  in  the  driven  tube  is  less  than  1  micron/minute,  thus  ensuring 
a  high  degree  of  purity  in  the  test  gas  sample. 

It  is  well  known  (Ref.  1)  that  conditions  in  the  quiescent 
gas  behind  the  reflected  shock  can  be  accurately  computed,  given 
the  initial  test  gas  conditions  (composition  and  pressure)  and 
incident  shock  wave  velocity.  Three  stations  (12  in.  apart)  are 
provided  for  measuring  the  initial  shock  velocity,  and  the  last 
station  is  located  3  in.  from  the  test  section  diaphragm  station. 
Shock  wave  arrival  at  each  of  these  stations  is  detected  by  one 
of  two  types  of  pressure  transducers  that  were  used  in  this  study, 
and  the  output  is  recorded  on  both  a  dual-beam  oscilloscope  and 
a  paralleled  Beckman  Time  Interval  Meter.  The  pressure  transducers 
used  for  detecting  shock  arrival  were  either  the  Piezotronics 
PCB113A21  gage  with  a  model  480A  power  supply  or  a  Kistler  type 
603A  gage  with  a  model  504  charge  amplifier.  Two  separate  time 
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intervals  were  recorded  with  the  three  pickup  stations  used. 

Hence,  shock  velocity  attenuation  could  be  measured  if,  in  fact, 
it  did  exist.  In  all  of  this  work,  time  intervals  accurate  to 
within  tl  Lsec  were  obtained,  and,  within  this  accuracy,  the 
shock  velocity  was  constant  over  the  span  for  which  it  was  mea¬ 
sured.  The  last  shock  velocity  station  (located  3  in.  from  the 
end  wall)  could  also  be  used  as  a  check  on  the  theoretically  pre¬ 
dicted  equilibrium  pressure  behind  the  reflected  shock.  Thin- 
film  heat  transfer  gages  were  also  used  as  shock  arrival  detectors 
in  the  early  phase  of  this  program,  but  were  never  completely 
satisfactory  because  the  fluorine  attacked  and  destroyed  the  sen¬ 
sitive  elements,  usually  after  six  to  eight  test  shots. 

Initial  test  gas  mixtures  containing  50%  (by  volume)  fluorine 
plus  helium  were  used  in  most  of  these  test  conditions,  and  these 
mixtures  were  prepared  and  allowed  to  stand  for  at  least  2  hr 
before  their  use  to  ensure  complete  mixing.  Some  of  the  first 
test  gas  mixtures  used  in  this  study  were  prepared  by  simply 
separately  adding  each  test  gas  to  the  driven  tube  section.  The 
resulting  shock  velocities  could  never  be  predicted  for  mixtures 
prepared  in  this  way,  and  it  was  concluded  that  improper  mixing 
was  responsible  for  mis  situation.  The  final  procedure  where 
the  gases  were,  preni  red  proved  much  more  satisfactory,  and  the 
resulting  measured  shock  velocities  were  in  excellent  agreement 
with  predicted  values. 

As  previously  described,  hydrogen  (helium  in  some  cases)  was 
injected  into  the  supersonic  fluorine  carrying  streams  from  ports 
located  at  the  downstream  end  of  the  splitter  plate.  Two  different 
splitter  plates,  one  containing  18  orifices  (0.020  in.  in  diameter) 
and  the  other  containing  a  series  of  slots  (0.0028  in.  wide  and 
covering  90%  of  the  available  exit  length)  were  used  in  this  study. 
A  detailed  sketch  o'  the  orifice  injector  splitter  plate  is  shown 
in  Figure  7,  and  Figure  8  is  a  photograph  of  this  splitter  plate 
looking  upstream  toward  the  nozzle  throat.  Figure  9  is  a  sketch 
and  photograph  of  the  slotted  hydrogen  injector  splitter  plate  with 
the  trailing  edge  tilted  to  expose  the  interior  construction.  The 
injectant  gas  flow  through  each  of  these  splitter  plates  is  acti¬ 
vated  by  a  microswitch  on  the  harpoon  mechanism  shaft  that  breaks 
the  driver  diaphragm.  Adjusting  this  switch  position  permits  the 
gas  flow  to  the  injector  to  be  timed  so  that  the  injection  flow 
starts  approximately  1  msec  before  start  of  fluorine  nozzle  flow. 
Pressure  on  the  injection  line  was  varied  from  50  to  150  psi  under 
no-flow  conditions. 
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0.235  i  0.001 


Figure  7.  Orifice  Splitter  Plate  Schematic 


Figure  8.  Orifice  Splitter  Plate 
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a.  Schematic 


b.  Photograph 


Figure  9.  Slotted  Splitter  Plate 


2.  PARALLEL  MIXING  SHOCK  TUNNEL 


This  facility  is  very  similar  to  the  basic  injection  shock 
tube/shock  tunnel  facility  except  that  part  of  the  driven  tube 
is  divided  into  two  parallel  chambers.  The  high-pressure,  high- 
temperature  gas  produced  behind  the  reflected  shock  wave  in  each 
chamber  is  then  expanded  through  two-dimensional  supersonic  nozzles, 
and  the  resulting  parallel  supersonic  streams  are  then  allowed 
to  mix.  The  m: xing  zone  can  be  observed  through  a  pair  of  optical 
quality  quartz  windows.  The  principal  dimensions  of  this  facility 
are  shown  in  Figure  10,  where  '  he  driver  and  dump  tank  chambers 
are  identical  with  those  described  for  the  injection  mixing  fa¬ 
cility. 

As  shown  in  Figure  10,  the  split  tube  assembly  consists  of 
two  sections  that  total  16  ft  long.  The  square  (not  divided) 
10-ft-long  buffer  chamber  is  used  as  a  conventional  shock  Lube 
section,  and  the  initial  shock  wave  generated  in  this  section 
simply  passes  through  the  thin  Teflon  diaphragm  (0.002  in.)  at 
Scation  A  almost  undisurbed.  This  shock  wave  then  simultaneously 
initiates  two  separate  shock  waves  in  the  divided  tube  section. 

Argon  has  been  used  in  the  buffer  section  of  this  facility  both 
for  convenience  and  because  a  heavy  monatomic  gas  provides  a 
higher  degree  of  stability  for  the  driver  and  driven  gas  inter¬ 
face  (Ref.  1). 

Figure  11  is  a  schematic  of  the  two  full  nozzles  used  for  these 
tests.  Each  of  these  nozzles  is  designed  for  an  ideal  gas  Mach 
number  of  4.5,  and  the  contours  are  different  because  of  the  dif¬ 
ferent  value  of  .■  (specific  heat  ratio  Cp/Cvj  for  the  two  different 

gases.  A  sketch  of  the  field  of  view  provided  by  the  optical 
window  is  shown  superimposed  on  this  figure.  A  photograph  of  the 
assembled  nozzle  blocks  and  test  section  is  shown  in  Figure  12, 
with  one  side  of  the  test  section  removed. 

The  gas  handling  procedure,  shock  arrival  stations,  and  instru¬ 
mentation  fox  this  facility  are  identical  with  those  described 
for  the  injection  mixing  facility*  Shock  arrival  stations  are 
located  in  both  chambers  used  as  supply  reservoirs  for  the  nozzles 
because  the  initial  load  conditions  and  shock  Mach  numbers  will 
generally  be  different  on  each  side  of  the  divided  chamber. 

The  divided  driven  sections  are  constructed  of  aluminum  by  cut¬ 
ting  the  5x5-in.'  tube  along  its  full  length  and  then  inserting  a 
1/4  in.  aluminum  plate  between  the  two.  A  continuous  weld  is  used 
to  seal  this  divider  plate,  and  great  care  is  taken  to  ensure  that 
no  leaks  exist  between  tubes. 
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Mach  4.5  Nozzle  Blocks  F„  Throat  —  0.62x5.000  in 
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Figure  11.  Parallel  Mixing  Test  Section  Schematic 


2 


Figure  12.  Parallel  Mixing  Test  Section 


3.  THEORETICAL  PERFORMANCE  ANALYSIS  OF  SHOCK  TUBE/SHOCK  TUNNEL 

The  basic  shock  tube  performance  analysis  in  both  the  incident 
and  reflected  shock  regions  was  based  on  treatment  of  the  shock 
heated  gases  as  equilibrium  mixtures  (Ref.  2).  The  nozzle  ex¬ 
pansion  processes  were  computed  assuming  an  equilibrium  sudden- 
freeze  model  (Ref.  3)  to  predict  the  exit  plane  thermodynamic  con¬ 
ditions.  A  summary  description  of  the  analytical  procedures  and 
typical  results  are  presented  in  this  subsection. 

To  be  assured  that  the  degree  of  ionization  can  be  neglected 
for  typical  conditions  behind  the  reflected  shock,  an  equilibrium 
composition  calculation  for  a  mixture  of  fluorine  molecules  and 
atoms  was  performed.  The  results  of  this  calculation,  assuming 
a  2000°K  gas  temperature  and  using  data  from  JANAF  tables  (Ref.  4), 
indicated  that  the  degree  of  ionization  is  much  less  than  0.17o. 
Therefore,  no  charged  particle  species  need  be  considered  in  the 
ana  lysis . 

Furthermore,  it  was  assumed  that  for  the  frozen  expansion  case 
considered  here,  all  thermodynamic  states  of  the  F2  and  F  could 
be  treated  by  ideal  gas  theory.  Thus,  for  the  equilibrium  reaction 

F^2F  (1) 


the  equilibrium  constant, 


y 

K  =  y1-  P 


K  ,  is  given  by 
P 


(2) 


where  is  the  mole  fraction  of  fluorine  atoms,  X^  is  the  mole 

fraction  of  fluorine  molecules,  and  P  is  the  total  pressure,  in 
atmospheres,  of  the  mixture.  The  equilibrium  constant  as  a  func¬ 
tion  of  temperature  for  Equation  (1)  is  obtained  from  reference 
4. 


For  the  mixture  we  can  write 
XF+XF 


(3) 


Substituting  Equation  (3)  into  Equation  (2)  and  solving  for  X^,  yields 


*F 


-  (%/p)  +  VftAf  +  *(VP) 


(4) 
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This 


The  positive  root  is  always  selected  because  0  <_  X^,  <_  1. 

calculation  proceeds  in  a  straightforward  manner  once  tempera¬ 
ture  and  pressure  are  specified.  The  helium  atoms  of  the  initial 
gas  mixture  will  always  be  in  the  ground  state,  and  there  will  be 
no  ionization  of  helium  at  this  temperature  level.  For  an  initial 
ideal  gas  mixture  containing  one  mole  of  F2  and  3  moles  of  helium, 
the  moles  of  helium  can  be  expressed  as 


(5) 


where  X  is  the  mole  fraction  of  helium  in  the  initial  fluorine/ 

helium  mixture.  The  material  balance  for  the  equilibrium  calcu¬ 
lation  given  by  Equation  (1)  may  be  determined  as  follows: 

Basis:  1  mole  of  F2  gas  (initially  undissociated) 

S  =  moles  of  He/mole  of  F2 
a  =  moles  of  F2  dissociated 

Table  I  summarizes  the  expressions  for  the  resultant  composition 
of  the  initial  mixture,  as  well  as  the  final  equilibrium  mixture. 


Table  I  EQUILIBRIUM  COMPOSITION  OF  F2/He  MIXTURES 


Component 

No.  of  Moles 

Mole  Fraction 

Initial 

Equilibrium 

He 

U 

E 

V  =  8 

nie  1  +  a  +  E 

f2 

1 

1  -  a 

Y  _  - L_I_i - 

^2  1  +  a  +  B 

2a 

F 

0 

2  a 

^  1  +  a  +  3 

Thus,  the  equilibrium  constant  for  the  fluorine  dissociation 
reaction  is  given  by 

X2  (2a)2 

V  X^  P  =  (1  +  a  +  S)(l  -  a)  P 


(6) 
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(7) 


This  expression  leads  to  the  quadratic  equation, 

/  K  \  6K  K 

\4  +  f-j  u2  +  a  -  ^  (1  +  3)  =  0 

The  value  of  a  is,  again,  obtained  from  the  positive-root  solution, 


Note  that  when  2*0,  Equation  (6)  reduces  to  the  result  ob¬ 
tained  for  a  pure  fluorine  gas.  The  mole  fraction  of  each  compon¬ 
ent  in  the  mixture  is  obtained  by  substituting  a  from  Equation  (8) 
for  each  value  of  K^,  which  is  determined  by  the  temperature,  the 

pressure,  and  the  initial  composition,  into  the  appropriate  expres¬ 
sions  from  Table  I.  Application  of  this  analysis  over  the  temper¬ 
ature  range  (300  to  2000°K)  and  pressure  range  (2  to  5  atmospheres) 
appropriate  to  the  experimental  study  demonstrates  that  the  impor¬ 
tant  chemical  species  in  the  gas  models  are  F 2,  F,  and  He  for  both 
the  reflected  shock  region  of  the  injection  shock  tube  and  the 
fluorine  side  of  the  parallel  mixing  shock  tube. 

Perhaps  the  most  significant  result  of  these  considerations 
is  that  the  flow  fields  in  the  injection  shock  tube  and  the  par¬ 
allel  mixing  shock  tube  can  be  treated  by  ideal  dissociating  gas 
theory  in  the  reflected  shock  region  before  the  nozzle  expansion. 

This  means  that  the  equilibrium  thermochemical  state  of  the  mix¬ 
tures  in  the  reflected  shock  region  can  be  determined  by  solving 
the  well-known  shock  crossing  relations  (Ref.  1)  for  conservation 
of  mass,  momentum,  and  energy,  along  with  satisfaction  of  the  law 
of  mass  action  for  an  ideal  dissociating  gas  and  the  appropriate 
thermal  and  caloric  equations  of  state.  These  equations  were 
solved  numerically  for  a  range  of  initial  driven  tube  conditions 
and  shock  strengths  encompassing  the  desired  experimental  test  con¬ 
ditions  to  give  quantitative  information  about  the  thermo chemical 
state  of  the  reflected  shock  region.  The  general  forms  of  the 
equations  that  were  solved  for  the  postshock  thermochemical  state 
of  the  gas  are 

Conservation  of  mass 

caua  "  Db%  (9) 
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Conservation  of  momentum 


C  U  *  +  P  =  0,0/  +  P, 
a  a  a  bo  b 

Conservation  of  energy 


ha  +  V  *  h,  +  V 
Thermal  equation  of  state 
P  -  (1  +  a)  oRT 
Caloric  equation  of  state 

Tb  -  [\  -  *b D) 1  (4  +  “b)  B1 


From  the  law  of  mass  action  (Ref.  2),  the  dissociation  fraction, 
a,  is  given  by 

'6d 


(10) 


(11) 


(12) 


(13) 


(14) 


where  6^  is  the  characteristic  dissociation  temperature  for  the 
diatomic  species  (in  this  case  F2)  and  the  constant  p  is  given  by 

/n  c  \  0 


where 

m  =  mass  of  the  atomic  species 

6^  =  characteristic  rotational  temperature 

9  =  characteristic  vibrational  temperature 

QgJc  =  electronic  partition  function 
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k  ■  Boltzmann's  constant 

h  ■  Planck's  constant 

The  subscripts,  a  and  b,  denote  the  equilibrium  states  ahead  of 
and  behind  the  shock  wave,  respectively.  For  the  incident  shock, 
the  values  of  a  and  b  are  1  and  2,  respectively,  and  for  the  re¬ 
flected  shock  a  and  b  are  2  and  5,  respectively. 

A  computer  program  was  developed  to  solve  these  equations 
iteratively , for  both  the  incident  shock  and  the  reflected  shock 
regions.  The  results  of  these  solutions  were  used  to  model  the 
nozzle  expansions  and  the  thermochemical  state  of  the  expanded 
helium/fluorine  mixtures.  Figure  13  shows  the  fluorine  dissociation 
level  in  the  reflected  shock  region  for  different  initial  driven 
tube  pressures  as  a  function  of  initial  shock  Mach  number.  From 
these  results,  it  is  clear  that  the  dissociation  level  is  higher 
for  lower  initial  pressure  at  any  given  value  of  initial  shock 
Mach  number.  Table  II  summarizes  the  nozzle  dimensions  used  in 
this  study. 


Table  II  SUMMARY  OF  AREA  RATIOS 


Nozzle 

Throat 

Dimension,  in. 

Exit 

Dimensions,  in. 

A/A* 

Injection 

Facility 

0.264  high  by 
0.500  wide 

3.439  high  by 
0.500  wide 

13.10 

Parallel 

Mixing 

Facility 

H,  nozzle 

0.050  high  by 
5.00  wide 

nozzle 

0.794  high  by 
5.00  wide 

15.88 

F2  nozzle 

0.062  high  by 
5.00  wide 

F?  nozzle 

0.652  high  by 
5.00  wide 

10.51 

Note:  "A"  denotes  the  area  of  the  nozzle,  exit  plane 
and  "A*"  denotes  the  area  of  the  nozzle 
throat. 
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Figure  13.  Fluorine  Dissociation  Fraction  vs  Incident  Shock  Mach  Number 


Figure  14  shows  the  computed  nozzle  exit  pressures  as  a  func¬ 
tion  of  initial  driven  tube  conditions  and  shock  strength  for  a 
nozzle  with  A/A*  =  13.1.  Again,  for  a  fixed  initial  shock  Mach 
number,  the  nozzle  exit  pressure  increases  with  pi .  The  results 
given  in  Figure  14  are  for  the  injection  mixing  nozzle  and  a  SOX 
by  volume  test  gas  mixture  of  He  and  F^.  Figure  15  shows  typical 
nozzle  exit  static  pressures  for  the  parallel  mixing  hydrogen 
nozzle  (A/A*  =  15.88).  A  complete  set  of  these  data  and  several 
cross  plots  were  used  to  design  the  run  conditions  for  this  study. 

The  nozzle  expansion  processes  were  modeled  using  a  similar 
set  of  conservation  equations  and  thermodynamic  state  equations 
as  used  for  the  reflected  shock  analysis,  together  with  suitable 
modifications  for  the  boundary  layer  formation  along  the  nozzle 
walls.  The  procedure  discussed  here  was  applied  to  the  perfor¬ 
mance  analysis  of  both  the  injection  shock  tube  nozzles  and  the 
parallel  mixing  shock  tube  nozzles  whose  designs  and  geometries 
have  already  been  discussed.  In  both  cases,  the  reflected  shock 
region  represents  a  stagnation  reservoir  and  the  real-gas  expansion 
processes  were  carried  out  to  determine  the  thermochemical  changes 
at  the  nozzle  exit  planes.  The  conservation  of  mass  eauation  for 
a  converging/ diverging  nozzle  with  adiabatic  inviscid  flow  is 
used,  along  with  conservation  of  momentum  and  energy  and  the  thermal 
and  caloric  equations  of  state.  Again,  the  equations  were  solved 
by  iteration  to  give  exit  plane  values  of  the  flow  velocity  and 
the  thermodynamic  properties  of  the  gas  mixtures,  as  a  function 
of  the  nozzle  area  expansion  ratio. 


The  governing  equations  for  the  nozzle  flow  field  are: 

Conservation  of  mass 

r.uA  =  C]  =  constant 

Conservation  of  energy 

h  +  h  u^  =  C<_  =  constant 

Conservation  of  momentum 

du  ,  1  dP  n 
U  dx  p  dx  ' 


(16) 


(17) 


(18) 


Thermal  equation  of  state 
P  =  (1  +  a)  oRT 


(19) 
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Test  t.as  -  +  He ,  50 L  by  volume 

p  =  Initial  Driver  Section  Pressure 


•on  qoBK  q^oqs  iuappuj 
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Figure  14.  Nozzle  Exit  Static  Pressure  vs  Incident  Shock  Mach  Number 


6.0 
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3.5 


Caloric  equation  of  state 
h  =  (4  +  <)  T  + 


(20) 


The  law  of  mass  action  is  given  by: 


D 


where  a  is  taken  as  constant  at  the  value  computed  for  reflected 
shock  conditions.  This  system  of  equations  was  solved  by  iter¬ 
ation  to  give  thermodynamic  conditions  at  the  nozzle  exit  planes. 

In  all  cases,  the  experimentally  measured  static  pressures  at 
the  nozzle  exits  are  higher  than  the  theoretically  predicted  exit 
plane  pressures.  This  discrepancy  was  not  discovered  until  late 
in  the  experimental  program  and,  hence,  could  not  be  corrected  in 
time  to  rerun  the  series  of  interferometric  tests.  Two  possible 
causes  are  being  investigated  to  explain  this  difference.  First, 
nozzle  boundary  layers  are  being  reexamined  theoretically  in  more 
detail  to  determine  if  they  are  really  sufficiently  large  to 
chanee  the  effective  nozzle  area  ratio.  Second,  reflected  shock/ 
boundary  layer  interactions  unique  to  reactive  diatomic  gases 
that  appear  to  result  in  higher  reflected  shock  pressures  are  being 
examined , 


The  initial  boundary  layer  thickness  estimates  at  the  nozzle 
exit  planes  were  based  on  conventional  simplified  techniques  (Ref. 
6)  that  use  the  local  flow-field  Reynolds  number  and  the  length 
of  the  nozzle  wall  as  the  characteristic  dimension,  X.  The  exit 
plane  Reynolds  number  can  be  expressed  as: 


Re  = 


u 

e  e  X 
'e 


(22) 


where  the  subscript  e  denotes  properties  evaluated  at  the  nozzle 
exit  plane.  The  laminar  boundary  layer  thickness  is  then  given 

by 


which  can  be  rewritten  as 
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These  equations  give  typical  boundary  layer  thicknesses  of 
0.05  in.  for  the  parallel  mixing  hydrogen  carrying  nozzle  and 
0.03  in.  for  the  fluorine  carrying  nozzle.  Based  on  these  num¬ 
bers,  it  was  originally  concluded  that  boundary  layers  would  not 
alter  the  effective  nozzle  area  ratio  in  a  significant  manner. 

It  is  possible  that  a  turbulent  boundary  layer  is,  in  fact,  more 
nearly  correct  for  the  flow  condition  here.  The  criterion  for 
determining  the  nature  of  the  wall  boundary  layer  is  the  relative 
magnitude  of  the  momentum  thickness  Reynolds  number  Re  .  That 

is,  if  Re^  (based  on  freestream  conditions)  is  less  than  500, 

the  boundary  layer  can  be  assumed  laminar,  while  values  greater 
than  500  indicate  the  presence  of  a  turbulent  boundary  layer. 

These  calculations  and  consideration  of  increased  reflected  shock 
pressures  are  in  progress. 

4.  DIAGNOSTIC  INSTRUMENTATION 

Three  basically  different  experimental  techniques  were  employed 
in  this  program  to  measure  the  variations  in  flow-field  properties 
within  the  mixing  region  downstream  of  the  nozzle  exit.  The  three 
different  techniques  include  a  Mach-Zehnder  interferometer  (MZI) 
designed  to  map  the  mixing  layer,  a  rapid  scanning  spectro¬ 
meter  to  map  the  infrared  emission  intensity  within  the  mix- 
.  lg  layer,  and  a  series  of  probe  measurements  tailored  to  map¬ 
ping  the  temperature  and  static  and  pitot  pressures  at  the  nozzle 
exit  region  and  in  the  mixing  layers.  The  actual  instrumentation 
and  the  basic  techniques  are  interchangable  for  both  the  injection 
mixing  and  the  parallel  shear  flow  experiments  with  only  minor 
variations.  The  remainder  of  this  subsection  will  be  devoted  to 
a  detailed  discussion  of  these  diagnostic  techniques  and  the  gen¬ 
eralized  data  reduction  procedure  that  has  been  used  in  this  study. 

a.  Mach-Zehnder  Interferometer  -  All  optical  examinations  of  com¬ 
pressible  fluid  flows  are  based  on  the  fundamental  premise  that 
variations  in  density  across  the  transluminated  medium  give  rise 
to  variations  in  the  absolute  local  index  of  refraction  that,  in 
turn,  affect  the  traversing  light  rays.  Appropriate  optical 
equipment  is  then  used  to  translate  the  resulting  effect  into 
variations  of  intensity  of  illumination  on  a  screen  or  photographic 
plate.  The  basic  physical  laws  used  in  interpreting  these  effects 
are  Fermat's  law  of  stationary  transit  time  and  Snell's  law  of 
refraction  at  surfaces  of  discontinuity  of  refractive  index, 
n  (x,y , z) . 
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The  principal  feature  of  all  of  the  well-known  optical  sys¬ 
tems  can  best  be  described  in  the  following  manner.  Consider  a 
disturbance  in  a  flow  field  that  is  illuminated  by  a  parallel 
beam  of  light  that  can  be  observed  on  a  screen  (Fig.  16). 


P*,  t* 


Figure  16.  Generalized  Ray  Diagram 

Consider  next  what  happens  to  a  single  ray,  like  R,  as  it  trans- 
verses  the  distuibance.  In  the  absence  of  any  disturbance,  the 
ray  would  arrive  at  the  point  P  on  the  screen  at  a  time  t.  How¬ 
ever,  the  presence  of  the  disturbance  alters  the  path  of  the  ray 
R,  and  it  actually  arrives  at  the  screen  at  point  P*  at  time  t* 
and  from  a  direction  characterized  by  6. 

The  insertion  of  appropriate  optical  equipment  into  the  ex¬ 
isting  light  path  gives  a  record  of  either  the  phase  lag  (t*-  t), 
the  deflection  (0),  the  displacement  (P*-P),  or,  in  general,  some 
combination  of  all  three  of  these  quantities  at  the  position  of 
the  screen.  The  phase  lag  determination  is  the  physical  basis 
for  interferometry,  while  measuring  the  deflection  angle  forms  the 
basis  of  schlieren  methods  and  the  displacement  is  measured  by  the 
familiar  shadowgraph  technique.  The  theory  of  each  of  these  sys¬ 
tems  gives  a  unique  correspondence  between  the  intensity  of  illu¬ 
mination  at  the  screen  and  n  (x,y,z).  Once  n  (x,y,z)  is  known, 
an  analysis  b  sed  on  electromagnetic  theory  gives  the  gas  intensity 
in  terms  of  n  and  a  constant,  K  (Gladstone-Dale) .  This  relationship 
is 


n-1  =  K.  (25) 

where  K  is  a  function  of  the  molecular  weight  and  the  polarizability 
per  atomic  particle.  In  general,  the  polarizability  is  a  function 
of  wavelength  but  this  dependency  is  very  weak  and  can  usually  be 
ignored. 
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In  optical  systems  where  either  the  deflection  (0)  or  the 
phase  lag  (t*-t)  is  to  be  measured,  the  displacement  (P*- P)  must 
be  eliminated.  To  accomplish  this,  an  objective  or  field  lens  is 
used  to  focus  a  particular  plane  of  the  test  section  on  the  ob¬ 
servation  screen. 

As  mentioned,  the  basic  quantity  measured  in  interferometry 
is  the  phase  lag  associated  with  a  reference  beam  passing  through 
an  undisturbed  gas  as  compared  to  its  traversing  a  disturbance 
under  investigation.  To  understand  how  the  interference  technique 
works,  consider  a  Mach-Zehnder  interferometer  (MZI)  shown  in  Fig¬ 
ure  17  where  one  assumes  that  all  of  the  optical  components  are 
perfect.  First,  note  that  the  division  of  amplitude  at  the  beam 
splitter  produces  two  coherent  waves.  Now  if  the  traversal  time 
for  the  two  waves  to  reach  the  second  beam  splitter  via  two  dif¬ 
ferent  paths  is  exactly  the  same,  the  two  waves  will  reconstruct 
the  original  wave  front.  Otherwise,  we  may  expect  to  obtain  some 
interference  phenomena  that  will  be  evidenced  by  a  pattern  of 
fringes.  The  possibility  of  this  interference  occurring  can  be 
explained  in  any  of  three  ways:  the  four  optical  plates  may  not 
be  exactly  parallel  to  each  other;  the  linear  distances  along  the 
two  paths  may  not  be  identical;  and  a  change  in  refractivity  due 
to  a  disturbance  in  the  test  section  may  change  the  phase  velocity 
of  the  test  beam  relative  to  the  reference  beam. 

The  ray  diagram  shown  in  Figure  17  could  have  been  drawn  for 
the  case  where  two  distinct  sets  of  fringes  are  visible  (broken 
lines).  Generally  this  'lost':  set  of  fringes  is  not  used,  as 
the  quality  (a  result  of  imperfect  beam  splitters)  is  not  as 
good  as  that  of  the  fringes  formed  by  the  primary  beam.  As  will 
be  discussed  in  Section  III,  Martin  Marietta  used  the  secondary 
fringe  pattern  to  observe  the  total  integrated  radiation  intensity 
from  the  mixing  region. 

If  one  enclosed  the  entire  instrument  with  a  box  and  viewed  the 
exit  beam  only,  there  would  be  no  way  to  describe  what  optical 
elements  were  inside  the  box.  Doing  this,  one  could  actually  see 
two  apparent  sources  and  the  light  rays  intersecting  at  some  plane 
within  the  enclosure.  This  plane  then  is  the  location  at  which 
the  fringes  appear  to  be  formed  but,  of  course,  there  can  actually 
be  no  interference  until  the  two  beams  are  reunited  at  the  last 
beam  splitter.  This  is  the  old  idea  of  forming  a  real  image  (at 
the  film  plane)  of  a  virtual  object  (the  apparent  intersection  of 
the  test  and  reference  beams).  The  main  utility  of  the  MZI  is  the 
wide  range  of  choice  one  has  in  locating  the  plane  of  intersection 
(i.e.,  the  focal  plane  within  the  test  section)  and  adjusting  the 
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Lost  Beam 
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Figure  17.  Schematic  of  Basic  M2 I  Arrangement 


width  and  orientation  of  the  fringes.  The  fringe  width,  W,  from 
linear  theory  is  related  to  the  apparent  angle  of  intersection,  if, 
of  the  two  beams  by 

W  »  - 

where  X  is  the  wavelength  of  the  light  source.  Adjustment  of 
either  of  the  beam  splitters  allows  control  over  both  the  final 
fringe  width  and  fringe  orientation. 

The  actual  construction  of  the  MZ1  used  in  this  study  is  by 
now  considered  as  conventional  and,  hence,  will  not  be  discussed 
in  great  detail.  There  are  a  number  of  excellent  descriptions 
in  the  literature  of  the  problems  of  design,  construction,  and 
alignment  adjustments  of  this  type  of  interferometer  (Ref.  7 
thru  9).  Figure  18  is  an  isometric  sketch  of  the  associated  test 
hardware  and  the  MZI  mounted  with  its  plane  of  centers  perpendic¬ 
ular  to  the  axis  of  the  shock  tube.  The  lower  beam  passes  through 
two  windows  mounted  at  the  nozzle  exit  region  of  the  test  section, 
and  the  upper  beam  passes  through  a  matched  pair  of  windows  in 
the  compensating  beam.  Both  the  test  section  windows  and  the 
compensating  windows  are  made  of  optically  flat  (A/ 1C)  schlieren 
quality  quartz  with  faces  parallel  to  within  0.0005  in.  Figure 
18  also  shows  the  concept  of  a  radiation  stop  that  is  employed 
on  the  camera  side  of  the  MZI.  The  diffuse  burst  of  radiation 
from  the  mixing  layer  is  sufficiently  intense  to  overexpose  the 
interferogram  recording  film  unless  some  effective  radiation  stop 
is  employed.  As  illustrated  in  Figure  18,  the  field  lens  is 
some  distance  from  the  test  section,  and  the  collimated  beam 
from  the  light  source  is  completely  transmitted  through  the 
pin  hole  located  at  the  focal  point  of  the  lens  while  the  diffuse 
radiation  from  the  mixing  layer  is  blocked  at  the  focal  plane 
of  the  field  lens.  This  technique  has  been  used  in  this  study 
and  the  resulting  interf erograms  are  not  influenced  by  the  radia¬ 
tion  resulting  from  either  the  mixing  region  or  the  driver  gases 
arriving  at  the  end  of  the  uniform  test  flow. 
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Figure  18,  Schematic  of  MZI  Installation  Showing  Shock 
Tube  and  Radiation  Stop 


As  previously  mentioned,  a  great  many  practical  details  such 
as  alignment,  light  source  size,  imperfections  in  the  optical 
components,  and  proper  choice  of  fringe  width  and  orientation 
must  be  considered  in  obtaining  high-quality  interferograms.  With¬ 
out  intending  to  underestimate  the  importance  of  these  practical 
aspects,  which  are  generally  a  matter  of  experience  and  practice, 
a  generalized  discussion  of  the  evaluation  of  an  interferogram  is 
presented. 

The  concept  fundamental  to  the  interferometric  determination 
of  the  density  distribution  in  gas  flows  is  that  if  a  change  in 
the  test  section  density  by  an  amount  Ap  occurs,  then  the  resulting 
phase  change  produces  a  fringe  shift,  5  (nondimensional  number), 
given  by 


(n-i) : 
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where  K  is  the  Gladstone-Dale  constant  (-0.22  cm;/gn  for  air)  and 
L  is  the  geometrical  path  length  over  which  the  change  occurs. 
Knowing  K,  L,  and  '•  enables  one  to  calculate  ip  from  a  measurement 
of  c.  The  above  equation  is  written  for  a  change  in  optical  path 
(or  density)  which  is  the  same  for  all  light  rays  and,  hence, 
results  in  a  uniform  translation  without  distortion  of  the  entire 
fringe  pattern.  In  a  more  realistic  case,  we  must  assume  that 
the  test  section  density  is  a  function  of  the  spatial  coordinates. 
Thus,  in  a  two-dimensional  flow  field,  the  optical  path  of  the 
light  rays  will  depend  on  the  particular  (x,y)  coordinates  at 
which  the  Z  traversal  is  made  (see  Fig.  19). 


Y 
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Figure  19.  Generalized  Light  Ray  Path 


The  optical  path  is  defined  by  the  following  integral  equation 
over  the  path,  s 

L 


optical  path 


1  " 
•a 


(x,y,z)  ds 


(28) 


where  ds  is  an  element  of  distance  along  the  actual  path  a  light 
ray  travels.  Consider  now  a  situation  where  the  refractive  index 
is  uniform  in  the  reference  leg  and  a  function  of  (x,y,z)  in  the 
test  section.  Then,  if  there  is  some  portion  of  the  flow  field 
in  the  test  section  that  is  known  to  be  uniform  (say  xo>yo>zo), 

the  fringe  shift  between  the  nonuniform  portion  of  the  field  and 
the  uniform  (reference  position  in  this  example)  part  is  given  by: 


(x,y) 


1 


L 


n  (x,y,z)  -  n 


ds 


(29) 
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The  important  point  to  notice  here  is  that  the  reference  beam 
need  not  actually  be  in  a  leg  of  the  instrument  that  is  exterior 
to  the  test  section.  In  fact,  for  the  example  situation  described 
above,  one  does  not  require  a  "no-flow"  interferogram  in  addition 
to  the  flow  interferogram.  That  is,  the  change  in  density  is 
measured  relative  to  some  portion  of  the  flow  field  where  the 
density  is  already  known  and  the  fringe  shift  of  interest  is  shown 
on  a  single  flow-field  interferogram.  In  this  study,  a  no-flow 
interferogram  was  used  in  conjunction  with  the  flow  interferogram 
to  establish  the  overall  density  in  the  undisturbed  nozzle  exit 
region  before  mixing. 

Returning  now  to  the  generalized  integral  equation  for  5  (x,y), 
if  the  path  deviation  of  the  actual  light  ray  from  a  straight 
line  is  not  large,  one  can  approximate  ds  (actual  path  segment)  by 
dz  (ideal  path  segment).  Doing  this,  and  using  the  relationship 
between  n  and  .  ,  the  generalized  expression  for  fringe  shift  be¬ 
comes 


(x.y)  =  y 

r  L 

/ 

L 

(x,y,z)  dz  -  I  ( 

(x  ,y  ,z  )  dz 
\  o,Jo'  0 !  0 

L  0 

r 

0 

_ 

(30) 


Further  simplification  of  this  result  is  possible  if  one  assumes 
that  the  flow-field  variables  are  two-dimensional  and,  as  such, 
depend  only  on  (x,y).  In  this  case,  the  basic  evaluation  equation 
becomes 


—  ;  (x,y)  =  r  (x,y)  -  cq  (31) 

One  useful  scheme  for  gaining  confidence  in  the  assumption 
that  the  deviation  of  a  light  ray  is  small  (i.e.,  ds  -  dz)  is  to 
very  carefully  focus  on  a  grid  network  (most  often  etched  on  the 
exit  window  of  the  test  section)  and  photograph  it  under  flow 
and  no-flow  conditijns.  If  one  cannot  detect  a  distortion  of  this 
grid  work  under  actual  flow  conditions,  then  the  assumption  of 
approximating  ds  by  dz  is  for  all  practical  purposes  a  valid  one. 

In  this  same  regard,  when  the  interferometer  is  used  for  examining 
regions,  such  as  boundary  layers,  where  there  are  large  density 
gradients,  it  is  necessary  to  consider  the  possibility  of  an  error 
in  the  measured  fringe  shift  that  occurs  as  a  result  of  refraction 
of  the  light  rays.  That  is,  the  gradient  in  density  will  cause 
a  light  ray  to  bend  in  the  direction  of  increasing  density  and, 
hence,  follow  a  curved  path  through  the  test  section  rather  than 
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the  straight  path  that  is  assumed  in  the  above  theoretical  dis¬ 
cussion.  This  phenomenon  can  result  in  two  main  effects:  First, 
each  ray  in  the  strong  gradient  region  will  travel  over  a  longer 
optical  path;  second,  along  this  path,  the  refractive  index  (and, 
hence,  density)  is  constantly  changing.  Several  authors  have 
addressed  this  problem,  including  Howes  and  Buchele  (Ref.  10) 
and  Bunting  and  DeVoto  (Ref.  11).  The  main  conclusion  drawn  by 
each  of  these  authors  is  that  the  fringes  should  be  focused  at 
the  2/3  span  of  the  test  section  rather  than  center  span  to  min¬ 
imize  the  refractive  error.  This  study  always  used  the  2/3  span 
location  for  focusing  the  fringes  even  though  the  gradients  here 
are  not  nearly  as  large  as  those  in,  for  example,  a  thermal 
boundary  layer. 

The  above  discussion  has  been  a  rather  generalized  treatment 
of  the  basic  principles  of  interferometry.  A  more  specialized 
treatment  starting  from  these  general  results  is  presented  in 
the  following  paragraphs.  This  division  of  discussion  is  necessary 
because  several  bench  experiments  were  required  in  this  study 
to  generate  refractive  index  variations  for  F2  and  HF,  which 
had  not  previously  been  reported  in  the  literature.  Thus,  the 
details  of  measuring  refractive  indices  of  gases  and  the  appli¬ 
cation  of  interferometric  techniques  to  multicomponent  gases  will 
now  be  discussed  starting  with  the  results  of  the  above  generalized 
discussion. 

As  shown  above,  the  absolute  fringe  shift,  6,  for  monochromatic 
illumination  at  wavelength,  X,  for  a  two-dimensional  flow  field  is 
given  as 


where  n  is  th*  refractive  index  for  the  test  case  under  consider¬ 
ation,  n^  is  the  refractive  index  for  some  initially  known  condi¬ 
tions,  and  L  is  the  physical  pathlength  of  the  test  cell.  The 
reference  conditions  for  measuring  the  refractive  indices  of  HF  and 

F  are  an  initially  evacuated  test  cell  where  n  can  be  taken  as 

o 

unity.  Hence,  using  this  criterion  as  a  basis,  the  expression 
for  can  be  written  as 


;  „  (n-1)  L 

) 


(33) 
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Recalling  now  the  relationships  between  refractive  index  and 
density,  for  the  present  situation  using  an  initially  evacuated 
test  cell  as  reference,  the  absolute  fringe  shift  can  be  writ¬ 
ten  as 


(n-1)  L  Kt'L 

,  S  ■■  "  "  ■  — -  =  • . — 1 

X  X 


(34) 


Because  all  of  the  refractive  index  data  to  be  used  in  this 
study  are  referenced  to  standard  conditions  of  temperature  and 
pressure,  the  above  relationship  can  be  written  as 


'  STP 


(35) 


where  the  ratio 

(n-1)  =  Kc _ 

|nSTp-l\  K.  STp 


(36) 


has  been  used  along  with  the  fact  that  K  remains  only  a  slightly 
varying  function  of  X. 

.Vow  the  retractive  index  of  any  given  chemical  species  i  can 
be  expressed  i»-  terms  of  the  usual  two-term  Cauchy  formula 

B. 

(n-1) .  =  A .  | ?  v  j  (37) 

X 

where  and  B .  are  uniquely  determined  for  each  of  the  chemical 

species  involved.  These  constants  have  been  experimentally  de¬ 
termined  for  f  and  HF,  and  these  data  will  be  presented  later 

in  this  section. 

Continuing  now  with  the  above  expression  for  the  fringe  shift 
of  a  single  .omponent  gas,  one  can  write 

(Vr  <w 

w'ne  re 

( nSTP" ' ) 
ci  “  , 

’  STP 


40 


and  Xj  denotes  the  different  wavelengths  at  which  interferograms 
were  recorded. 


If  the  test  section  contains  several  chemical  species,  then 
their  respective  concentrations  at  any  given  point  in  the  flow 
field  can  be  deduced  from  interferograms  by  noting  that  the 
total  composite  fringe  shift,  6  at  a  point  is  just  the 

superposition  of  the  fringe  shifts,  6^  N  ,  due  to  each  of 

the  chemical  species  i.  Here  one  should  note  that  a  different 
wavelength  interferogram  is  required  for  each  different  chemical 
species  to  be  determined.  Expanding  the  above  single  species 
equation  using  self-explanatory  notation  gives 

5  N  ■  l  si  (h) 

or 

6  I*1)  =  C1HF  PHF  +  C1F2 

5  (A2)  =  C2HF  PHF  +  C2F2 

5  (X3)  =  C3HF  PHF  +  C3F2 
where 

c..  -  ("STP~  1)t  L 

J  '  1  PSTP 


pf2  +  cih2  ph2 

PF2  +  C2H2  PH2 
PF2  +  C3H2  PH2 


which  are  the  three  simultaneous  equations  that  can  be  solved 

for 


and  P„  once  the  C..  coefficients  are  known,  and 
H2  ji 


the 


"HF ’  "F2  ’ 

6  are  measured. 


Before  proceeding  further  with  this  discussion,  the  results 
of  bench  experiments  conducted  on  a  separate  IRAD  task  (Ref.  12) 
to  obtain  the  A  and  B  coefficients  in  the  Cauchy  representation 
of  the  refractive  index  of  F2  and  HF  will  be  discussed. 


The  refractive  indices  of  argon,  fluorine,  and  hydrogen  flu¬ 
oride  have  been  measured  directly  and,  along  with  values  quoted 
by  Allen  (Ref.  13),  are  plotted  in  Figure  20.  The  method  used 
here  is  based  on  the  fringe  shift  observed  with  an  MZI  as  a 
small  test  cell  is  slowly  filled  with  the  gas  to  be  observed. 
Figure  21  shows  a  diagrammatic  sketch  of  the  test  setup. 


As  has  been  shown  in  this  section,  the  fringe  shift,  <5,  observed 
as  the  chamber  pressure  is  changed  from  a  value  where  density  is 
denoted  by  to  a  final  pressure  with  density  corresponding  to 

.  can  be  expressed  as 


nSTP~  1 
'  STP 


(43) 


For  the  bench  condition  established  for  this  purpose,  is  iden¬ 
tically  zero,  so  the  observed  fringe  shift  will  be  given  by 

.  HsTP-1  L 

•'  STP  2 


(44) 


or  expressed  in  terms  of  experimental  observables  (assuming  a  per¬ 
fect  gas) 


nSTP-1 


PSTP 

T 

STP 


L 


(45) 


The  technique  employed  here  was  to  measure  ;  as  the  test  cell 
pressure  "as  varied  very  slowly.  Increasing  the  test  cell  pressure 
rapidly  resulted  in  a  nonisothermal  expansion  of  the  gas  and  thus 
a  changing  temperature.  In  the  experiments  with  HF  as  a  test  gas, 
it  was  necessary  to  monitor  the  test  cell  temperature  in  addition 
to  the  pressure  because  the  windows  in  the  cell  were  destroyed  if 
the  pressure  was  allowed  to  rise  sufficiently  slowly  to  maintain 
an  isothermal  expansion.  Fringe-shift  data  for  each  of  the  three 
test  gaseu  considered  (A,  F  ,  HF)  were  recorded  as  a  function  of 
temperature  and  pressure  at  three  separate  wavelengths.  Data  points 
for  pressure  and  fringe  shift  were  plotted  for  each  different  wave¬ 
length  (see  Fig.  22),  and  the  experimental  value  of  n^p-1  was  de¬ 
termined  from  the  measured  slope  using  Equation  45.  Because  of 
the  uncertainty  in  both  pressure  and  temperature,  an  error  bar  as 
shown  in  Figure  20  was  assigned  to  each  of  these  values.  Doing 
this  for  each  of  the  three  different  observed  wavelengths  allows 
one  to  produce  plots  of  (ngTp"l)  vs  k  such  as  those  shown  in  Figure 

20.  This  variation  of  n-1  with  *  is  then  used  to  produce  the  A 
and  11  constants  in  the  Cauchy  formula,  Equation  37,  with  the  results 
as  shown  in  Table  III. 
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Figure  22.  Fringe  Shift  vs  Pressure  for  Refractive  Index  Measurement  of  HF 


Table  III 


EXPERIMENTAL  VALUES  OF  A  AND  B  FOR  VARIOUS  SPECIES 


Species 

A 

B, 

HF 

12.62  x  10“5 

16.40  x  10-3 

f2 

19.42  x  10“5 

19.24  x  10"3 

A 

27.92  x  10"5 

5.60  x  10-3 

H, 

13.58  x  10-8 

7.52  x  10-3 

Returning  now  to  the  set  of  three  simultaneous  equations  for 
solving  for  .  „  ,  r>  and  cu  ,  one  can  use  the  results  of  Table 

f-i  hr  K' 

<■  *■ 

III  to  compute  the  C  coetficients  in  these  equations.  Using  the 

values  from  Table  III,  the  actual  test  section  span  of  5  in.,  and 
Loschmidts  number  for  "g^p,  the  coefficients  corresponding  to 

the  three  wavelengths  used  for  recording  the  interferometric  re¬ 
sults  were  computed  and  are  given  in  Table  IV. 


Table  IV 


CO'IPUTED  VALUES  OF  COEFFICIENTS 


j 

0 

X ,  A 

Cj HF’  Cm3 

CjF:»  cnr 

c jh2 ’  cm3 

1 

4358 

1.487  x  10"16 

2.319  x  10~13 

1.531  x  10"18 

2 

4880 

1.306  x  10" 18 

2.033  x  10"1B 

1.357  x  10~18 

3 

6100 

1.021  x  10" :  " 

1.583  x  10-18 

1.073  x  10~18 

By  inserting  each  of  these  nine  C  coefficients  into  the  pre¬ 
viously  developed  set  of  simultaneous  equations,  one  can,  in  prin¬ 
ciple,  obtain  absolute  number  densities  through  the  mixing  layer 
by  measuring  corresponding  fringe  shifts  at  three  different  wave¬ 
lengths.  As  it  turns  out,  this  ambitious  task  is  not  possible 
because  the  solution  set  to  these  equations  is  extremely  sensitive 
to  the  precision  to  which  the  resulting  fringe  shifts  can  be  meas¬ 
ured.  There  are  two  different  ways  of  seeing  that  this  is  indeed 
the  case  in  the  present  situation.  First,  one  can  consider  each 
of  the  elements,  C^,  in  a  row  of  the  matrix  as  being  a  component 

of  a  three-dimensional  vector.  Doing  this,  normalizing  the  vectors, 
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and  taking  the  dot  product  of  each  vector  with  the  others  results 
in  dot  products  extremely  close  to  unity.  This  means  that  the 
vectors  are  almost  colinear  and,  hence,  finding  a  solution  set 
involves  finding  the  difference  between  very  small  vectors  that 
consist  of  the  closure  vector  between  any  two  pairs.  Another  way 
of  looking  at  this  overall  problem  is  to  note  that  the  three  plots 
of  nSTp-l  vs  X  are  very  nearly  parallel, even  though  displaced  in 

absolute  value.  After  a  careful  analysis  of  the  values  involved, 
it  turns  out  that  one  needs  to  know  the  absolute  fringe  shifts 
for  each  of  the  different  wavelengths  to  better  than  0.001  fringe. 
This  type  of  precision  is  beyond  the  present  capability  and,  hence, 
the  three  wavelength  solutions  originally  planned  for  determining 
absolute  number  densities  are  not  possible.  These  number  densities 
are  obtained  by  another  technique  involving  certain  additional  as¬ 
sumptions  that  will  be  explained  in  Section  III. 

The  scheme  used  to  reduce  the  interferometric  data  involved 
a  Scherr  Tumico  optical  comparator.  This  comparator  has  a  very 
fine  X-Y  translational  table  and  an  enlarging  viewer  with  cross¬ 
hairs.  A  negative  of  the  interferogram  is  placed  on  the  trans¬ 
lational  table,  and  the  enlarged  fringes  plus  crosshairs  are 
viewed  on  a  screen.  An  X-Y  scan  is  performed  on  both  the  flow 
and  no-flow  negatives  with  corresponding  fringe  maxima  and  minima 
recorded.  These  raw  data  are  then  reduced  in  terms  of  absolute 
fringe  shift  (i.e.,  fringe  movement  between  no-flow  and  flow 
negatives  divided  by  original  fringe  spacing)  at  preselected  down¬ 
stream  stations  measured  from  the  edge  of  the  splitter  plate.  Re¬ 
sults  of  these  efforts  and  a  summary  of  the  data  obtained  will  be 
presented  in  Section  III. 

b.  Rapid  Scanning  Spectrometer  -  An  OCLI  Model  501  rapid  scanning 
spectrometer  was  used  to  monitor  the  absolute  emission  intensity 
between  1.7  and  4.7  pm  within  the  mixing  zone  of  both  the  H2 
injection  and  the  parallel  shear  flow  experiments.  Figure  23 
is  a  schematic  representation  of  the  optical  system.  Figure  24 
is  a  photograph  of  the  spectrometer  located  at  the  shock  tube  test 
section.  Translation  supports  mounted  to  the  basic  instrument 
cover  provide  a  means  for  varying  the  focal  point  location  to  any 
given  station  within  the  mixing  layer. 
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Figure  23.  Rapid  Scanning  Spectrometer  Optical  Schematic 


1 


The  Model  501  collects  radiation  fror.  any  given  zone  of  the 
reaction  layer  through  a  variable  focus  Cassegrainian  optical 
system.  This  radiation  then  passes  through  a  double  pass  grating 
monochromator  optical  system  onto  a  double-exit  slit  and  detector 
configuration.  A  scan  wheel  equipped  with  corner  mirrors  causes 
the  second  pass  radiation  to  sweep  in  wavelength  across  the  fixed 
exit  slits.  The  shorter  wavelength  half  of  the  spectrum  is  scan¬ 
ned  across  one  exit  slit  at  the  same  time  that  the  longer  wave¬ 
length  portion  is  scanned  across  the  other  exit  slit.  The  avail¬ 
ability  of  two  exit  slits  allows  coverage  over  a  broader  spectral 
region  with  higher  resolution  and  provides  a  wavelength  overlap 
useful  in  checking  system  accuracy.  The  spectral  region  to  be 
scanned  is  determined  by  the  gratings,  filters,  and  detectors 
used.  The  detectors  used  in  this  study  are  indium  arsenide  (InAs) 
and  indium  antimonide  (InSb)  cooled  to  LN  temperature.  The 
signal  resulting  from  the  detected  radiation  is  displayed  on  two 
separate  oscilloscopes  and  recorded  photographically.  The  Model 
501  has  a  built-in  secondary  standard  of  radiation  that  can  be 
used,  when  fully  calibrated,  to  convert  observed  radiation  signals 
to  units  of  absolute  radiance.  This  instrument  is  capable  of  a 
full  spectral  scan  between  1.7  and  4.7  urn  in  1  msec,  with 
a  repetition  time  of  1.25  msec. 

The  data  reduction  procedure  used  here  for  deducing  the  spe¬ 
cific  gas  intensity  from  measured  spectrometer  signals  is  similar 
to  that  outlined  by  Menard  and  Thomas  (Ref.  14).  The  portion  of 
their  analysis  that  follows  serves  as  an  introduction  to  the 
specific  details  of  the  present  calibration. 


Assuming  an  optically  thin  gas  and  an  isotropic  radiation 
source,  the  radiant  power  (w/jm)  to  the  detector,  at  a  specific 
wavelength  ■,  is  given  as 


P  = 
g 


dAdz 


and  the  power  from  the  standard  calibration  source  is  given  as 


P  = 
c 


(46) 


(47) 


where  £,  is  a  combined  parameter  containing  all  transmission  losses, 
monochromator  efficiency,  and  detector  sensitivity,  G.  is  the  spe¬ 
cific  radiation  intensity  (w/'cn^  -jm  -sr)  of  the  gas,  '  is  the 
solid  angle  where  A  is  the  cross-sectional  area  of  the  field  of 
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view  rays  as  defined  by  the  entrance  optics,  and  z  is  the  coor¬ 
dinate  along  the  optical  path.  N  ,  the  intensity  (w/cm‘  -urn  -sr) 

of  the  blackbody  calibration  source,  is  given  in  terms  of  the 
absolute  temperature  of  the  cavity  by  Planck's  radiation  law, 


Cj  and  Cz  are  the  first  and  second  radiation  constants  given  by 
Stair,  Johnston,  and  Halbach  (Ref.  15)  as 

Cj  =  1.19088  x  10- 12  w-cm2/sr 

C ?  =  1.4380  cm-°K 

Over  the  linear  response  range  of  the  photo  detectors,  the 
output  voltage,  E,  is  proportional  to  the  input  power,  P.  There¬ 
fore,  by  taking  the  ratio  of  Equations  (46)  and  (47),  the  specific 
intensity  emitted  by  the  gas  at  wavelength  A  is  given  as 

E 

G>.  *  '“a  /  (49) 

c 


where  K,  the  view  factor  coefficient,  is  defined  by 
K  „  dA 

JJJd.  dAdz 


An  analytical  solution  of  Equation  (50)  would  be  very  dif¬ 
ficult  for  the  particular  geometry  used  in  this  experiment,  so 
an  alternative  method  is  employed  to  determine  K  experimentally. 
Integration  of  Equation  (47)  over  the  path  length,  L  (shock  tube 
width) ,  yields 


(50) 


(51) 


Combining  Equation  (51)  and  Equation  (46) ,  once  again  assuming 
linear  detector  response,  is  given  as 


E 


(52) 
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The  integral  J  E^z)  dz  is  obtained  experimentally  by  moving  the 
o 

calibration  source  along  the  optical  path  from  0  to  L  and  inte¬ 
grating  a  plot  of  resulting  outputs. 


Equating  Equation  (52)  and  Equation  (49),  the  constant,  K,  is 
given  by 


(z) 


E  (2)  dz 
c 


(53) 


Although  E^  is  a  function  of  A  and  z,  K  is  only  a  function  of  z 
as  can  be  seen  from  Equation  (50).  To  determine  K  for  any  con¬ 
venient  location  such  as  z^,  the  r  E^(z)  dz  need  be  measured  only 

Jo  C 

once  at  any  desired  wavelength,  A^.  Thus,  at  z^, 


(54) 


After  K  has  been  determined  as  above,  G  can  be  computed  from 

A 

Equation  (49)  once  E^  (z^)  is  known  as  a  function  of  A.  The 
measurements  of  (z^)  vs  A  are  obtained  by  placing  the  calibra¬ 
tion  source  at  z^  and  measuring  the  wavelength  dependence  of  the 
output  signal.  The  only  restriction  on  z^  is  that  it  be  within 

the  optical  path.  For  convenience,  the  centerline  of  the  shock 
tube  where  the  entrance  slit  of  the  spectrometer  is  focused  i„ 
always  used.  .'Jote,  at  this  point,  that  the  interferograms  were 
taken  at  the  2/3  span  of  the  test  section  rather  than  the  center- 
line  span  solely  to  minimize  refractive  errors. 

Wavelength  calibration  of  the  spectrometer  was  accomplished 
by  using  a  series  of  calibrated  narrow-band  infrared  filters.  These 
filters  were  inserted  between  the  calibrated  blackbody  and  the 
entrance  optics  of  the  spectrometer.  A  total  of  10  separate  fil¬ 
ters  were  used,  and  a  resulting  average  dispersion  of  0.223  pm/cm 
was  determined. 
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The  absolute  intensity  calibration  was  accomplished  using  a 
bench  setup  with  a  calibrated  blackbody  cavity  located  at  the 
simulated  shock  tube  centerline.  The  centerline  position  can  be 
located  quite  carefully  by  backlighting  the  entrance  slit  of  the 
spectrometer  with  a  strong  diffuse  light  source  and  then  placing 
the  exit  of  the  blackbody  cavity  at  the  image  of  the  entrance 
slit.  The  entrance  optics  are  adjusted  to  give  a  one-to-one 
magnification  of  the  entrance  slit  as  viewed  from  the  shock  tube 
centerline. 

As  a  precaution  against  possible  damage  in  shipment,  the 
linearity  of  the  slit  opening  was  checked  before  the  calibration 
tests.  This  check  was  accomplished  by  measuring  the  spectrometer 
output,  at  a  fixed  wavelength,  for  various  physical  slit  widths. 
Eecause  the  output  intensity  varies  as  the  square  of  the  slit 
width,  a  plot  of  slit  width  vs  square  root  of  intensity  should  de¬ 
fine  a  straight  line  passing  through  the  origin.  Any  deviation 
i'om  this  straight  line  represents  an  irregularity  in  the  slit 
opening.  Failure  to  pass  through  the  origin  represents  an  error 
in  the  zero  setting  of  the  slit  width  dial.  This  check  was  per¬ 
formed,  and  the  straight  line  defined  by  the  data  points  passed 
through  the  origin  of  the  plot  well  within  experimental  error 
limitations. 

The  blackbody  source  used  for  this  calibration  was  operated 
at  a  temperature  of  1000  r  2°C.  Before  each  calibration  run,  the 
cavity  was  allowed  90  minutes  to  stablize,  and  all  tests  were 
accomplished  without  interruption.  Data  points  for  the  curves 

of  E  ( ■ )  vs  >  were  taken  with  z.  corresponding  to  the  centerline 

c  1 

location  and  a  0.500  mm  slit  width.  Results  from  this  calibration 
test  for  the  injection  mixing  (employing  a  2-in. -diameter  sapphire 
viewing  window)  are  shown  in  Figure  25  where  one  can  notice  the 
atmospheric  HO  and  CO^  absorption  bands  between  2.65  and  2.90  juu. 
Figure  26  shows  the  resulting  calibration  curve  of  E^jz^jvs  z^  at 

a  fixed  wavelength  (>  =  3.06  um)  when  the  calibration  souce  is 
moved  along  the  z^  coordinate.  Integrating  the  area  under  this 

curve  provides  the  necessary  input  to  evaluate  the  viewfactor 
constant,  K,  which  for  this  case  turns  out  to  be  0.111  cm-!. 
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Figure  25.  Ec(X)  vs  X  for  the  Injection-Mixing  Repid  Scanning  Spectrometer  Calibration 


c.  Pitot  Probe  -  A  sketch  and  a  photograph  of  the  actual  pitot 
probe  rake  is  shown  in  Figure  27.  Two  pitot  tubes  were  employed 
for  each  run,  and  the  rake  position  could  be  located  in  any  of 
seven  different  basic  positions  shown  in  Figure  27.  At  each 
different  basic  position,  the  probes  could  be  located  in  any  of 
13  different  vertical  positions.  Kulite  strain-gage  type  pres¬ 
sure  transducers  were  used  in  these  probes,  and  each  transducer 
was  calibrated  in  place  before  starting  any  given  survey.  Typ¬ 
ical  calibration  constants  for  these  transducers  are  0.15  mv/torr. 
The  natural  damped  ringing  frequency  of  the  tube  geometry  and 
mounting  cavity  was  computed  on  the  basis  of  the  Helmholtz  res¬ 
onator  formulation  and  a  useful  recording  frequency  of  1.22  kc 
was  indicated.  The  pitot  probe  data  were  all  recorded  using  a 

1  kc  RC  filter  network,  and  the  data  were  reduced  using  the 
Rayleigh  pitot-tube  formulation  given  by  Shapiro  (Ref.  16)  . 

d.  Static  Pressure  Probes  -  Figure  28  is  a  photograph  of  the 
side-wall  mounting  plate  with  the  static  pressure  locations  in¬ 
dicated.  Quartz  crystal  PCB  Model  113-821  transducers  were 
mounted  in  each  of  the  positions  indicated  in  the  photograph. 

A  complete  series  of  static  pressures  was  recorded  by  succes¬ 
sively  repeating  the  run  conditions.  As  will  be  discussed 
in  Section  III,  this  side-wall  static  pressure  is  used  in  con¬ 
junction  with  the  pitot  probe  pressure  to  predict  the  Mach  number 
in  the  exit  region  and  the  mixing  layer  from  experimental  results. 

e.  Temperature  Probe  -  Figure  29  is  a  photograph  of  the  1  mil 
tungsten/tungsten-26%  rhenium  thermocouple  wire.  This  probe 
was  mounted  in  the  same  holder  that  was  used  for  the  pitot  probe 
survey  and  it  could  be  moved  over  a  wide  range  of  locations,  both 
in  the  nozzle  exit  region  and  the  mixing  layer.  Results  from 
each  of  these  three  types  of  probes  and  further  discussion  will 
be  presented  in  Section  III,  when  the  actual  data  are  given. 
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Figure  28.  Side-Wall  Static  Pressure  Transducer  Locations 


Figure  29.  Thermocouple  Probe 


SECTION  III 


EXPERIMENTAL  RESULTS 


As  described  in  Section  II,  two  basically  different  experi¬ 
mental  methods  were  used  to  investigate  the  flow-field  parameters 
within  the  mixing  layer  during  the  time  interval  that  the  overall 
H2  +  F2  +  M  reaction  was  taking  place.  The  experimental  data  from 
these  different  techniques  will  now  be  presented  by  proceeding 
first  with  a  summary  and  selected  details  of  the  H2  injection  re¬ 
sults,  followed  by  a  summary  and  selected  details  of  the  parallel 
shear  flow  mixing  experiments. 

1.  SUMMARY  OF  AVERAGE  DENSITY  PROFILES 

a.  Injection  Mixing  -  A  series  of  experiments  was  conducted 
in  the  injection  shock  tunnel  facility,  where  the  basic  supersonic 
fluorine  carrying  stream  was  unchanged  and  the  hydrogen  injection 
splitter  plate  was  varied  to  determine  its  influence  on  the  over¬ 
all  extent  of  the  mixing  layer.  As  described  in  Section  II,  two 
splitter  plates  were  used  as  nozzle  dividers.  One  of  these  con¬ 
tained  a  series  of  sonic  orifices  for  injecting  hydrogen  into  the 
freestream.  The  other  was  constructed  with  a  series  of  two-di¬ 
mensional  slots  from  which  the  hydrogen  was  injected  into  the 
fluorine  stream. 


Table  V  is  a  summary  of  the  slotted  hydrogen  injection  mixing 
experiments  and  Table  VI  is  a  summary  of  the  orifice  hydrogen  in¬ 
jection  mixing  experiments.  Each  of  these  tables  contains  the 
coded  run  number  (year,  mouth,  day,  shot),  initial  test  gas  mix¬ 
ture,  initial  test  gas  loading  pressure  (Pi),  incident  shock  Mach 
number  |Mgj,  injectant  gas  type,  injectant  gas  pressure 


test  section  static  pressure  ( Pe]  >  degree  of  F;>  dissociation 

|  j,  and  freestream  exit-plane  Reynolds  number  (Re/x).  These 

experiments  were  designed  to  produce  nominal  test  section  static 
pressures  of  5,  10,  and  20  torr  while  maintaining  a  dissociation 
level  of  70%  for  the  fluorine  behind  the  reflected  shock.  As  in¬ 
dicated  in  Tables  V  and  VI,  five  different  injection  conditions 
were  run  for  each  set  of  nominal  conditions.  That  is,  the  mani¬ 
fold  pressures  feeding  the  different  hydi.ien  injection  ports  were 
varied  from  50  to  150  psi  to  distinguish  this  type  of  influence 
on  the  mixing  layer,  and  a  150  psi  helium  injection  case  was  run 
to  distinguish  the  influence  of  a  no-reaction  case  compared  to 
the  reaction  case. 
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Table  V 

INTERFEROGRAM  SUMMARY  FOR  SLOTTED  H2  INJECTION 


Run  No. 

Test  Gas 
Mixture 
by 

Volume 

P1 

torr 

M 

s 

Injectant 

B9 

H 

Pe» 

torr 

“Fa* 

% 

Re/x 

cm"  ^ 

72071802 

507oHe/50%F 

12 

5.58 

None 

- 

15.0 

82.0 

72071803 

507*He/507.F 

12 

5.58 

H2 

150 

15.0 

82.0 

72071901 

507oHe/507,F2 

12 

5.54 

»2 

100 

15.0 

80.0 

72071902 

507He/507.F2 

12 

5.54 

H2 

50 

15.0 

80.0 

1.786x10^ 

72071903 

507„He/507,F 

12 

5.58 

None 

- 

15.0 

82.0 

72071904 

50%He/507<,F2 

12 

5.58 

H2 

150 

15.0 

82.0 

72072101 

507„He/507oF2 

8 

5.34 

H2 

150 

9.8 

75.5 

72072104 

507oHe/507,F2 

8 

5.30 

*2 

100 

9.6 

75.0 

A 

72072106 

50731c /507.F2 

8 

5.30 

H2 

50 

9.6 

75.0 

1.233xl04 

72072401 

507,He/507oF2 

8 

5.34 

None 

- 

9.8 

75.5 

72072402 

507,I!c/507.F2 

8 

5.34 

He 

150 

9.8 

75.5 

72072503 

50/3Ic/507oF2 

n 

5.30 

H2 

150 

4.99 

76.0 

72072504 

507oHc/507„F2 

H 

5.34 

H2 

100 

5.05 

77.0 

0.6449x10** 

72072701 

5  0731  e/  507i’2 

5.30 

H2 

50 

4.99 

76.0 

72072601 

507,He/507,F2 

H 

5.30 

He 

150 

4.99 

76.0 

72072603 

507oHe/507.F2 

u 

5.42 

None 

- 

5.20 

80.0 

Table  VI 


INTERFEROGRAM  SUMMARY  FOR  ORIFICE  H2  INJECTION 


Run  No. 

Test  Gas 
Mixture 
by 

Volume 

P1 

torr 

M 

s 

Injectant 

1 

Pe> 

torr 

V  » 
f2 

% 

Re/x 

-1 

cm 

72071701 

507oHe/50%F 

12 

5.54 

H2 

150 

15.0 

80.0 

72071702 

50%He/50%F2 

12 

5.58 

None 

- 

15.0 

81.0 

72071703 

50YoHe/50%F2 

12 

5.54 

H2 

100 

15.0 

80.0 

1.829xl04 

72071704 

50%He/50YoF2 

12 

5.54 

H2 

50 

15.0 

80.0 

72071801 

50%He/507oF2 

12 

5.54 

He 

150 

15.0 

80.0 

72052503 

F2 

3 

7.18 

H2 

100 

72061502 

507oHe/507„F2 

10 

5.08 

»2 

150 

10.8 

65.0 

72061503 

507oHe/507„F2 

10 

5.35 

■* 

100 

12.2 

74.5 

A 

72061504 

507oHe/507„F2 

10 

5.39 

H2 

50 

12.5 

76.0 

1.381xl04 

72061602 

507=He/507oF 

10 

4.95 

None 

- 

10.0 

61.0 

72061603 

507oHe/507.F2 

10 

5.11 

He 

150 

11.0 

66.0 

72060804 

507oHe/507oF2 

5 

5.19 

H2 

150 

5.8 

71.0 

72060901 

507,He/507„F2 

5 

5.08 

H2 

100 

5.5 

68.0 

72061301 

507oHe/507oF2 

5 

5.27 

H2 

50 

6.0 

74.0 

7.749x10 

72061302 

507.He/50%F2 

5 

5.27 

None 

- 

6.0 

74.0 

72061303 

50%He/507oF2 

5 

5.05 

He 

150 

5.4 

66.0 
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Injection  pressures  quoted  here  were  the  values  measured  at 
the  solenoid  valve  inlet  located  approximately  3  in.  from  the 
outside  wall  of  the  test  section  under  no-flow  conditions.  The 
feed  line  supplying  both  the  orifice  injection  ports  and  slotted 
injection  ports  has  a  cross-sectional  area  more  than  twice  that  of 
the  total  injection  area.  This  area  ratio  should  be  sufficient  to 
ensure  a  minimum  pressure  drop  (not  over  5X)  between  the  large  bot¬ 
tle  supply  source  and  the  small  chamber  orifice  feed  manifold. 

Several  test  cases  indicated  in  Tables  V  and  VI  are  for  runs 
with  no  injection  to  allow  comparisons  of  the  mixing  regions  with 
and  without  the  disturbances  that  are  necessarily  generated  by 
the  injection  process. 

Figure  30  sjjows  six  monochromatic  interferograms  that  were  re¬ 
corded  at  4358  A  for  the  slotted  hydrogen  injection  mixing  case. 

The  run  code  number  is  indicated  on  each  figure  to  allow  a  more 
detailed  identification  of  the  shot  conditions  by  referring  to 
Table  V.  Initial  driver  and  driven  tube  load  pressures  for  each 
of  these  separate  runs  were  held  constant  to  within  experimental 
uncertainties.  As  Table  V  indicates,  incident  shock  Mach  numbers 
for  a  given  series  all  agree  to  with  +  0.04;  hence,  resulting 
test  section  conditions  can  be  assumed  to  be  identical  for  all 
practical  purposes.  The  first  four  interferograms  snown  in  Fig¬ 
ure  30  correspond  to  the  various  different  injection  pressures, 
including  the  150-psi  helium  injection. 

Figure  30e  shows  an  interf erogran  of  the  flow  field  with  no 
injection,  and  Figure  3Qf  is  the  corresponding  zero  field  inter- 
ferogram  taken  immediately  before  running  the  test.  As  can  be 
noted  in  Figure  30e,  the  no-injection  flow  interferogram  shows 
almost  no  fringe  distortion  except  just  downstream  in  the  wake 
region  of  the  splitter  plate.  Even  in  this  wake  region,  the 
fringe  distortion  is  extremely  small  (<0.1  fringe).  Notice  that 
the  term  "fringe  distortion"  has  been  used  here  because  there  is 
in  fact  an  actual  overall  uniform  fringe  shift  between  the  no¬ 
flow  field  (Figure  30f,  test  section  evacuated)  and  the  flow  field. 
Under  these  test  conditions,  fringe  shift  at  the  nozzle  exit  lo¬ 
cation  and  at  all  undisturbed  regions  of  the  flow  field  is  approx¬ 
imately  2.5  fringes,  as  will  be  shown  when  the  final  reduced  re¬ 
sults  are  presented.  This  kind  of  qualitative  result  would  indi¬ 
cate  that  any  overall  fringe  distortion  (as  appears  in  the  first 
interferograms  of  Figure  30)  would  necessarily  result  from  the 
injection  process  and/or  the  chemical  reaction. 
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a)  72061502  150-psi  H2  Injection 


b)  72061503  100-psi  Ho  Injection 


d)  72061603  150-psi  He  Injection 


c)  72061504  50-psi  Hj  Injection 


f)  72061602  No  Flow 


72061602 


Figure  31.  Comparison  of  Orifice  Injection  Snapshot  Interferograms , 


Consider  now  the  interferogram  of  Figure  3Qd,  where  150-psi 
helium  is  used  as  an  injectant.  Here,  it  is  well  known  that  the 
helium  will  not  react  with  the  fluorine  in  the  supersonic  free- 
stream;  hence,  this  record  visually  displays  the  net  result  of 
the  flow-field  characteristics  when  a  simple  nonreacting  injection 
process  occurs.  Notice  in  this  record  that,  just  downstream  of 
the  injection  location,  a  very  small  fringe  distortion  (Fringes 
bending  toward  the  splinter  plate  indicate  an  increase  in  density.) 
does  occur,  but  overall  >.i  effect  is  small  and  is  very  limited 
in  the  direction  nornnj  o  che  plane  of  the  injection  points. 

This  type  of  record  also  provides  direct  experimental  evidence 
that  at  least  no  strong  shock-wave-type  disturbances  result  from 
the  basic  gas  dynamic  injection  process. 

Moving  now  to  the  first  three  interf erograms  in  Figure  30, 
where  the  results  for  the  different  H2  injection  feed  pressures 
are  presented,  a  casual  glance  at  these  records  would  lead  one  to 
conclude  that  the  overall  plume  boundary  downstream  of  the  50-psi 
F  injection  case  is  greater  in  extent  than  is  the  plume  for  the 
150-psi  H2  injection  case.  In  fact,  a  careful  examination  of 
these  records  will  show  this  to  be  true.  Here  we  are  identifying 
the  apparent  plume  boundaries  with  the  rather  sharp  fringe  shift 
that  is  quite  visible  in  Figure  30a.  This  same  type  of  boundary 
can  be  identified  in  each  of  these  interf  erograms,  but  the  close 
fringe  spacing  makes  it  a  little  more  difficult  to  visualize. 

In  general,  one  must  look  at  the  final  reduced  fringe-shift  data 
to  see  this  effect.  In  each  of  the  records,  a  fringe  disturbance 
that  bends  a  given  fringe  toward  the  splitter  plate  indicates  an 
overall  increase  in  the  number  density  relative  to  the  undisturbed 
fringe  condition  that  exists  in  the  freestream  region.  Again,  it 
is  quite  apparent  from  these  records  that  a  larger  fringe  shift 
occurs  within  the  mixing  region  for  the  150-psi  H2  injection  case 
than  occurs  for  either  the  100-  or  the  50-psi  injection  case.  It 
is  also  apparent  from  these  two  cases  that  the  extent  of  the  appar¬ 
ent  reaction  region  is  continued  further  downstream  for  the  150-psi 
H2  injection  case  than  is  true  for  either  of  the  other  two  injection 
pressures.  Detailed  fringe-shift  profiles,  mixing  layer  growth, 
and  density  profiles  will  be  presented  for  the  slotted  injection 
test  cases  following  a  brief  discussion  of  the  orifice  H2  injection 
interferograms. 

We  proceed  next  with  the  interferometric  results  for  the  ori¬ 
fice  H2  injection  splitter  plate.  The  total  area  of  the  series  of 
sonic  orifices  is  less  than  that  available  in  the  slotted  injection 
/  case  by  a  factor  of  two;  hence,  the  comparison  must  be  considered 
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with  this  fact  in  mind.  Using  the  same  general  approach  to  the 
discussion  here  as  was  followed  for  the  slotted  injection  case, 

Figure  31  shows  six  interferograms  with  varying  hydrogen  injection 
pressure,  helium  injection,  and  no-injection  conditions.  Here 
again,  more  complete  details  for  the  run  conditions  are  given  in 
Table  VI  with  identification  made  through  the  run  code  numbers. 

As  was  the  case  for  the  slotted  injection  scheme,  the  no-injection 
interferogram  of  Figure  31e  indicates  that  only  a  very  small  fringe 
disturbance  is  visible  in  the  wake  region  of  the  orifice.  Here 
again,  the  uniform  fringe  shift  between  the  flow  and  no-flow  inter¬ 
ferograms  is  approximately  2.5  fringes.  Considering  next  the  150- 
psi  helium  injection  interferogram  shown  in  Figure  31d,  it  is  clear 
that  no  large  disturbances  are  visible  as  a  result  of  the  gas  dy¬ 
namic  simulation  of  the  injection  process  without  chemical  reactions. 
In  all  these  interferograms,  a  uniform  fringe  shift  or  a  fringe  dis¬ 
tortion  that  causes  a  fringe  to  bend  toward  the  splitter  plate  in¬ 
dicates  an  increase  in  density. 

The  first  three  interferograms  of  Figure  31  show  the  results 
of  varying  hydrogen  injection  pressure  from  150  to  50  psi.  Here 
again,  as  was  the  case  for  the  slotted  injection  scheme,  the  most 
noticeable  effect  Lhat  can  be  observed  from  a  casual  inspection 
is  the  increased  downstream  extent  of  the  mixing  region  character¬ 
istic  of  the  higher  injection  pressure.  In  the  three  cases  shown 
here,  growth  of  the  mixing  region  does  not  appear  to  be  strongly 
influenced  by  changing  injection  pressures,  and  the  plume  bound¬ 
aries  appear  to  be  slightly  larger  for  the  higher  hydrogen  injec¬ 
tion  pressure.  Before  presenting  the  actual  reduced  interfero¬ 
grams  for  these  two  different  hydrogen  injection  techniques,  it 
should  be  pointer  out  that  a  Mach  disc  is  almost  always  visible 
in  these  interferograms  just  downstream  of  the  hydrogen  injection 
location.  This  effect  is  especially  visible  in  Figures  30a,  31a, 
and  31c. 

As  the  previous  discussion  indicates,  a  simple  visual  obser¬ 
vation  of  tne  snapshot  interferograms  does  yield  valuable  quali¬ 
tative  information  about  the  two  different  injection  techniques. 
However,  co  extract  the  maximum  amount  of  information,  a  detailed 
examination  of  the  reduced  interferograms  is  obviously  required. 

Discussion  and  comparison  of  selected  results  from  the  two 
different  injection  cases  will  now  be  presented.  Figures  32,  33, 
and  34  show  the  reduced  interferometric  results  for  six  of  the 
150-psi  hydrogen  injection  cases.  Each  figure  contains  fringe- 
shift  profiles  as  a  function  of  test  section  coordinates  (x,y) 
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Flgvra  32.  Intcrf troMCr&c  Rcaulta.  <.um  72072503  and  72040604 
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>t«r  Result  • ,  Runs  72071803  and  72071701 


for  identical  test  section  conditions  using  the  orifice  and  slot¬ 
ted  11^  injector  splitter  plates.  Again  the  reader  is  reminded 
that  the  maximum  value  of  fringe  shift  (6)  corresponds  to  thu  max¬ 
imum  total  density  level.  The  results  presented  in  Figure  J3  were 
taken  from  the  interf erograms  that  were  presented  and  discussed 
in  connection  with  Figures  30a  and  31a.  The  first  point  to  note 
here  is  that  the  fringe-shift  profiles  are  presented  for  only  one 
half  of  the  mixing  region,  as  the  injection  mixing  experiments  are 
symmetrical  about  the  plane  of  the  injection  ports.  The  nominal 
designed  test  section  static  pressures  in  Figures  32,  33,  and  34 
are  5,  10,  and  15  torr,  respectively,  and,  as  indicated  in  each 
of  these  figures,  the  maximum  fringe  shift  occurs  along  the  center- 
line.  This  indicates  that  the  peak  density  always  occurs  along  a 
plane  containing  the  injection  ports.  The  reduction  in  fringe 
shift  along  this  centerline  plane  moving  downstream  from  the  split¬ 
ter  plate  indicates  that  the  peak  density  is  gradually  decaying  as 
the  diffusive  mixing  increases.  Another  important  feature  appar¬ 
ent  in  Figures  32,  33,  and  34  is  the  increase  in  total  fringe  shift 
in  the  undisturbed  region,  indicative  of  the  increase  in  test  sec¬ 
tion  static  pressure.  This  effect  is  expected  because  the  abso¬ 
lute  fringe  shift  is  a  measure  of  the  density  change  and  each  of 
these  fringe  shifts  is  referenced  to  an  initially  evacuated  test 
chamber.  Before  proceeding  with  a  discussion  of  these  results 
in  terms  of  number  density  profiles,  we  will  next  present  compar¬ 
ison  fringe-shift  results  for  the  50-psi  hydrogen  injection  and 
the  150-psi  hydrogen  injection  experiments. 

Figure  35  shows  the  fringe-shift  profiles  from  12  of  the  ex¬ 
perimental  runs  given  in  Tables  V  and  VI.  Each  of  the  a  through 
f  segments  of  Figure  35  contains  comparable  fringe-shift  profiles 
showing  the  difference  between  150-  and  50-psi  hydrogen  injection 
for  a  given  set  of  test  section  conditions.  Figures  34a,  b,  and 
c  are  results  for  the  slotted  injection  splitter  plate  for  nominal 
test  section  static  pressures  of  5,  10,  and  15  torr,  respectively, 
while  the  c,  d,  and  e  segments  of  Figure  35  are  results  for  the 
orifice  injector  splitter  plate  at  5,  10,  and  15  torr,  respectively. 
The  most  obvious  feature  in  this  comparison  is  that  the  downstream 
influence  of  the  50-psi  hydrogen  injection  is  decreased  as  com¬ 
pared  with  the  150-psi  injection.  That  is,  the  fringe  distortion 
decreases  faster  for  the  lower  injection  pressure  and  for  the  5- 
torr  static  pressure  cases  (Figures  35a  and  d);  the  measurable 
effect  of  the  50-psi  injection  does  not  extend  beyond  1.8  cm  down¬ 
stream  of  the  injection  point,  while,  for  the  150-psi  case,  this 
effect  can  still  be  measured  up  to  3.5  cm  downstream  of  the  injection 
point.  The  other  noticeable  effect  here  is  that  the  overall  widths 
of  the  mixing  regions  for  the  orifice  injection  cases  are  generally 
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Figure  35*.  Comparison  of  Frlnga-Shlft  Profiles  for  150-  and  50-pal  Slotted  H:  InjacCioo  at  5  Torr 
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Flgura  35b.  Comparison  of  Fr  lnga-Sfelf t  Frofllaa  for  150-  and  50-pal  Slottad  N2  Injactlon  at  I  Torr 
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Figure  354.  QapulHi  of  Frli»e-»lft  Frofllee  for  150-  end  50-pel  Orifice  B;  Injection  or 
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slightly  greater  than  the  widths  for  the  slotted  injection  cases. 

This  effect  is  to  be  expected  because  the  pluming  of  a  two-di¬ 
mensional  sonic  slot  is  less  than  that  for  a  three-dimensional 
sonic  orifice. 

A  quantitative  measure  of  the  reaction  zone  thickness  can  be 
obtained  from  the  interferometric  results  by  measuring  the  dis¬ 
tance  between  fringe-shift  minima  at  a  fixed  distance  downstream 
from  the  point  of  hydrogen  injection.  Figure  36  shows  the  results 
obtained  from  the  twelve  interferograms  recorded  at  three  dif¬ 
ferent  injection  pressures  for  the  10-  and  15-torr  nominal  test 
section  pressures  with  both  the  slotted  and  orifice  injector  schemes. 
As  indicated  in  the  inserted  sketch,  the  thickness  "d"  is  measured 
1.5  cm  downstream  of  the  point  of  injection.  As  these  results  in¬ 
dicate,  the  mixing  layer  thickness  increases  with  increased  test 
section  static  pressure,  and,  for  the  orifice  injection  case,  the 
mixing  layer  increases  with  increasing  injection  pressure.  The 
mixing  layer  width  for  the  slotted  injection  case  increases  only 
slightly  with  increasing  injection  pressure  for  the  15-torr  test 
section  pressure  and  shows  an  opposite  trend  for  the  10-torr  test 
section  pressure. 

Averagjd  density  profiles  have  been  computed  at  three  dif¬ 
ferent  down.-1  ream  stations  for  four  of  the  hydrogen  injection  mix¬ 
ing  cases.  L'wo  slotted  and  two  orifice  injection  runs  with  50- 
auc  150-psi  injection,  respectively,  were  chosen  as  representative 
for  this  comparison.  As  explained  in  Section  II,  it  was  impossible 
to  coi.pute  species  number  densities  with  the  multiple  wavelength 
technique;  hence,  to  accomplish  this  objective,  it  was  necessary 
to  make  certain  assumptions  in  this  procedure.  Figure  37  shows 
a  typical  fringe-shift  profile  and  an  indication  of  the  assump¬ 
tions  that  were  made  in  this  data  reduction  procedure.  That  is, 
it  was  assumed  that  the  fringe  shift  along  the  centerline  plane 
was  due  to  hydrogen  molecules;  the  minimum  fringe  shift  within 
the  mixing  layer  was  due  solely  to  hydrogen-fluoride  molecules; 
the  fringe  shift  within  the  nozzle  exit  region  (external  to  the 
mixing  layer)  was  due  to  the  combined  contribution  of  fluorine 
atoms,  fluorine  molecules,  and  helium;  and  the  fringe  shift  at 
intermediate  points  was  due  to  an  averaged  value  of  the  consti¬ 
tuents.  The  averaged  number  density  profiles  presented  for  the 
four  runs  in  Figure  38  were  computed  using  the  theoretical  de¬ 
velopment  of  Section  II  and  the  curves  are  smoothed  values  com¬ 
puted  with  the  assumptions  listed  above. 
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Figure  36  Mixing  Layer  Width  vs  Injection  Pressure 


Figure  37  Fringe  Shift  Showing  Assumed  Species  Distribution 

Jsed  for  i.. imputing  Averaged  Number  Densities 
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As  t he  profiles  of  Figure  38  indicate,  the  hydrogen  number 
density  moving  downstream  along  the  plane  of  injection  decreases, 
thus  indicating  an  increase  in  the  diffusive  mixing  of  the  hydro¬ 
gen.  A  comparison  between  the  slotted  hydrogen  injection  and  the 
orifice  hydrogen  injection  also  shows  the  higher  number  densities 
along  the  centerline  for  the  slotted  injection  scheme.  This  is 
an  indication  of  the  fact  that  the  mass  flow  through  the  slotted 
injector  was  greater  (by  a  factor  of  2)  than  that  through  the 
orifice  injector.  The  HF  number  densities  (minimum  values  in  each 
of  these  profiles)  are  almost  identical  for  each  of  the  different 
injection  pressures,  but  differ  again  by  roughly  a  factor  of  two 
between  the  slotted  and  the  orifice  cases.  The  fact  that  the  HF 
number  densities  for  the  different  injection  pressures  are  almost 
identical  is  also  shown  in  the  infrared  emission  results  that  will 
be  presented  later. 

b.  Parallel  Mixing  -  To  examine  the  resulting  chemically  re¬ 
acting  mixing  region,  the  buffer-driven,  split  shock  tube/shock 
tunnel  facility  described  in  Section  II  was  used  to  generate  two 
adjacent  parallel  supersonic  streams.  The  data  to  be  presented 
in  this  discussion  were  all  taken  in  that  facility  under  a  variety 
of  conditions  designed  to  increase  understanding  of  the  influence 
of  the  exothermic  chemical  reaction  on  the  diffusive  mixing  re¬ 
gions  under  conditions  closely  simulating  the  large  continuous- 
flow  HF  chemical  laser  devices.  Table  VII  is  a  complete 
listing  of  the  interferometric  records  obtained  in  this  effort. 
Specific  experimental  results  presented  later  in  this  discussion 
can  be  identified  in  more  detail  by  referring  to  the  run  code 
numbers,  which  are  shown  with  the  results  in  Table  VII. 

Three  interf erograms  for  the  parallel  mixing  cases  are  shown 
in  Figure  39.  These  particular  records  were  chosen  because  they 
represent  extreme  flow-field  variations  between  the  hydrogen  car¬ 
rying  stream  and  the  fluorine  carrying  stream.  Figure  39a  shows 
a  case  in  which  the  experimentally  measured  nozzle  exit  static 
pressures  are  very  nearly  matched.  Note  that  the  design  nozzle 
exit  static  pressures  and  the  measured  nozzle  exit  pressures  are 
in  considerable  disagreement.  This  point  has  been  discussed  in 
Section  II.  In  this  record,  a  fringe  bending  away  from  the  nozzle 
wall  znd  indicates  a  decrease  in  average  density  compared  to  the 
undisturbed  fieestream  density.  The  record  of  Figure  39a  clearly 
indicates  that  the  average  density  just  downstream  of  the  splitter 
plate  is  much  less  than  the  density  in  both  of  the  nozzle  exits. 
Proceeding  further  downstream  in  the  mixing  layer,  this  density 
profile  tends  to  smooth  out,  and  a  very  weak  pressure  adjustment 
wave  is  clearly  visible  on  the  F2  +  He  nozzle  side. 
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b)  Run  72111602 


c)  Run  72111303 


Figure  39.  Comparison  of  Parallel  Mixing  Interferograms 


85 


Figure  39c  is  a  run  with  conditions  close  to  those  of  the  run 
presented  in  Figure  39a,  except  that  the  hydrogen  has  been  replaced 
by  helium.  This  interferogram  is  then  a  representation  of  the  mix¬ 
ing  region  without  the  chemical  reaction  effects.  It  is  quite  clear 
from  a  comparison  of  these  two  cases  that  the  chemical  reaction  does 
indeed  markedly  modify  the  mixing  region  both  in  extent  and  overall 
density  profiles. 

Figure  39b  shows  the  results  for  a  mixing  layer  in  which  pure 
hydrogen  is  used  as  a  supply  for  one  of  the  nozzles.  In  this  record, 
the  pressure  adjustment  waves  are  very  clear,  and  it  would  appear 
that  much  more  combustion  is  now  occurring.  These  general  remarks 
will  be  extended  and  explained  in  more  detail  in  Section  V,  where 
the  detailed  results  of  the  parallel  mixing  interferometric  dati  are 
compared  with  the  theoretical  results  of  the  mixing  code  to  be  de¬ 
scribed  in  Section  IV. 

The  data  from  these  interf erograms  were  reduced  in  the  same 
manner  as  described  for  the  injection  mixing,  except  that,  in  this 
situation,  there  will  never  be  a  plane  of  symmetry,  and  thus  the 
entire  record  must  be  reduced.  This  data  reduction  has  been  accom¬ 
plished  for  all  shots  listed  in  Table  VII  and  the  fringe-shift  pro¬ 
files  for  six  of  the  cases  are  presented  in  Figures  40,  41,  and  42. 

The  reduced  interferometric  results  presented  in  Figure  40  cor¬ 
respond  to  the  three  interferograms  in  Figure  39.  These  results 
are  plotted  with  a  common  downstream  axis  and  individual  fringe- 
shift  profiles  at  various  distances  downstream  from  the  common  noz¬ 
zle  end  wall.  Additional  details  of  the  flow-field  conditions  can 
be  obtained  by  referring  to  Table  VII  and  the  run  code  number. 

The  fringe  distortion  for  run  72091301  occurs  primarily  in  the 
wake  region  downstream  of  the  common  nozzle  end  wall.  Moving  down¬ 
stream  to  approximately  3  cm,  this  distortion  has  almost  completely 
disappeared.  We  tentatively  identify  this  result  as  indicating 
that  no  appreciable  reaction  has  occurred  in  this  case.  The  fringe- 
shift  profiles  for  run  72111602,  where  a  large  AV  between  the  two 
streams  and  pure  hydrogen  are  used,  do  indicate  that  a  strong  re¬ 
action  has  occurred,  as  shown  by  the  large  fringe  distortion  in  the 
center  region. 

The  fringe-shift  profiles  for  six  additional  parallel  mixing 
experiments  are  shown  in  Figures  41  and  42.  The  most  important 
aspect  of  the  results  for  these  six  cases  is  that  there  is  only  lim¬ 
ited  apparent  reaction  occurring.  This  conclusion  is  based  on  a 
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figur*  40.  Co^mUoq  of  FrlAgo-Shlft  Proflloo  for  Thrtt  Firtllcl  Klxlng  Exporlaonto 
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fi|»n  41.  Ca^MltM  of  Pt  la** -541 ft  Profllaa  for  Thro*  Parallel  MliU|  Crprrlatan 


Figure  42.  Conparlsos  of  Fringe-Shift  Profllu  for  Throe  Parallel  Nixing  Experiments 


comparison  of  the  first  five  cases  with  the  results  for  the  case 
(run  72111303)  where  helium  replaces  hydrogen  in  the  initial  test 
gas  mixture.  That  is,  all  the  results  in  Figure  41  and  42  are 
essentially  identical  and,  since  one  of  these  cases  is  clearly 
a  nonreacting  case,  it  would  appear  that  only  limited  reaction 
could  be  occurring  in  either  of  the  other  experiments.  Additional 
comments  and  comparisons  of  the  parallel  mixing  experimental  re¬ 
sults  with  theoretical  prediction  will  be  presented  in  Section  V. 

2.  SUMMARY  OF  RAPID-SCAN  RESULTS 

A  rapid  scanning  infrared  spectrometer  was  used  in  both  the 
injection  and  parallel  mixing  facilities  to  monitor  the  radiation 
emission  from  excited  HF  molecules  created  in  the  mixing  layer. 

The  12-mm-long  spectrometer  entrance  slit  was  focused  at  the  center- 
line  span  of  the  test  section  with  a  magnification  factor  of  unity. 
Slit  widths  of  0.5  and  1  mm  were  used  in  various  experiments. 

To  eliminate  the  absorption  characteristic  of  quartz,  a  sapphire 
window  was  used  for  a  viewing  port  with  the  spectrometer.  The 
quartz  test  section  windows  used  for  the  interferograms  could  not 
be  used  here  because  they  contain  water  that  absorbs  radiation  in 
the  spectral  region  of  interest. 

In  the  injection  mixing  facility,  a  series  of  separate  experi¬ 
ments,  holding  test  section  flow  conditions  essentially  constant, 
were  performed  as  the  spectrometer  was  translated  in  a  downstream 
direction  moving  away  from  the  hydrogen  injection  ports.  Figure 
43  is  an  interferogram  of  the  test  conditions  (run  72072101)  for 
which  an  extensive  series  of  infrared  emission  scans  was  recorded. 

The  12-mm  field  of  view  is  shown  superimposed  over  the  mixing 
region.  The  slotted  injector  splitter  plate  was  used  in  this 
series  and  25  stations,  varying  from  0.318  cm  to  7.62  cm  down¬ 
stream  of  the  splitter  plate  edge,  were  observed. 

The  exact  starting  time  for  any  given  spectral  scan  could  not 
be  controlled  in  this  experiment  to  within  1.25  msec  due  to  the 
completely  random  starting  signals  derived  from  the  incident  shock 
wave.  That  is,  the  spectral  scans  can  only  be  initiated  when  the 
corner-mirror  rotating  wheel  is  properly  positioned,  and  this  occurs 
every  1.25  msec.  To  ensure  that  a  region  of  uniform  flow  was  ob¬ 
served  by  the  rapid  scanning  spectrometer,  a  separate  time-resolved 
observation  of  the  radiation  emission  was  recorded  using  an  indium 
antimonide  detector  in  conjunction  with  a  2.7-um  narrow-bandpass  fil¬ 
ter.  This  detector  was  on  the  opposite  side  of  the  shock  tunnel 
test  section,  facing  the  entrance  slit  of  the  spectrometer,  with 
a  view  angle  covering  the  entire  mixing  region.  An  oscilloscope 
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record  of  this  2.7-pm  radiation  is  shown  in  the  top  trace  of  Fig¬ 
ure  44,  along  with  a  time  mark  (lower  trace)  from  the  rapid  scanning 
spectrometer  that  allows  one  to  identify  the  particular  scan  re¬ 
corded  during  the  initial  phase  of  the  mixing.  Each  of  the  four 
time  marks  represents  a  scan  from  the  spectrometer,  with  the  left 
time  mark  corresponding  to  the  lower  trace  in  the  radiation  re¬ 
cord  of  Figure  45.  The  marker  placed  on  the  second  trace  from  the 
bottom  of  Figure  45  corresponds  to  the  intensity  versus  wavelength 
record  that  is  consistent  with  the  most  uniform  flow  time,  as  in¬ 
dicated  in  Figure  44.  This  particular  infrared  scan  is  considered 
in  the  data  reduction. 


Figure  43.  Run  72072101  Showing  12-mm  Slit  Height 


91 


Time  — ► 

Figure  45.  Rapid  Scanning  Infrared  Spectrometer  Traces 
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The  radiation  observed  in  Figure  45  is  emitted  chiefly  from 
transitions  in  the  1-0  and  2-1  vibrational  bandr .  The  radiation 
from  2.5  to  3.1  urn  consists  of  the  overlapping  I’-branches  of  the 
two  bands.  The  region  below  2.6  pm  consists  of  radiation  from  the 
partial  overlap  of  the  R-branches  of  the  respective  bands.  Figure 
46  is  a  schematic  representation  of  an  intensity-versus-wavelength 
plot  of  hydrogen  fluoride  for  temperatures  at  which  higher  rotational 
states  are  populated.  This  would  correspond  to  plume  conditions 
where  combustion  takes  place.  The  spectrometer  thus  measures  a 
composite  of  the  radiation  from  several  bands.  There  undoubtedly 
is  some  radiation  contributed  by  the  3-2  band;  however,  the  third 
vibrational  level  is  not  likely  to  be  highly  populated  for  present 
test  conditions. 

As  a  standard  or  figure  of  merit,  the  specific  radiation  in¬ 
tensity,  ,  at  2.795  urn  has  been  used  as  a  means  for  comparing  the 

experimental  infrared  emission  results  corresponding  to  50-  and 
150-psi  hydrogen  injection  pressures  for  the  slotted  splitter 
plate.  This  particular  frequency  was  selected  because  it  corre¬ 
sponds  to  one  of  the  major  lasing  transitions  in  the  2-1  band  of 
HF.  The  two  curves,  shown  on  separate  plots  in  Figure  47,  repre¬ 
sent  a  least-squares  fit  of  the  measured  values  of  (evaluated 

at  2.795  urn)  as  a  function  of  distance  downstream  from  the  splitter 
plate.  It  is  evident  from  Figure  47  that  a  variation  of  injector 
pressure  by  a  factor  of  three  has  virtually  no  effect  on  the  ra¬ 
diation  at  corresponding  downstream  positions.  A  number  of  no¬ 
injection  test  runs  were  interspersed  with  these  runs  at  various 
downstream  locations,  and  in  all  cases  there  was  no  indication 
whatsoever  of  any  radiation  between  1.7  and  4.7  pm.  This  type  of 
check  was  necessary  because  in  some  of  our  earlier  results  we  did 
see  very  unusual  radiation  bands  that  were  later  traced  to  oil 
contamination  in  the  hydrogen  injection  lines.  It  should  also 
be  pointed  out  that  a  concentrated  effort  was  made  to  determine 
if  any  radiation  could  be  observed  at  the  higher  wavelengths. 

These  efforts  did  not  indicate  the  presence  of  radiation  beyond 
the  3.1-^m  level  under  any  conditions. 

The  remainder  of  this  section  will  be  devoted  to  a  general 
discussion  of  the  infrared  emission  radiation  measured  with  the 
parallel  mixing  facility.  The  general  experimental  technique 
used  with  the  parallel  mixing  facility  was  identical  to  that  for 
the  injection  mixing  facility  except  that  the  downstream  survey 
was  not  as  extensive  due  to  the  limited  extent  of  the  sapphire 
viewing  window.  Table  VIII  summarizes  the  experimental  condi¬ 
tions  that  were  used  for  the  parallel  mixing  infrared  spectral 
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Figure  47.  Specific  Radiation  Intensity  vs  Downstream  Distance 
for  Two  Different  Injection  Pressures 


Table  VIII 

PARALLEL  MIXING  INFRARED  SURVEY  DATA 
(ALL  MIXTURES  50%  BY  VOLUME) 


Run  No . 

Matching 

Interferogram 

P2  lower 

H2  +  A, 
torr 

Pi  upper 

He  +  F2, 
torr 

M 

s 

H2  +  A 

M 

s 

He  +  F2 

Focus*, 

in. 

72101101 

20 

6 

4.46 

5.34 

7/32 

72101102 

20 

6 

4.43 

5.26 

5/32 

72101201 

72110903 

20 

6 

4.38 

5.22 

3/32 

72101202 

20 

6 

4.40 

5.34 

9/32 

72101701 

20 

6 

4.27 

5.22 

11/32 

72102001 

17.5 

6 

4.51 

5.34 

11/32 

72102301 

17.5 

6 

4.38 

5.44 

9/32 

72102401 

72110902 

17.5 

6 

4.49 

5.34 

7/32 

72102402 

17.5 

6 

4.49 

5.38 

5/32 

72102403 

17.5 

6 

4.46 

5.34 

3/32 

72102601 

20 

10 

4.48 

4.97 

3/32 

72102602 

20 

10 

4.35 

4.94 

5/32 

72102701 

72110901 

20 

10 

4.35 

5.04 

7/32 

72103001 

20 

10 

4.3r 

5.00 

9/32 

72103002 

20 

10 

4.25 

4.94 

11/32 

*Focus  distance  measured  downstream  from  termination  point  of  the  two  adjacent 
nozzles . 


surveys.  This  table  contains  identification  data  for  15  separate 
runs  recorded  for  three  different  basic  sets  of  test  conditions. 
Three  matching  interf erogram  test  conditions  have  been  reproduced 
in  this  effort,  and  the  run  identification  number  for  the  corre¬ 
sponding  interferogram  is  listed  with  each  of  the  series  of  five 
infrared  spectral  emission  runs.  Here  again,  the  same  figure  of 
merit  (i.e.,  G  at  2.795  urn)  is  used  as  a  means  for  comparing  the 

A 

infrared  emission  as  a  function  of  downstream  location.  Measured 
values  of  G.  at  five  different  downstream  locations  are  shown  in 

A 

Figure  48.  Emission  data  for  each  of  the  three  different  basic 
sets  of  test  conditions  are  shown  on  separate  axes,  and  the  match¬ 
ing  interferogram  run  numbers  are  also  identified.  It  would  appear 
from  this  limited  survey  that  conditions  corresponding  to  inter¬ 
ferogram  72110901  produced  the  best  reacting  mixture  about  1  cm 
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downstream  of  the  splitter  plate.  The  other  two  experimental 
conditions  produced  considerably  less  radiation  and  thus  less  re¬ 
acting  mixture  at  all  positions.  It  should  also  be  noted  that 
no  radiation  was  ever  observed  in  the  parallel  mixing  experiments 
for  wavelengths  greater  than  3.1  um. 

To  gain  some  idea  of  the  combustion  efficiency  in  the  mixing 
layer,  the  total  integrated  radiation  intensity  from  2.6  to  3.1  um 
was  evaluated  for  each  of  the  runs  listed  in  Table  VIII.  The 
radiation  power  density  in  watts  per  cubic  centimeter  was  obtained 
by  evaluating  the  area  within  the  experimentally  obtained  curve 
of  G  vs  X  for  X  values  between  2.6  and  3.1  pm.  Table  IX  9umma- 

A 

rizes  total  radiation  output  determined  in  this  manner  for  each 
of  the  parallel  mixing  experiments.  The  number  of  radiative 
transitions  per  second  from  the  upper  to  the  lower  state  can 
also  be  evaluated  by  assuming  that  the  average  energy  per  tran¬ 
sition  is  that  associated  with  the  2.795-um  line.  This  line  repre¬ 
sents  a  transition  with  an  energy  increment  of  7.1066  x  10-20 
Joules.  Dividing  the  power  density  obtained  experimentally  by 
this  energy  value  then  gives  the  number  of  transitions  in  emis¬ 
sion  occurring  in  a  cubic  centimeter  of  radiating  HF  gas  per 
second.  Self-absorption  was  not  taken  into  account  because  the 
gas  is  assumed  to  be  optically  thin.  The  last  column  of  Table  IX 
lists  the  averaged  number  of  radiative  transitions  obtained  from 
the  measured  radiation  emission  levels. 

3.  SUMMARY  OF  PITOT  PROBE  RESULTS 

A  discrete-point  pitot  probe  survey  of  the  nozzle  exit  and 
mixing  layer  region  was  obtained  using  the  rake  assembly  and  pitot 
probe  geometry  described  in  Section  II.  Two  pitot  tubes  were  em¬ 
ployed  with  every  run,  and  the  basic  rake  position  was  moved  into 
each  of  seven  different  locations.  At  each  separate  basic  posi¬ 
tion,  the  probes  were  located  in  13  different  stations.  This 
type  of  survey  was  conducted  for  the  slotted  injection  splitter 
plate  under  conditions  essentially  identical  to  those  used  for 
the  interferogram  of  run  72072101.  A  corresponding  case  was  also 
run  under  conditions  with  no  injection.  The  data  that  were  taken 
with  the  pitot  tubes  in  a  fixed  position  (with  the  base  being  moved 
across  the  width  of  the  test  section)  indicated  that  the  flow  field 
was  indeed  uniform  to  within  the  accuracy  of  these  measurements. 
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Figures  49  and  50  show  the  actual  data  points  and  a  smoothed 
curve  connecting  them  for  the  150-psi  hydrogen  injection  case  and 
the  zero-injection  case  for  this  extensive  series  of  runs  in  which 
the  basic  conditions  remained  unchanged.  As  indicated  in  Figure  49 
for  the  150-psi  hydrogen  injection  case,  the  minimum  pitot  pressure 
is  observed  along  the  centerline  plane  of  injection  and  then  in¬ 
creases  to  a  peak  value  within  the  mixing  layer  before  decaying  to 
a  freestream  value  out  near  the  nozzle  walls.  This  same  type  of 
^asic  variation  in  pitot  pressure  can  be  observed  for  the  no-in¬ 
jection  case,  but  the  variations  are  much  smaller.  The  average 
value  of  the  pitot  pressure  near  the  wall  in  each  different  rake 
position  is  approximately  16  psi.  Table  V  gives  more  complete  de¬ 
tails  of  the  reference  conditions  (i.e.,  run  number  72072101)  for 
which  this  survey  was  recorded.  As  will  be  discussed  later,  an 
extensive  static  pressure  survey  for  these  same  test  section  con¬ 
ditions  was  conducted  and  an  average  test  section  static  pressure 
of  30  torr  has  been  measured.  Using  these  two  values  with  the 
classical  Rayleigh  pitot  tube  formula  and  the  theoretical  value 
of  y  for  the  reflected  shock  region,  one  obtains  an  exit  Mach 
number  of  4.5  for  this  nozzle.  The  actual  design  Mach  number  for 
this  nozzle  was  5,  so  this  type  of  agreement  is  considered  to  be 
good . 


Table  IX 

TOTAL  INTEGRATED  RADIATION  FROM  2.6  TO  3.1  pm 


Run 

No. 

Interferogram 

Power 

w/cm3 

Transitions/cm3  sec 

72102601 

64.5  x  10~4 

9.06  x  1016 

72102602 

55.5  x  10"4 

7.80  x  1016 

72102701 

72110901 

79.5  x  10~4 

11.18  x  1016 

72103001 

85.6  x  10'4 

12.03  x  1016 

72103002 

122.2  x  IQ*4 

17.20  x  1016 

72102001 

69.3  x  IQ"4 

9.75  x  1016 

72102301 

75.0  x  10"4 

10.53  x  1016 

72102401 

72110902 

57.6  x  10-4 

8.10  x  1016 

72102402 

86.8  x  10-4 

12.05  x  1016 

72102403 

63.0  x  10"4 

8.85  x  1016 

72101101 

65.5  x  10-4 

9.20  x  1016 

72101102 

56.9  x  10'4 

8.00  x  1016 

72101201 

72110903 

54.0  x  10-4 

7.60  x  1016 

72101202 

57.4  x  10-4 

8.06  x  1016 

72101701 

73.8  x  10“4 

10.38  x  1016 

99 


y  (inches) 


Nozzle  Wall 


4.  SUMMARY  OF  STATIC  PRESSURES 


A  complete  series  of  side-wail  static  pressure  measurements 
has  been  recorded  for  one  basic  set  of  test  conditions  in  the  in¬ 
jection  mixing  facility  where  the  nozzle  exit  conditions  and  hydro¬ 
gen  injection  pressures  corresponding  to  interferometric  data  run 
72072101  were  reproduced.  This  complete  static  pressure  survey 
was  provided  for  this  particular  set  of  conditions  because  they 
also  correspond  to  the  conditions  for  which  the  extensive  IR  scans 
and  pitot  probe  surveys  were  obtained.  That  is,  we  have  chosen 
a  single  set  of  run  conditions  (viz,  72072101)  and  then  proceeded 
with  the  most  extensive  possible  surveys  with  each  of  the  basic 
experimental  techniques.  Static  pressures  have  also  been  re¬ 
corded  at  only  one  station  in  the  nozzle  exit  region  for  each 
of  the  other  basic  sets  of  run  conditions  in  both  the  injection 
mixing  facility  and  the  parallel  shear  flow  facility. 

Figure  51  is  a  sketch  of  the  injection  mixing  nozzle  blocks 
showing  the  relati.  •  positions  of  the  static  pressure  ports  at 
which  the  full  static  pressure  survey  was  recorded.  Ten  separate 
stations,  as  shown,  were  used  to  obtain  these  data.  The  results 
of  a  no-injection  and  a  150-psi  H2  injection  are  shown  in  Figures 
52  and  53.  Each  of  these  figures  shows  the  average  value  of  at 
least  three  different  runs  for  all  positions  recorded.  At  the 
time  of  this  particular  survey,  we  were  using  pressure  transducers 
that  had  been  in  service  for  a  long  time  and  the  noise  level  for 
these  gages  was  poor.  New  transducers  were  installed  after  these 
data  were  obtained,  and  they  provide  a  much  better  signal-to-noise 
ratio. 
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Figure  51.  Schematic  of  Injection  Mixing  Nozzle  Blocks  Shoving 
Location  of  Static  Pressure  Ports 
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Figure  52  is  a  typical  scope  record  obtained  with  these  new 
gages  where  two  transducer  locations  on  either  side  of  the  in¬ 
jector  plate  are  displayed  in  the  two  separate  traces.  These  re¬ 
cords  indicate  that  the  static  pressures  at  stations  D  and  E  for  the 
injection  nozzle  exits  are  indeed  quite  uniform. 


Time  - ► 

Figure  52.  Static  Side-Wall  Pressure  Measurements 


Figure  53  shows  the  results  obtained  from  a  series  of  test  *uns 
for  the  initial  8-torr  pressure  of  a  50%  by  volume  mixture  of 
He  +  F2  with  no  injection  and  with  an  average  initial  shock  Mach 
number  of  5.4.  As  the  values  shown  here  indicate,  the  average  test 
section  static  pressure  is  approximately  30  torr  at  all  stations 
except  C2,  which  is  located  just  downstream  of  the  splitter  plate. 
As  indicated  and  discussed  in  Section  II,  the  nominal  test  section 
pressure  computed  for  these  run  conditions  is  actually  10  torr; 
hence,  the  measured  value  is  some  three  times  higher  than  we  had 
anticipated.  The  specific  reasons  for  this  disagreement  are  dis¬ 
cussed  in  Section  II. 
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Figure  54  shows  the  measured  static  pressure  values  for  these 
same  run  conditions  with  a  150-psi  hydrogen  injection  case.  As 
the  values  indicate,  the  static  pressure  downstream  of  the  noz¬ 
zle  exit  is  increased  by  factor  of  1.8  over  the  zero-injection 
case.  Values  presented  for  this  case  are  averaged  results  of  two 
different  repeated  runs,  and,  here  again,  the  noise  level  was  not 
as  good  as  we  can  now  expect  to  obtain  with  the  new  gages.  Values 
for  the  static  pressure  measurements  obtained  with  parallel  mixing 
configurations  are  given  in  Table  VII,  along  with  the  summary 
conditions  for  the  interferometric  data. 

5.  SUMMARY  OF  TEMPERATURE  MEASUREMENTS 

Two  custom-made  tungsten/ tungsten-rhenium  thermocouple  probes, 
described  in  Section  II,  were  obtained  with  the  intent  of  conduct¬ 
ing  a  complete  total  temperature  survey  of  the  mixing  region.  Due 
to  delayed  delivery  and  the  extremely  long  rise  time  (^4  msec)  of 
the  2-mil-diameter  probe,  only  a  limited  number  of  data  points 
were  obtained.  The  test  conditions  for  this  survey  were  identical 
to  those  for  the  interferometric  data  run  72072101,  and  a  hydrogen 
injection  pressure  of  150-psi  was  used.  The  results  of  this  sur¬ 
vey  presented  in  Figure  55  are  certainly  questionable  in  terms 
oi  the  final  absolute  value  because  the  probe  was  not  shielded 
and  no  attempt  was  made  to  account  for  radiation.  These  results 
do  indicate  a  very  symmetrical  temperature  profile  over  the  mix¬ 
ing  region,  with  the  peak  measured  temperature  occurring  at  a 
station  1  in.  above  the  centerline  plane. 

Four  different  runs  with  the  1-mil-diameter  probe  have  been 
obtained,  and,  as  anticipated,  the  rise  time  for  this  configuration 
is  approximately  2  msec,  which  is  well  within  usable  limits.  This 
particular  probe  survived  four  runs  before  being  destroyed,  and 
because  we  nave  just  identified  a  new  source  for  replacement  junc¬ 
tions  a",  less  than  $10  per  junction,  we  feel  that  a  very  accurate 
temperature  survey  can  now  be  obtained  with  small  additional  effort 
to  provide  an  effective  shield. 
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Probe  2  in.  Downstream 

2-mil  tungsten/ tungsten-267,,  rhenium 


Figure  55.  Temperature  Survey  with  150-psi  Slotted  Injection 
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SECTION  IV 


THEORY  AND  ANALYSIS  OF  PARALLEL  TURBULENT  MIXING  LAYERS 


1.  TURBULENT  TRANSPORT  MODELS 

Historically,  turbulence  in  compressible  flow  has  been  analyzed 
by  employing  extensions  of  concepts  and  methods  obtained  from  ex¬ 
perience  with  turbulence  in  incompressible  flow.  Advances  in  the 
understanding  of  turbulence  in  supersonic  flows  have  stemmed  almost 
entirely  from  the  pioneering  work  of  Kovasznay  (Ref.  17)  and  Mor- 
kovin  (Ref.  18)  in  the  development  of  the  hot  wire  anemometer  tech¬ 
nique  as  a  reliable  method  for  quantitative  observation  of  compres¬ 
sible  turbulence. 

Interpretation  of  experimental  results,  such  as  those  obtained 
by  Morkovin  (Ref.  19),  Kovasznay  (Ref.  20),  Kistler  (Ref.  21),  and 
Demetriades  (Ref.  22),  has  established  that  in  nonreacting  compres¬ 
sible  turbulent  flow,  the  basic  mechanisms  for  turbulent  transport 
are  essentially  the  same  as  those  for  incompressible  turbulent  flow, 
as  long  as  the  fluctuation  velocities  are  subsonic.  Such  results 
lead  to  the  conclusion  that  density  fluctuations  and  pressure  fluctu¬ 
ations  in  nonreacting  compressible  turbulent  mixing  layers  do  not 
contribute  significantly  to  the  turbulent  transport,  and  compressi¬ 
bility  is  felt  only  through  the  temperature  and  compression  effects 
on  the  mean  density. 

On  the  other  hand,  i:  has  been  demonstrated  by  Glass  (Ref.  23) 
that  intense  sound  impingement  can  drastically  change  the  character¬ 
istics  of  turbulent  jets.  Sanderson  and  Steel  (Ref.  24)  have  shown 
that  exothermic  reactions  may  induce  a  dramatic  increase  in  turbu¬ 
lent  mixing  layer  growth  rate.  Experimental  observations  of  com¬ 
bustion  in  turbulent  boundary  layers  by  Wooldridge  and  Muzzy  (Ref. 

25)  show  a  similar  effect,  in  contradiction  to  earlier  results  by 
Kulgein  (Ref.  26). 

Engineering  analyses  of  reacting  compressible  turbulent  mixing 
have  traditionally  employed  phenomenological  models  of  turbulent 
transport  properties  adapted  from  nonreacting  incompressible  ex¬ 
perience,  without  attempt  to  account  for  a  direct  influence  of 
exothermic  energy  on  the  transport  processes.  An  extensive  review 
and  comparison  of  such  transport  models  as  employed  in  nonreacting 
turbulent  flows  has  been  published  by  Harsha  (Ref.  27). 
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Such  phenomenological  models  of  turbulent  transport  may  be 
termed  "eddy  viscosity"  models,  in  that  they  attempt  to  relate 
the  turbulent  shear  stress  to  local  flow  conditions  in  a  manner 
that  is  analogous  to  the  dependence  of  laminar  viscosity  on  local 
flow  conditions.  It  is  well  known  that  the  mechanics  of  the  gen¬ 
eration  of  turbulence  or  change  in  the  character  of  turbulence  can¬ 
not  be  defined  uniquely  by  mean  local  flow  conditions  alone.  None¬ 
theless,  such  approaches  have  exhibited  much  success  in  predicting 
the  gross  behavior  of  mixing  flows,  particularly  when  sufficient 
experimental  data  exist  to  permit  empirical  correlations  of  the 
transport  models. 

Over  the  past  20  years  or  so  many  workers  have  attempted  to 
construct  models  of  turbulent  flow  based  on  closure  of  the  system 
of  equations  of  turbulent  motion  at  a  higher  order  than  that 
treated  by  the  phenomenological  eddy  viscosity  approach.  Recently, 
such  approaches  have  yielded  excellent  results  in  the  prediction 
of  compressible  turbulent  boundary  layer  phenomena.  The  work  of 
Donaldson  (Ref.  28),  for  example,  has  proved  particularly  fruitful. 
He  closes  the  system  of  equations  by  mathematically  modeling  the 
appropriate  terms  in  the  equations  for  the  second  order  correla¬ 
tions  of  fluctuating  quantities  so  that  the  terms  possess  the 
correct  tensor  qualities  but  are  tractable  to  computation.  How¬ 
ever,  such  methods  are  not  sufficiently  well  in  hand  to  permit 
their  application  to  the  general  problem  of  mixing  with  chemical 
reactions. 

Consequently,  in  the  development  of  a  practical  analytical 
method  for  the  routine  engineering  prediction  of  turbulent  mixing 
flow  with  chemical  reactions,  one  has  little  choice  but  to  employ 
one  of  the  phenomenological  or  eddy  viscosity  models.  Our  approach 
is  to  select  a  reasonable  eddy  viscosity  model  and  to  experimentally 
determine  the  adjustable  parameters  of  that  model  over  a  range  of 
conditions  that  are  pertinent  to  the  hydrogen-fluoride  chemical 
laser.  This  approach  is  a  continuation  of  work  that  began  some 
years  ago  with  the  supersonic  mixing  of  hydrogen  and  oxygen  (see 
Ref.  24,  29  thru  37).  We  have  selected  Prandtl’s  third  eddy 
viscosity  model  [see  for  instance  Hinze  (Ref.  38)].  For  compres¬ 
sible  flow  we  find  that,  with  the  boundary  layer  approximation, 
the  gradient  of  the  turbulent  shear  stress  may  be  written 
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(55) 


where,  from  Prandtl's  model 

<  *  k  b  (  u  -  u  .  )  (56) 

\  max  rain ) 

Here,  b  is  the  mixing  layer  width  and  k  is  a  constant.  This  formu¬ 
lation  yields  excellent  theoretical  results  for  two-dimensional, 
two-stream  mixing  layers  and  co-axial  stream  mixing  layers  [see 
Harsha  (Ref.  27)]. 


The  available  experimental  results  for  two-stream  mixing  layers 
with  chemical  reactions  indicate  that,  in  terms  of  gross  features, 
such  layers  exhibit  similarity.  In  that  case,  the  two-dimensional 
equations  of  motion  may  be  written  in  terms  of  a  single  independent 
variable,  -  =  -•  y/x,  where  o  is  a  "universal"  length  scale  or  mixing 
layer  growth  rate.  At  the  mixing  layer  boundary 


=  —  =  constant 
b  x 


(57) 


and 


b  = 


=  cx 


The  eddy  viscosity  model  then  may  be  written  as 
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*Note  that  in  laminar  flow,  t 

’  xy 
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It  may  be  shown  that  this  formulation  is  a  "natural"  definition 
for  the  eddy  viscosity  with  either  incompressible  or  compressible 
mixing  layers,  as  long  as  the  similarity  assumption  is  valid. 

The  eddy  viscosity  is  employed  in  a  numerical  mixing  analysis 
program  to  predict  contours  of  the  dependent  variables  across  the 
mixing  layers.  This  procedure  is  discussed  in  detail  in  Subsection 
2  of  this  section.  The  parameter,  c,  is  adjusted  until  the  pre¬ 
dicted  and  observed  contours  are  in  reasonable  agreement.  The  de¬ 
pendent  variable  most  easily  observed  for  two-dimensional,  super¬ 
sonic,  reacting  mixing  layers  has  been  the  index  of  refraction 
for  the  gas  mixture. 

The  study  attempted  to  find  reasonable  correlations  for  the 
value  of  z  thus  determined  in  terms  of  some  easily  computed  param¬ 
eter.  In  the  case  of  co-axial  compressible  mixing,  the  momentum 
deficit  across  the  layer,  or  pu  ratio,  is  an  important  parameter. 
Clearly  other  factors,  such  as  exit  Reynolds  number,  also  must  be 
involved  [see  for  instance  Harsha  (Ref.  27)].  Such  correlation 
with  the  parallel  stream  data  has  not  been  found. 

The  amount  of  compressible  stream  data  available  is  not  large; 
however,  there  is  a  considerable  amount  of  low-speed  parallel  stream 
data  in  the  literature  with  a  variation  of  the  ratio  of  densities 
between  the  streams  as  large  as  15.  Such  data  for  o  correlate 
reasonably  well  with  a  velocity  parameter  defined  as 


u  -  u 
max  min 


u 

max 


+  u 


min 
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particularly  for  values  of  8  greater  than  0.3.  The  wide  scatter  of 
o  values  for  values  of  8  less  than  0.3  indicates  that  as  the  stream 
velocity  difference  decreases,  the  boundary  layers  on  the  stream 
separating  surface  become  dominant  in  controlling  the  mixing. 

The  values  of  a  obtained  in  previous  work,  both  with  and 
without  combustion,  are  shown  in  Figure  56  as  a  function  of  8, 
together  with  low-speed  data  by  other  laboratories.  Recent  re¬ 
sults  by  other  laboratories  are  plotted  for  compressible  reacting 
flows.  The  values  of  o  derived  from  this  latter  data  are  not  as 
reliable  as  those  obtained  in  Martin  Marietta's  previous  work, 
because  only  limited  information  was  available  on  the  specific 
experimental  conditions  and  results  from  these  laboratories. 
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Mixing  Parameter 
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The  value  of  that  best  fitted  a  similar  profile  to  the  available 
data  was  determined.  Although  the  points  show  severe  scatter, 
there  is  strong  evidence  that  at  low  values  of  the  velocity  param¬ 
eter  (,;  0.3),  the  mixing  layer  growth  rate  is  strongly  affected 

by  exothermic  reactions.  In  such  cases,  a  a  value  near  20  is  used 
rather  than  the  noncombustion  value  that  is  nearer  40.  The  large 
scatter  of  noncombustion  results  for  low  values  of  6  is  felt  to  be 
a  result  of  initial  boundary  layer  dominance. 

Note  the  Drewry  (Ref.  39)  and  Shackleford  (Ref.  40)  data  are 
related  to  the  observations  of  a  local  species  mass  fraction  or 
number  density,  for  which  the  transport  coefficient  is  the  eddy 
diffusivity  or  conductivity.  The  ratio  of  eddy  viscosity  to  eddy 
conductivity  is  the  turbulent  Prandtl  number,  Pr^.  Observed  values 

of  the  turbulent  Prandtl  number  range  from  about  0.6  to  about  1. 

The  Drewry  and  Shackleford  data  shown  on  Figure  56  have  not  been 
corrected  for  turbulent  Prandtl  number,  which  is  equivalent  to 
setting  it  to  one.  If  the  value  of  Prt  is  taken  as  different  from 

unity,  the  values  of  c  will  increase  by  the  factor  l/yJPr^ . 

Martin  Marietta's  previous  work  (Ref.  24)  did  not  numerically 
produce  realistic  index  of  refraction  contours  without  accounting 
for  the  boundary  layers  produced  on  the  stream  separating  surface 
(see  Fig.  57).  Note  that  most  of  Martin  Marietta's  previous  data 
were  obtained  for  values  of  the  velocity  parameter,  3,  less  than 
0.3,  which  implies  a  boundary  layer  influenced  mixing  layer.  The 
past  work  simply  employed  initial  dependent  variable  contours  that 
were  consistent  with  observed  boundary  layer  thicknesses.  However, 
this  procedure  does  not  account  for  the  transition  from  boundary 
layer  controlled  mixing  to  fully  formed  mixing  layer  flow. 

In  future  work,  Martin  Marietta  plans  to  account  for  this 
transition  in  the  following  way.  Transition  to  turbulent  mixing 
is  assumed  if  the  momentum  thickness  Reynolds  number  exceeds  150 
at  the  initiation  of  a  free  shear  layer  [see  Dryden  (Ref  41)  or 
Landau  and  Lifshitz  (Ref.  42)].  The  velocity  defect  employed  in 
the  transport  properties,  however,  is  taken  initially  as  that 
across  the  boundary  layer  profile.  As  the  mixing  progresses,  the 
velocity  minimum  resulting  from  the  nozzle  wall  gradually  increases. 
Once  this  defect  becomes  equal  to  the  velocity  difference  across 
the  mixing  layer,  the  usual  velocity  difference  is  employed. 
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The  correlation  given  in  Figure  56  appears  to  serve  adequately 
for  nonreacting  flows,  but  the  data  on  reacting  flows  are  insuf¬ 
ficient  to  allow  a  curve  to  be  drawn  for  that  category  of  mixing 
layers.  As  was  stated  earlier,  the  data  indicate  a  strong  effect 
of  combustion  energy  release  under  conditions  where  this  energy 
is  significant  in  relation  to  the  flow  kinetic  energy.  In  view  of 
this  fact,  it  seems  reasonable  to  expect  that  any  complete  corre¬ 
lation  must  include  a  parameter  to  express  the  relative  magnitudes 
of  these  two  energy  sources,  as  well  as  a  parameter  to  measure  the 
upstream  boundary  layer  effects.  Martin  Marietta  has  been  able 
to  correlate  the  MMC  and  TRW  combustion  data  (Ref.  36  and  40)  for 
6  less  than  0.3  with  the  following  equation: 

0  =  _ 1-35  (Re/5)a _ 

1  +  0.002  [(1  t  C  •  Re3)  (1  -  43)]1*35 

where  C  *  ratio  of  available  chemical  energy  to  the  fluid  kinetic 
energy , 

Re  =  maximum  nozzle  exit  Reynolds  number, 

a  *=  1/2  for  laminar  boundary  layers  and  1/5  for  turbulent 
boundary  layers, 

i  =  velocity  parameter. 

The  correlation  is  incomplete,  however,  because  it  does  not  readily 
extend  to  nonreacting  mixing  layers  (C  ■  0)  as  it  should.  More 
reacting  mixing  layer  data  will  be  required  to  develop  a  complete 
correlation. 

2.  ANALYTICAL  TECHNIQUES 

The  finite  difference,  finite  chemical  rate  analysis  employed 
in  this  work  is  an  adaptation  of  one  developed  by  Edelman  and 
Fortune  (Ref.  43).  This  analysis  and  its  associated  computer  code 
were  originally  written  for  the  mixing  and  burning  of  an  RP-l/lox 
system.  Modifications  that  were  made  involve  addition  of  an 
equilibrium  chemistry  option  for  hydrogen/fluorine  combustion, 
alteration  of  the  existing  nonequilibrium  chemistry  subroutines 
to  reflect  the  hydrogen/fluorine  composition  of  interest,  and 
incorporation  of  certain  features  to  facilitate  consideration 
of  boundary  layers  that  form  on  surfaces  upstream  of  the  mixing 
region.  Further  discussion  of  the  mixing  analysis  will  be  divided 
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into  three  parts:  (1)  formulation  of  the  governing  equations  in 
finite  difference  form,  (2)  the  nonequilibrium  chemistry  computa¬ 
tion  scheme,  and  (3)  the  computer  program  that  employs  (1)  and 
(2)  to  achieve  the  desired  solution. 

a.  Finite  Difference  Analysis  of  the  Flow  Field  -  Certain 
basic  assumptions,  important  to  the  analysis  described  here,  are: 

1)  Flow  geometry  is  restricted  to  two-stream,  parallel  mix¬ 
ing,  either  axisymmetric  or  two-dimensional; 

2)  The  flow  field  is  fully  turbulent;  molecular  transport 
mechanisms  are  assumed  negligible  in  comparison  to  their 
turbulent  counterparts; 

3)  Turbulent  Prandtl  and  Lewis  numbers  are  constant; 

4)  Turbulent  transport  may  be  expressed  in  terms  of  an  eddy 
viscosity/mixing  length  concept; 

5)  The  perfect  gas  equation  of  state  is  adopted  for  both 
individual  species  and  mixtures; 

6)  The  boundary  layer  form  of  the  Navier-Stokes  equations  is 
assumed  applicable,  pressure  gradient  across  the  mixing 
region  is  assumed  to  be  negligible; 

7)  Transfer  of  mass  by  thermal  or  pressure  diffusion  is 
neglected . 

By  applying  these  simplifying  restrictions,  the  governing 
differential  conservation  equations  for  the  mixing  region  become: 

Mass : 
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Energy: 


where  N  =  0  for  two-dimenaional  geometry  and  N  **  1  for  axisymmetric 
geometry.  The  appropriate  initial  and  boundary  conditions  are: 


x  -  0 

y  -  0 


u  -  u(y) ,  -  ai(y),  T  -  T(y) 

—  »  —  m  -  0 

9y  9y  9y 


y  -  °°  u  *  Vig,  -  a±  ,  T  *  Tg  (67) 

e 

With  the  inclusion  of  the  equation  of  state,  these  equations  repre¬ 
sent  a  complete  set,  with  the  exception  of  the  chemical  production 
terms  that  are  calculated  separately. 


Solution  of  these  equations  in  the  physical  (x,y)  plane  may  be 
undertaken  directly.  However, for  purposes  of  this  finite  differ¬ 
ence  solution,  it  is  convenient  to  employ  a  stream  function  trans¬ 
formation  to  the  (x,i|>)  plane,  ip  being  the  stream  function  defined 
by 
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(68) 
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This  transformation  satisfies  the  mass  conservation  equation  iden¬ 
tically  and  is  straightforward  in  application.  Note  that  under 
this  transformation 


(69) 


Substitution  of  the  stream  function  into  the  momentum,  energy,  and 
species  conservation  equations  yields: 


Momentum: 
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bx  CU  bx  N 
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Energy: 
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Species : 
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A  finite  difference  scheme  is  employed  to  formulate  these 
parabolic  differential  equations  in  a  form  suitable  for  a  numer¬ 
ical  initial/boundary  value  marching  solution.  The  finite  dif¬ 
ference  mesh  is  composed  of  streamlines  spaced  at  a  specified 
Interval  and  normal  surfaces  spaced  at  varying  intervals,  Ax,  as 
determined  by  a  stability  criterion.  This  stability  criterion 
arises  from  the  choice  of  an  explicit  difference  scheme  which 
allows  solution  in  each  streamwlse  step  without  need  for  itera¬ 
tion.  Forward  differences  are  employed  for  x  derivatives  and 
central  differences  for  ip  derivatives.  For  a  dummy  dependent 
variable  F,  these  derivative  forms  are 
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where 
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and  the  mesh  is  as  follows 
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The  explicit  form  of  differencing  has  the  disadvantage  of  relat¬ 
ing  Ax  to  A1;.  This  may  lead  to  long  solution  times  when  the  total 
streamwise  distance  to  be  traversed  is  much  larger  than  a  charac¬ 
teristic  dimension  in  the  cross  stream  (i|/  or  y)  direction,  e.g., 
a  boundary  layer  thickness. 

The  final  finite  difference  equations  for  u  ,,  ,  T  . ,  ,  and 

n  n+l,m  n+l,m’ 

u.  are  lengthy  and  will  not  be  derived  here.  However,  they 

1n+l , m 

are  easily  obtained  by  direct  application  of  the  above  derivative 
forms  to  the  governing  equations.  It  is  important  to  note,  how¬ 
ever,  that  these  equations  are  uncoupled  and  present  no  difficul¬ 
ties  for  machine  solution  once  the  chemical  species  production 
rates  have  been  obtained.  This  phase  of  the  analysis  will  be 
discussed  next. 

b.  Nonequilibrium  Chemistry  Analysis  -  Finite  rate  analysis 
is  the  most  realistic,  but  also  the  most  difficult  of  the  various 
chemical  analysis  options.  One  must  provide  a  model  for  the 
chemical  processes  that  can  be  used  repeatedly  in  the  finite  dif¬ 
ference  procedure  without  causing  computer  run  times  to  become 
prohibitive.  To  achieve  this,  the  numerous  reactions  that  form  the 
actual  (and  often  undefined)  chemical  chains  must  be  replaced  by 
a  limited  set  of  important  controlling  reactions  and  appropriate 
reaction  rate  constants.  Such  approximations  meet  with  success 
or  failure  depending  on  how  well  the  controlling  reactions  are 
selected  and  how  well  their  rate  constants  are  known. 

For  analysis  of  the  present  hydrogen/fluorine  system,  the 
species  selected  for  consideration  are  H,  H2,  F,  F2,  and  HF.  Pro¬ 
vision  is  also  made  for  the  presence  of  nonreacting  species  He  and 
A  acting  as  diluents.  The  reactions  selected  as  representative 
of  the  important  chemical  processes  and  their  respective  rate  con¬ 
stants  are  shown  in  Table  X.  These  reactions  were  selected  from 
those  given  by  Cohen  (Ref.  44),  with  the  production  of  excited 
HF  specifically  neglected. 

The  solution  procedure  adopted  for  our  analysis  follows  that 
used  by  Edelman  and  Fortune  (Ref.  43),  which  derives  in  turn  from 
the  work  of  Moretti  (Ref.  45),  and  DeGroat  and  Abbett  (Ref.  46). 

An  outline  of  the  major  steps  in  this  procedure  will  be  given  here, 
with  further  details  available  from  References  45  and  46. 
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Table  X 


CHEMICAL  REACTION  AND  RATE-CONSTANT  EQUATIONS 


Reaction 

Rate  Constant 

F2  +  M  ~  •  2F  +  M 

kf  »  5  x  1013  •  10-7740/T 

H2  +  M  _ _ 1  2H  +  M 

i 

•-* 

o 

»— • 

CD 

• 

H 

1 

t— • 

HF  +  M  ,  '  H  +  F  +  M 

kf  -  1.2  x  1018  •  T-1  •  10_32180/T 

F  +  H2  “  *  HF  +  H 

kf  -  2.69  x  1013  •  10_350/T 

H  +  F?  ^  *  HF  +  F 

kf  -  4.8  x  1012  •  10"525/T 

The  following  thermodynamic  relations  for  gaseous  mixtures 
are  assumed: 

State: 


P  -  p  R 


(76) 


Enthalpy; 

k 

h  «  £  «1hi  (77) 

i-1 

Specific  Heat: 


k 


As  a  matter  of  convenience  for  the  subsequent  development,  a  defi¬ 
nition  is  presented  of  a  molar  concentration  of  species  i 


(78) 
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(79) 


and  a  third-body  concentration 
k 

y  y,  yi  (80) 

i-1 


where  the  summation  includes  all  species  present  in  the  gaseous 
mixture. 


The  various  chemical  reactions  of  Table  X  can  be  written  in 
the  general  form 


k 

E 

i-1 


v!  M, 


k 

E 

i-1 


M. 


(81) 


where  the  v!  are  the  stoichiometric  coefficients  of  the  reactants. 

1 

the  vV  represent  the  corresponding  coefficients  for  the  reaction 

products,  and  the  represent  the  various  chemical  species.  The 

v!^  and  vV  are  either  integers  or  zero,  the  latter  indicating  that 

a  species  does  not  occur  as  a  reactant  or  product.  At  any  instant, 
the  time  rate  of  change  of  the  concentration  of  species  i  in  reaction 
j  is  given  by 


(82) 


Summing  over  all  reactions  produces  the  total  rate  of  change  of 
concentration  of  species  i: 


y 


i 


j 


(83) 
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The  next  step  in  the  Moretti  procedure  is  to  linearize  the 
equations  for  y^.  The  following  assumptions  are  made: 

1)  The  reaction  rate  coefficients  k^  and  are  constant 

throughout  a  step  (i.e.,  the  change  in  temperature  through¬ 
out  a  step  is  not  so  large  as  to  affect  their  values); 

2)  The  sum  of  concentrations  of  all  species,  Y,  is  constant; 

3)  Terms  involving  the  product  of  species  concentrations, 

y  y  ,  can  be  replaced  by  the  approximation 
ran 

V«*  "Vn  f  V.  +  Vn 

where  the  superscript  o  denotes  values  at  the  beginning 
of  a  step. 

Substitution  of  the  expression  from  assumption  3)  results  in  the 
following  set  of  linear,  nonhomogeneous,  first  order  differential 
equations : 


k' 


h  '  £  yt  +  ci 


£-1 


(85) 


where  k'  is  the  number  of  chemically  active  species. 

Moretti  solves  these  linear  equations  by  obtaining  the  eigen¬ 
values  of  the  matrix  a  .  This  procedure  becomes  unwieldy  for 

chemical  systems  involving  many  species.  DeGroat  and  Abbett  have 
developed  an  alternative  technique  based  on  the  assumption  of  a 
truncated  power  series  solution  of  the  form 

2 

*i  *  E  di£  t£  <86> 

£-0 
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Substitution  of  this  solution  into  the  general  differential  equa¬ 
tion  yields  residue  terms: 


2 

yc>  ■  E 


E  dJo  - 


j-1 


(87) 


This  set  of  equations  contains  two  unknowns,  d^  and  for  each 

species  i.  Two  conditions  are  needed  to  solve  for  these  unknowns. 
These  may  be  obtained  by  specifying  that 


(88) 


Evaluation  of  these  integrals,  employing  the  expression  for  R^(t), 

leads  to  a  system  of  linear,  nonhomogeneous,  algebraic  equations 
(10  unknowns  in  10  equations).  These  algebraic  equations  are  then 
solved  for  the  di  by  triangularizing  the  coefficient  matrix.  The 

d^.  being  known,  an  explicit  solution  is  given  for  the  molar  con¬ 
centration  of  each  species  at  the  end  of  the  time  interval  con¬ 
sidered. 


c.  Computer  Program  -  The  computer  program  developed  for 
analyzing  hydrogen/ fluorine  mixing  layers  can  be  divided  into 
three  parts.  The  first  of  these  involves  the  main  program  and 
several  subroutines  used  for  performing  various  initialization 
tasks  (see  Fig.  58).  The  other  two  parts  are  designed  to  compute 
various  flow  properties  and  solve  the  finite  difference  equations 
and  are  driven  by  subroutines  0NE  and  TW0,  respectively.  The 
finite  difference  solution  procedure  is  composed  of  two  parts, 
corresponding  to  the  gas  dynamic  and  chemical  analyses  previously 
described.  The  necessary  operations  are  designed  to  take  place 
in  sequence,  and  the  procedure  marches  downstream  in  small  steps 
as  the  flow  solution  in  the  mixing  region  is  determined. 
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PRINT 

Figure  58.  Interrelationship  of  Subroutines  in  the 
Mixing  Analysis  Computer  Code 
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For  convenience,  the  computer  code  is  divided  into  several 
subroutines.  These  subroutines  and  their  interrelationship  are 
shown  in  Figure  58.  The  function  of  each  subroutine  is  described 
in  the  paragraphs  that  follow. 

(1)  Main  Program  -  This  is  the  driver  for  the  entire  pro¬ 
gram.  The  input  data  are  read  into  the  computer  through  this 
routine  and  supplied  to  the  various  subroutines  by  means  of  C0MM0N 
statements  and  parameters  in  the  CALL  statements.  In  addition  to 
various  quantities  used  for  program  control,  the  required  input 
data  include  the  turbulent  Lewis  and  Prandtl  numbers,  the  turbu¬ 
lent  mixing  parameter,  static  pressure,  and  initial  profiles  of 
velocity,  temperature,  and  species  mass  fractions.  This  program 
also  writes  the  input  data  and,  in  the  case  of  a  power  profile 
boundary  layer,  computes  the  initial  velocity  and  temperature 
profiles.  After  completion  of  these  initialization  tasks,  the 
main  program  acts  as  the  driver  by  calling  subroutines  0NE  and 

TW0  alternately  until  the  designated  end  of  the  problem  is  reached 
or  until  an  error  occurs.  In  either  case,  control  is  returned  to 
the  beginning  of  the  program  so  as  to  run  a  following  case  if  data 
for  another  problem  are  provided. 

(2)  Subroutine  BEGIN  -  This  subroutine  supplies  informa¬ 
tion  of  a  basic  thermodynamic  nature.  Specifically,  it  provides: 

1)  Molecular  weights; 

2)  Coefficients  for  calculating  the  specific  heat 
at  constant  pressure  of  each  chemical  species; 

3)  Coefficients  for  calculating  the  enthalpy  of 
each  species. 

The  coefficients  for  specific  heat  and  enthalpy  result  from  curve 
fits  of  these  properties  as  functions  of  temperature  and  were 
derived  from  JANAF  data  (Ref.  47).  Presently,  one  set  of  coef¬ 
ficients  is  used  from  200  to  1000°K  and  a  second  set  from  1000  to 

3000°K. 


(3)  Subroutine  CHEM  -  CHEM  is  the  main  chemistry  subroutine 
and  acts  as  a  driver  for  this  phase  of  the  solution.  The  chemistry 
analysis  is  set  up  numerically  so  there  may  be  several  steps  of  the 
chemistry  solution  in  the  time  interval  corresponding  to  a  single 
fluid  dynamic  step.  CHEM  controls  the  chemistry  time  step  and  calls 
the  other  chemistry  subroutines,  REACT  and  JIM,  as  many  times  as 
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needed  to  reach  the  position  (or  time)  dictated  b"  the  fluid  dy¬ 
namic  step.  It  also  uses  the  species  mass  fractions  at  the  end 
of  each  step  to  compute  the  enthalpy  of  the  gaseous  mixture  and 
its  temperature. 

(4)  Subroutine  CLEM  -  The  program  solves  a  set  of  linear, 
nonhomogeneous ,  algebraic  equations.  The  coefficient  matrix  and 
the  constant  term  column  matrix  are  grouped  into  a  single  matrix, 
which  is  then  triangularized  by  reducing  all  the  elements  on  one 
side  of  the  diagonal  to  zero.  At  the  end,  one  of  the  unknowns  can 
be  obtained  directly  and  the  others  are  then  calculated  by  sub¬ 
stitution. 

(5)  Subroutine  DENSE  -  The  mixture  molecular  weight  is 
computed  from  the  species  mass  fractions  and  the  density  is  ob¬ 
tained  from  the  perfect  gas  equation  of  state. 

(6)  Subroutine  FSEQL  -  To  solve  the  mixing  layer  problem 
using  the  equilibrium  chemistry  option,  it  is  necessary  to  begin 
with  equilibrium  mixtures  in  the  two  freestreams.  Subroutine 
FSEQL  is  called  by  MAIN  as  part  of  the  initialization  process  to 
compute  the  required  equilibrium  mass  fractions  (from  subroutine 
HFEQL)  and  temperature  on  each  side  of  the  dividing  surface. 

(7)  Subroutine  HEAT  -  This  subroutine  uses  the  coefficients 
provided  in  subroutine  BEGIN  to  calculate  the  constant  pressure 
specific  heats  and  enthalpies  of  each  species.  Also  obtained  are 
the  mixture  specific  heat  and  enthalpy,  as  defined  by  the  mixture 
rules  previously  stated. 

(8)  Subroutine  HFEQL  -  HFEQL  is  a  special  purpose  routine 
designed  to  compute  the  species  composition  of  an  equilibrium  mix¬ 
ture  of  H,  F,  H2,  F2 ,  and  HF.  Temperature  and  density  of  the  mix¬ 
ture  are  assumed  known  for  this  solution.  Portions  of  HFEQL  are 
also  used  by  FSEQL  in  its  determination  of  freestream  composition 
and  by  REACT  to  obtain  equilibrium  constants. 

(9)  Subroutine  INGRT2  -  The  boundary  layer  input  option 
is  based  on  prescribing  the  behavior  of  the  velocity  and  temper¬ 
ature  profiles  (presently  these  are  simple  power  law  profiles) 
over  a  distance  specified  by  the  boundary  layer  thickness.  Since 
the  initial  profiles  must  be  provided  in  the  stream  function  plane, 
this  subroutine  is  used  to  perform  the  transformation  from  boundary 
layer  profiles  in  the  y  coordinate  system  to  profiles  in  the  ip 


126 


coordinate  system.  A  quadratic  interpolation  scheme  is  employed 
to  obtain  an  evenly  spaced  v  mesh  from  the  unevenly  spaced  mesh 
resulting  from  this  transformation. 

(10)  Subroutine  JIM  -  The  expressions  to  be  solved  for 

the  unknown  d.,  and  d . „  coefficients  form  a  system  of  10  linear, 
ll  i2 

nonhomogeneous,  algebraic  relations  in  10  unknowns.  The  coef¬ 
ficients  of  these  10  unknowns,  as  well  as  the  constant  terms  in 
the  equations,  are  computed  in  this  subroutine.  The  program  then 
calls  subroutine  CLEM  to  solve  the  resulting  matrix  and  uses  the 
results  to  calculate  the  post-reaction  mass  fractions. 

(11)  Subroutine  g)NE  -  This  program  calls  subroutines 
PRESS,  DENSE,  and  HEAT  to  calculate  the  pressure,  density,  specific 
heat,  and  enthalpy.  It  then  calculates  the  values  of  the  physical 
transverse  coordinate,  y,  from  the  values  of  the  stream  function, 

The  viscosity  is  obtained  by  calling  subroutine  VISC.  The 
program  keeps  track  of  the  number  of  axial  stations  and  checks  to 
see  if  the  desired  total  distance  downstream  has  been  reached. 

If  this  distance  has  been  covered,  subroutine  PRINT  is  called  to 
write  the  final  results. 

(12)  Subroutine  PRESS  -  The  input  data  provided  in  the 
main  program  contain  pressure  regions  and  a  pressure  polynomial 
for  each  region.  These  data  are  used  by  this  program  to  compute 
the  pressure  and  axial  pressure  gradient  for  a  given  value  of  the 
x  coordinate. 

(13)  Subroutine  PRINT  -  The  purpose  of  this  subroutine 
is  to  write  the  results  obtained  for  values  of  various  flow-field 
variables.  It  is  called  periodically  from  subroutine  TW0,  as  con¬ 
trolled  by  a  printing  interval  specified  by  input,  and  from  sub¬ 
routine  0NE  when  the  final  step  is  completed.  Variables  whose 
values  are  printed  include  the  y  and  ip  coordinates,  temperature, 
enthalpy,  density,  velocity,  fringe  shift,  Mach  number,  and  species 
mass  fractions  on  each  streamline  in  the  mixing  region. 

(14)  Subroutine  REACT  -  This  subroutine  calculates  the 
forward  and  backward  reaction  rate  constants  for  a  given  temper¬ 
ature.  Subroutine  HFEQL  is  called  to  provide  the  equilibrium  con¬ 
stants  needed  for  this  calculation.  Using  the  reaction  rate  con¬ 
stants,  the  coefficients  of  a  and  the  constants  C,  of  the  pre- 

viously  developed  differential  equations  for  y^  are  computed  for 
use  by  subroutine  JIM. 
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(15)  Subroutine  TW0  -  The  first  step  involved  in  this  sub¬ 
routine  is  the  establishment  of  an  axial  step  size  that  satisfies 
the  numerical  stability  requirements  of  the  chosen  explicit  finite 
difference  scheme.  The  program  has  options  for  frozen,  equilibrium, 
or  finite-rate  chemistry.  In  the  case  of  frozen  flow,  the  boundary 
layer  equations  can  be  solved  without  evaluating  the  species 
production  terms.  Equilibrium  and  finite  rate  problems  call  sub¬ 
routines  WD0TI2  and  WD0TI,  respectively,  to  obtain  the  updated 
information  regarding  chemical  species  concentrations  at  the  end 

of  a  given  interval.  Once  this  is  done,  the  program  is  ready  to 
solve  the  mixing  layer  equations.  The  program  tests  for  errors, 
with  particular  regard  for  negative  species  concentrations.  If 
an  error  condition  is  detected,  the  problem  will  ultimately  be 
terminated  after  printing  the  results.  Otherwise,  the  program 
sets  up  tests  to  determine  whether  there  is  any  need  for  increas¬ 
ing  the  width  of  the  mixing  zone.  These  tests  check  the  velocity, 
temperature,  and  species  mass  fraction  profiles  for  zero  slope  at 
the  upper  and  lower  edges.  Because  these  profiles  approach  free- 
stream  values  asymptotically,  a  small  finite  value  rather  than 
absolute  zero  is  used  as  a  minimum  slope  requirement.  The  program 
automatically  addc  streamlines  as  needed  as  it  marches  downstream. 

To  control  the  total  number  of  streamlines  being  carried  in  the 
mesh,  there  is  a  provision  for  doubling  the  mesh  size  by  discard¬ 
ing  every  other  streamline.  The  permissible  maximum  number  of 
streamlines  is  controlled  by  an  input  variable. 

(16)  Subroutine  VISC  -  This  subroutine  computes  the  tur¬ 
bulent  eddy  viscosity  according  to  the  model  previously  discussed. 
Options  for  other  viscosity  models  can  be  incorporated  in  this 
subprogram  or  another  specialized  version  can  readily  be  substi¬ 
tuted. 


(17)  Subroutine  WD0TI  -  The  initial  operations  carried 
out  in  this  program  involve  tests  to  see  if  subroutine  TW0  has 
either  doubled  the  mesh  size  or  added  streamlines  to  the  mixing 
zone  since  the  last  chemistry  calculation.  If  one  of  these  pro¬ 
cedures  has  occurred,  appropriate  adjustments  are  made  to  several 
parameters  affecting  the  chemistry  computations.  The  program  then 
computes  the  diffusion  time  step  (Ax/u)  on  each  streamline  and 
calls  subroutine  CHEM  to  calculate  the  new  temperature  and  species 
mass  fractions  at  the  end  of  this  time  step. 
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(18)  Subroutine  WD0TI2  -  This  program  performs  an  iter¬ 
ative  determination  of  the  equilibrium  temperature,  density,  and 
species  composition  at  the  end  of  each  fluid  dynamic  mixing  step. 
The  values  of  these  properties  resulting  from  the  mixing  process 
are  determined  in  subroutine  TW0,  and  are  modified  here  to  reflect 
a  new  equilibrium  state.  Subroutine  HFEQL  is  called  to  compute 
mass  fractions  for  given  temperature  and  density.  Because  the 
final  temperature  depends  on  enthalpy  (and,  thereby,  on  species 
mass  fractions)  and  the  final  density  depends  on  the  temperature 
and  mixture  molecular  weight,  an  iterative  process  is  required 
to  produce  a  solution  for  all  of  these  variables. 

3.  RESULTS  OF  ANALYSIS  OF  PARALLEL  TURBULENT  MIXING  LAYERS 

Three  cases  were  investigated  using  the  mixing  analysis  pro¬ 
gram.  The  data  for  these  cases  represent  two  types  of  experi¬ 
mental  runs  in  the  parallel  mixing  facility.  In  one  run  a  50% 
mixture  of  fluorine  and  helium  (by  volume)  was  used  for  one  stream, 
while  a  50%  mixture  of  hydrogen  and  argon  formed  the  other 
stream.  For  the  other  run  the  argon  was  eliminated  on  the  hydrogen 
side.  Both  of  these  situations  were  investigated  analytically 
using  finite-rate  chemistry,  and  the  former  run  was  also  com¬ 
puted  using  the  frozen  chemistry  program  option.  Initial  free- 
stream  conditions  used  for  these  calculations  are  provided  in 
Table  XI.  In  all  cases,  the  wall  temperature  of  the  splitter 
plate  was  assumed  to  be  295°K  and  the  boundary  layers  were  com¬ 
puted  using  1/7  as  the  power  for  the  power  law  profiles.  The 
mixing  parameter,  ' ,  was  taken  as  20  for  the  reacting  flows  and 
37  for  the  frozen  case. 

a.  Chemically  Frozen  Flow  with  Argon  -  Consider  first  the 
results  shown  in  Figures  59  thru  61  for  the  chemically  frozen 
flow  with  argon  present  (Run  72091301; .  Frozen  chemistry  does 
not  represent  a  realistic  situation  for  hydrogen/fluorine  mixing 
flows.  However,  it  is  useful  to  illustrate  the  results  of  mix¬ 
ing  occurring  without  chemical  reaction  before  proceeding  to  the 
reacting  case.  Note  that  the  ordinate  scale  of  these  figures  is 
greatly  expanded  relative  to  that  used  in  Figures  62  thru  67. 
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Table  XI 


INITIAL  CONDITIONS  FOR  MIXING  LAYER  ANALYSES 


Run  No. 

Species  Mass 
Fraction 

p, 

dynes/cm2 

u, 

cm/sec 

T, 

°K 

6bl, 

cm 

72091301 

F  0.4162 

F  0.4886 

He  0.0952 

6.306  x  104 

2.050  x  105 

264 

0.0825 

H2  0.0478 
Ar  0.9522 

6.306  x  104 

2.570  x  105 

450 

0.1375 

72111602 

F  0.417 

F2  0.488 

He  0.095 

4.099  x  104 

1.990  x  10s 

264 

0.0765 

h2  i.ooo 

4.099  x  104 

4.24  x  105 

159 

0.1275 

Four  temperature  profiles  are  shown  in  Figure  59  to  illustrate 
the  development.  The  initial  profile  indicates  the  effect  of  a  wall 
that  is  cold  relative  to  one  stream  and  hot  relative  to  the  other. 

As  the  flows  mix,  turbulent  viscous  heating  occurs,  causing  the 
temperature  to  rise  in  the  regions  of  maximum  shear.  Large  gra¬ 
dients  in  the  velocity  profile  due  to  boundary  layers  cause  highest 
local  shear  levels  in  profiles  closest  to  the  splitter  plate.  The 
temperature  profile  at  0.502  cm  shows  residual  effects  of  the  zero 
velocity  gradient  (and  consequent  zero  shear)  at  the  minimum  of 
the  velocity  profile.  As  mixing  continues,  the  velocity  gradients 
decrease  in  magnitude  and  shear  becomes  less  important.  At  the 
profile  position  4.466  cm  downstream,  the  temperature  peak  has 
decreased  to  about  300°K  as  the  mixing  process  dominates. 

Figure  60  shows  profiles  of  fringe  shift  computed  as  discussed 
in  Section  II.  These  profiles  reflect  the  local  density  varia¬ 
tion,  as  well  as  the  species  concentrations.  In  the  constant 
pressure  situation  that  holds  here,  the  maximum  in  the  initial 
profile  is  related  to  the  wall  temperature,  which  is  colder  than 
freestream  on  the  hydrogen  side.  As  mixing  proceeds,  the  density 
decreases  as  the  temperature  increases.  This  causes  the  maximum 
in  the  fringe-shift  profile  to  disappear  and  a  minimum  is  formed 
in  its  place.  Note  that  the  depth  of  this  minimum  is  somewhat 
smaller  at  the  4.466  cm  position  than  at  the  preceding  station. 
Apparently,  the  mixing  process  is  now  removing  thermal  energy 
from  the  central  portion  of  the  mixing  layer  at  a  faster  rate  than 
it  is  being  generated.  This  phenomenon  can  also  be  observed  in 
the  temperature  distributions  of  Figure  59. 
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Figure  59.  Computed  Temperature  Profiles  for  Run  72091301,  Frozen  Chemistry 
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Figure  60.  Computed  Fringe-Shift  Profiles  for  Run  72091301,  Frozen  Chemistry 
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Figure  61.  Computed  H2  and  F  Species  Profiles  for  Run  72091301, 
Frozen  Chemistry 
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Profiles  for  Run  72091301,  Nonequilibrium  Chemistry 
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Figure  65.  Computed  Temperature  Profiles  for  Run  72111602,  Nonequilibrium  Chemistry 
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Figure  66.  Computed  Fringe-Shift  Profiles  for  Run  72111602 »  Nonequilibrium  Chemistry 
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Figure  67.  Computed  HF  Species  Profiles  for  Run  72111602,  Nonequilibrium  Chemistry 


Figure  61  presents  the  mass  fraction  profiles  for  the  F  and  H2 
species  at  the  initial  station  and  4.466  cm  downstream.  These 
profiles  are  typical  of  the  shape  attained  when  no  chemical  re¬ 
actions  are  allowed  and  show  a  smooth  transition  from  freestream 
values  on  one  side  of  the  mixing  zone  to  zero  on  the  opposite  side. 

b.  Reacting  Flow  with  Argon  -  The  second  computation  is  repre¬ 
sented  by  Figures  62  thru  64.  The  initial  data  are  the  same  (Run 
72091301)  as  the  preceding  case,  but  finite-rate  chemical  reactions 
are  now  considered  in  the  analysis.  The  processes  that  occur  with 
frozen  chemistry  still  take  place  in  this  situation,  but  the  tem¬ 
perature  profiles  of  Figure  62  illustrate  that  they  are  greatly 
overshadowed  by  the  results  of  combustion  energy  release.  This 
energy  release  is  confined  to  a  relatively  narrow  region  of  the 
flow.  The  mixing  region  at  the  6.140  cm  station  has  a  total  width 
of  about  3  cm,  with  only  about  1  cm  showing  significant  temper¬ 
ature  rise.  This  result  is  influenced  by  the  chosen  eddy  viscosity 
model  and  the  value  of  the  associated  mixing  parameter,  a.  Smaller 
values  of  this  parameter  would  cause  higher  mixing  rates  and  some¬ 
what  faster  spreading  of  the  reaction  zone.  Our  experience  and 
the  data  from  Figure  56  indicate  that  the  value  o  *  20  is  appro¬ 
priate  for  reacting  flows  in  the  present  velocity  range.  However, 
until  substantiation  is  obtained  through  satisfactory  correlation 
with  experimental  data,  this  must  be  considered  as  an  assumption 
for  analysis  of  the  hydrogen/fluorine  mixing  process. 

The  fringe-shift  distributions  of  Figure  63  react  primarily 
to  the  density  and  show  the  expected  development  of  a  deep  mini¬ 
mum  in  the  high-temperature  region.  Because  the  temperature  rise 
between  the  stations  at  2.038  and  6.140  cm  is  not  large,  no  great 
change  is  found  in  the  value  of  the  minimum  in  this  region.  The 
initial  profile  is  identical  to  that  seen  in  the  frozen-chemistry 
case,  although  plotted  on  a  different  scale  than  in  Figure  60. 

Figure  64  shows  how  the  combustion  process  contributes  to  the 
potential  for  lasing  activity,  as  measured  by  the  fraction  of  the 
flow  that  forms  hydrogen  fluoride.  Because  no  HF  exists  until 
mixing  commences,  no  initial  profile  is  required  here.  As  would 
be  expected  from  the  associated  temperature  distributions,  the 
region  of  high  HF  mass  fraction  is  fairly  narrow  and  the  peak 
value  increases  slightly  as  the  flow  moves  downstream.  A  com¬ 
parison  of  mixing  layer  growth  between  this  case  and  the  preced¬ 
ing  frozen-chemistry  case  shows  that  the  influence  of  combustion 
is  dominant.  Chemical  processes  control  the  turbulence  level  and 
fluid  mechanical  effects,  e.g.,  velocity  differences  are  secondary. 
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As  a  result,  the  HF  species  mass  fraction  can  be  considered  as  a 
primary  measure  of  the  mixing  region  width,  while  the  growth  as 
indicated  by  fluid  mechanical  properties  is  constrained  to  follow 
the  trend  established  by  the  chemistry.  This  produces  nearly  iden¬ 
tical  mixing  layer  widths,  independent  of  the  property  chosen  as 
an  indicator. 

c.  Reacting  Flow  without  Argon  -  Figure  65  illustrates  the 
streamwise  development  of  the  temperature  distribution  for  the 
case  in  which  pure  hydrogen  forms  one  stream.  The  density  in 
this  stream  is  approximately  1/5  of  the  density  of  the  correspond¬ 
ing  flow  with  argon.  Additionally,  the  wall  is  now  "hot"  with 
*  ;spect  to  both  frees tream  temperatures.  Note  that  the  peak  of 
the  temperature  profile  moves  more  rapidly  toward  the  fluorine 
side  of  the  flow  than  was  the  case  in  Figure  62.  Abramovich  (Ref. 
48)  has  shown  that  a  mixing  layer  spreads  preferentially  into  the 
lower  velocity  stream  when  both  flows  have  the  same  density.  In 
the  present  cases  (both  Run  72091301  and  Run  7211160 2)  the  density 
on  the  fluorine  side  is  higher  than  on  the  hydrogen  side,  but  the 
same  trend  with  velocity  is  observed.  However,  comparison  of 
Figures  59  and  62  (which  differ  only  in  terms  of  combustion  ef¬ 
fects)  indicates  that  preferential  spreading  into  the  fluorine 
stream  occurs  more  rapidly  with  combustion  than  with  the  chem¬ 
ically  frozen  flow.  It  appears,  therefore,  that  other  parameters 
in  addition  to  velocity  are  important  in  determining  asymetrical 
growth  of  the  mixing  region.  Understanding  of  the  pertinent 
mechanisms  will  require  further  study. 

Figures  66  and  67  show  the  profiles  of  fringe  shift  and  HF 
mass  fraction  for  the  pure  hydrogen  case.  The  main  features  of 
these  distributions  are  reflections  of  corresponding  features  in 
the  temperature  profiles  discussed  in  the  preceding  paragraph. 

Of  course,  the  hydrogen-fluoride  mass  fraction  is  much  larger  than 
it  was  when  argon  was  present  as  a  diluent.  Note  also  that  the 
mixing  layer  at  any  station  is  substantially  wider  now  than  was 
the  case  previously.  This  is  largely  due  to  the  higher  velocity 
of  the  pure  hydrogen  stream. 


4.  EXTENSION  OF  ANALYSIS  TO  INJECTION  MIXING 


The  flow  conditions  in  the  immediate  region  of  injection  are 
quite  complex.  Figure  68  depicts  some  of  the  features  that  can 
be  expected  to  exist.  Because  the  splitter  plate  height  is  large 
in  comparison  to  the  injector  slot  height,  a  low-velocity  circu¬ 
lation  flow  will  be  present  in  the  base  region  of  the  plate.  As 
indicated,  the  circulation  will  consist  of  two  cells  which  mix 
gases  from  the  plume  and  external  flow,  with  the  possibility  of 
burning  at  low  velocities.  It  is  downstream  of  this  region  that 
the  two  supersonic  streams  impinge  on  one  another  and  supersonic 
turbulent  mixing  occurs.  The  flow  deflections  associated  with 
the  impingement  of  the  two  streams  introduce  a  shock  in  the  ex¬ 
ternal  flow,  along  with  a  new  shock  in  the  plume  flow.  The  new 
shock  reinforces  the  boundary  shock  that  already  exists  in  the 
plume  flow. 

An  analysis  of  the  flow  situation  depicted  in  Figure  68  has 
been  performed.  Plume  flow  fields  were  generated  for  conditions 
of  quiescent  surroundings  and  also  Mach  5  external  flow,  both 
with  an  ambient  pressure  of  30  torr.  Although  neither  of  these 
two  boundary  conditions  properly  describe  the  complex  condition 
that  exists,  the  results  do  provide  an  approximation  to  the  range 
of  pressure  and  velocity  conditions  that  exist  along  the  boundary. 
The  intent  of  this  analysis  was  to  determine  if  the  Lockheed 
method  of  characteristics  plume  program  (Ref.  49)  could  be  used 
to  establish  the  pressure  and  velocity  distributions.  These  dis¬ 
tributions  could  then  be  used  to  determine  the  characteristics  of 
the  mixing  and  combustion  region. 

Some  of  the  results  of  the  analysis  are  shown  in  Figure  69, 
where  plume  boundaries  and  internal  shock  structure  are  plotted 
for  the  two  external  flow  conditions.  In  both  cases  the  pressure 
and  Mach  numbers  along  the  boundary  are  nearly  constant.  For  the 
quiescent  conditions,  the  boundary  Mach  number  is  5.02  at  a  pres¬ 
sure  of  0.0394  atm  and  for  the  Mach  5  freestream  condition,  the 
boundary  Mach  number  is  2.62  at  a  pressure  of  0.9987  atm. 

Because  of  the  large  variation  in  pressure  between  these  two 
cases  and  the  strong  influence  of  external  flow,  it  was  decided 
that  the  analysis  must  be  refined  further  to  obtain  a  realistic 
description  of  the  boundary  conditions  between  the  two  streams. 

This  refinement  would  involve  making  minor  modifications  to  the 
method  of  characteristics  program  to  provide  a  capability  of  treat¬ 
ing  combined  boundary  conditions,  i.e.,  quiescent  conditions  that 
simulate  the  circulation  region  followed  by  a  supersonic  external 
flow  of  varying  Mach  number.  Unfortunately,  this  modification 
exceeds  the  scope  of  the  present  program. 
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Figure  68.  Injection  Mixing  Flow  Field 
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Figure  69.  Plume  Analysis  Results 
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Once  the  plume  flow  field  has  been  adequately  established,  the 
initial  approach  to  the  mixing  analysis  will  involve  superposition 
of  a  mixing  layer  on  the  plume  boundary.  The  required  profiles  of 
velocity  and  temperature  will  be  obtained  from  the  plume  and  free- 
stream  at  the  point  of  impingement.  Flow  direction  at  this  point 
should  apprcximately  satisfy  the  mixing  program  assumption  of  paral¬ 
lel  flow.  Of  course,  this  analysis  is  appropriate  only  while  the 
mixing  layer  remains  small  relative  to  the  plume,  so  that  its  in¬ 
fluence  on  the  overall  flow  is  also  small.  How  well  this  condition 
will  be  satisfied  can  only  be  determined  by  the  actual  computations 
and  remains  to  be  resolved. 

Further  refinement  of  the  injection  mixing  analysis  will  prob¬ 
ably  involve  reformulation  of  the  basic  equations  to  include  pro¬ 
vision  for  lateral  pres Jure  gradients.  The  effort  required  to 
accomplish  this  change  is  necessarily  much  greater  than  that  in¬ 
volved  in  implementing  the  relatively  straightforward  superposi¬ 
tion  approach.  Justification  for  such  an  effort  must  depend 
largely  on  the  success  or  failure  of  the  latter  procedure  in  cor¬ 
relating  the  injection  mixing  data. 
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SECTION  V 


DISCUSSION  OF  RESULTS 


As  earlier  discussions  have  indicated,  the  experimental  objec¬ 
tives  of  this  investigation  were  pursued  in  two  different  facili¬ 
ties  using  variations  of  a  basic  technique  developed  at  Martin 
Marietta  for  observing  the  mixing  and  chemical  reaction  of  hydrogen/ 
oxygen  streams.  The  primary  diagnostic  equipment  employed  with  each 
of  these  facilities  consisted  of  a  Mach-Zehnder  interferometer  for 
mapping  averaged  density  profiles  across  the  mixing  layer,  a  rapid 
scanning  infrared  spectrometer  for  mapping  radiative  emission  within 
the  mixing  layer,  pitot  probes  for  surveying  pressure  across  both 
the  mixing  layer  and  the  nozzle  exit  region,  and  static  side-wall 
pressure  taps  for  determining  nozzle  exit  static  pressure.  In  ad¬ 
dition,  attempts  were  made  to  measure  the  nozzle  exit  temperature 
with  an  unshielded  thermocouple  probe. 

The  major  emphasis  in  the  experimental  phase  of  this  study 
was  directed  towards  the  injection  mixing  experiments.  The  injec¬ 
tion  facility  also  was  used  to  perfect  the  diagnostic  techniques 
and  test  procedures  that  were  used  later  with  the  parallel  mixing 
•facility.  The  primary  emphasis  in  the  analytical  phase  of  the  study 
was  devoted  to  refining  and  modifying  an  existing  analytical  model 
of  the  reacting  turbulent  shear  layer.  The  remainder  of  this  section 
will  be  devoted  to  a  discussion  of  (1)  the  experimental  results  ob¬ 
tained  in  both  the  injection  and  the  parallel  mixing  facilities, 
and  (2)  a  comparison  of  the  theoretical  results  with  the  experimen¬ 
tal  results  obtained  in  the  parallel  mixing  experiments. 

1.  EXPERIMENTAL  RESULTS 

a.  Injection  Mixing 

Experimental  observations  within  the  test  section  region  of  the 
hydrogen  injection  facility  have  provided  quantitative  information 
about  the  overall  characteristics  of  the  mixing  layer.  In  partic¬ 
ular,  the  growth  rate,  radiative  emission  intensity,  and  averaged 
density  profiles  have  been  obtained  under  a  variety  of  test  section 
conditions  for  two  different  injection  splitter  plate  geometries. 

The  snapshot  interferograms  presented  in  Figures  30  and  31  of 
Section  III  show  the  overall  flow  characteristics  for  two  different 
injection  geometries  with  five  different  injection  conditions, 
including  a  no-reaction  helium  injection  case  and  a  no-injection 
nozzle  flow  case. 
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A  comparison  of  the  fringe-shift  profiles  between  the  150-psi  heli¬ 
um  injection  interferogram  and  the  150-psi  hydrogen  injection  inter-* 
ferogram  definitely  indicates  that  the  chemical  reaction  is  respon¬ 
sible  for  the  considerably  greater  fringe  distortion  that  occurred 
within  the  mixing  layer  of  the  latter  case.  The  other  general  con¬ 
clusion  that  can  be  drawn  from  a  simple  observation  of  these  inter¬ 
ferometric  records  is  that  the  downstream  extent  of  the  mixing  layer 
is  greater  for  the  150-psi  injection  cases  than  it  is  for  the  lower 
injection  pressures.  As  explained  in  detail  in  Section  III,  a  meas¬ 
ure  of  the  mixing  layer  width  can  be  obtained  from  the  interfero- 
grams.  The  results  presented  in  Figure  36  indicate  that,  for  the 
orifice  injector  geometry,  the  mixing  layer  width  increases  with 
increasing  injection  pressure.  A  decrease  in  mixing  layer  width 
is  recorded  for  the  slotted  injector  geometry,  especially  for  the 
10-torr  nominal  test  section  static  pressure.  In  both  cases,  the 
mixing  layer  widths  for  the  50-psi  injection  pressure  are  almost 
identical  for  both  the  10-torr  and  the  15-torr  nominal  background 
pressures . 

Perhaps  the  most  interesting  aspect  of  the  comparison  between 
the  slotted  and  the  orifice  injection  cases  is  shown  in  the  aver¬ 
aged  density  profiles  presented  in  Figure  38.  Several  assumptions 
were  required  in  arriving  at  these  averaged  density  profiles,  the 
most  significant  being  that  the  minimum  fringe  shift  was  identi¬ 
fied  as  due  primarily  to  HF  speciea.  This  assumption,  although 
seemingly  arbitrary,  is  based  upon  the  theoretical  results  obtained 
in  the  parallel  shear  flow  mixing  analysis.  As  indicated  in  Fig¬ 
ure  62  of  Section  IV,  the  peak  HF  species  number  density  always 
corresponds  to  the  minimum  fringe-shift  position.  The  averaged 
number  density  profiles  presented  in  Figure  38  show  the  gradual 
decrease  in  number  density  moving  downstream  along  the  plane  of 
injection.  A  comparison  between  results  for  the  slotted  hydrogen 
injection  and  the  orifice  hydrogen  injection  shows  higher  concen¬ 
tration  (by  a  factor  of  2)  of  hydrogen  along  the  centerline  plane 
for  the  slotted  injection  scheme.  This  is  a  direct  result  of  the 
higher  mass  flow  through  the  larger  exit  area  of  the  slotted  in¬ 
jector.  It  is  also  very  significant  to  note  that  the  indicated 
HF  number  densities  (corresponding  to  minimum  fringe-shift  posi¬ 
tion)  for  the  150-psi  and  the  50-psi  slotted  hydrogen  injection 
are  almost  identical.  This  observation  is  in  agreement  with  the 
measured  infrared  emission  intensities  shown  in  Figure  47.  That 
is,  the  specific  radiative  intensity  at  any  downstream  position 
does  not  change  with  increasing  injection  pressure,  at  least  be¬ 
tween  50  and  150  psi.  This  indicates  that  even  at  50  psi  the  mass 
flow  of  hydrogen  is  sufficiently  large  that  the  reaction  is  mixing 
controlled.  Hence,  changing  injection  pressure  above  50  psi  will 
not  change  the  degree  of  mixing  on  the  molecular  scale,  although 
it  does  change  the  geometry  of  the  mixing  region. 
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The  results  of  the  pitot  probe  survey  served  to  indicate  that 
the  flow  fields  were  indeed  uniform  within  the  nozzle  exit  region 
and  across  the  width  of  the  test  section.  A  limited  effort  was 
devoted  to  use  of  the  temperature  probe  because  of  late  delivery 
of  the  necessary  parts.  However,  results  were  obtained  that  very 
clearly  indicated  that  a  shielded,  vented  one-mil  thermocouple  probe 
using  a  tungsten/tungsten-rhenium  junction  will  adequately  measure 
the  temperature  within  reacting  mixing  layers  for  shock  tunnel  flows. 

b.  Parallel  Mixing 

Unlike  the  injection  mixing  case,  the  parallel  mixing  snapshot 
interf erograms  do  not  provide  a  significant  level  of  information 
about  the  mixing  zone  until  the  data  are  completely  reduced  and 
compared  with  theoretical  predictions.  Typical  interf erograms 
of  the  parallel  mixing  layers  are  shown  in  Figure  39  and  the  re¬ 
duced  fringe-shift  results  at  various  downstream  stations  are  shown 
in  Figures  40  and  41  of  Section  III.  Conditions  for  two  of  the  ex¬ 
perimental  cases  presented  in  Figure  40  have  been  used  as  input  for 
the  theoretical  model  discussed  in  Section  IV  and  a  detailed  com¬ 
parison  for  these  two  cases  will  now  be  presented  and  discussed. 

2 .  ANALYTICAL  RESULTS 

A  detailed  description  of  the  theoretical  model,  developed 
in  this  study,  for  predicting  flow-field  properties  within  a 
turbulent  shear  flow  mixing  layer  was  presented  in  Section  IV. 

As  discussed  therein,  the  criterion  for  ascertaining  the  exist¬ 
ence  of  turbulence  in  a  shear  layer  is  based  on  the  momentum 
thickness  Reynolds  number,  ReQ . 

A  problem  in  defining  the  state  of  the  mixing  layer  exists 
if  the  initial  boundary  layers  on  the  separating  surface  are 
laminar.  It  is  well  known  that  free  shear  layers  are  much  more 
unstable  than  boundary  layers.  The  critical  momentum  thickness 
Reynolds  number  for  boundary  layers  is  about  500,  while  for 
free  shear  layers  it  is  about  1/3  of  this  value.  Customarily 
the  free  shear  layer  Reynolds  number  is  defined  in  terms  of  the 
velocity  difference  across  the  layer,  |uj  -  u2 | .  However,  ex¬ 
amination  of  the  formation  of  a  free  shear  layer  from  initially 
separated  streams  reveals  a  region  of  low  energy  gas  that  orig¬ 
inated  in  the  wall  region  of  the  boundary  layer  and  continues  to 
exist  for  a  considerable  distance  downstream  in  the  free  shear 
layer.  Therefore,  the  shear  that  is  characteristic  of  the  initial 
stages  of  the  free  shear  layer  is  probably  more  closely  propor¬ 
tional  to  (uj  +  u2)/2  than  to  |uj  -  u2|.  This  reasoning,  coupled 
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with  the  fact  that  free  shear  layers  are  notoriously  unstable 
and  critically  subject  to  such  influences  as  free  stream  tur¬ 
bulence,  acoustic  environment,  and  separation  surface  roughness, 
leads  to  the  assumption  that  if  the  sum  of  the  boundary  layer 
momentum  thickness  Reynolds  numbers  is  near  the  critical  value  for 

free  shear  layers  (Re  150) ,  the  shear  layer  will  be  turbulent. 

0 

Such  an  assumption  has  led  to  good  correlations  between  theory 
and  experiment. 

Clearly,  a  transition  region  must  exist  between  the  end  of 
the  separating  surface  with  laminar  boundary  layers  and  a  fully 
formed  turbulent  free  shear  layer.  At  present,  however,  the 
sensitivity  of  free  shear  layer  stability  to  anomalous  disturb¬ 
ances  is  so  little  understood  that  it  is  impractical  to  account 
for  such  a  transition  region. 

Values  of  Re  for  each  nozzle  have  been  computed  for  the 

U 

parallel  mixing  cases  and  the  results  are  shown  in  Table  VII. 

In  all  cases,  the  summation  of  the  two  values  for  adjacent  nozzles 
does  exceed  the  critical  value  of  150  and  the  mixing  region  must 
be  considered  turbulent. 

A  theoretical  analysis  of  the  mixing  layer  for  two  different 
parallel  mixing  experiments  (Runs  72091301  and  72111602)  has  been 
completed.  Experimental  test  conditions  were  used  as  input  quanti¬ 
ties  for  the  computer  program,  along  with  computed  nozzle  wall 
boundary  layers  assuming  1/7  power  profiles.  A  value  for  the  mix¬ 
ing  layer  growth  rate  parameter,  a  (an  empirical  parameter  of  the 
eddy  viscosity  model),  was  fixed  at  20  for  these  first  calculations. 
This  is  a  very  reasonable  starting  value  for  o  based  on  earlier 
Martin  Marietta  results  that  are  also  in  agreement  with  results  ob¬ 
served  by  Shackleford  of  TRW  (Ref  40).  Figures  70  and  71  show  the 
comparisons  between  computed  fringe-shift  profiles  and  experimental 
data  points  for  the  two  sets  of  conditions.  Several  noteworthy 
features  are  evident  from  the  theoretical  results:  (1)  computed 
fringe-shift  profiles  show  the  mixing  proceeding  preferentially 
into  the  fluorine  side  of  the  stream,  as  indicated  by  the  minimum 
fringe-shift  position  (maximum  HF  number  density)  moving  upward 
into  the  fluorine  carrying  stream  with  distance  downstream;  (2) 
width  of  the  computed  mixing  region  very  nearly  matches  that  of 
the  experimental  data;  and  (3)  shapes  of  the  experimental  profiles 
are  reasonably  well  predicted  for  each  of  the  two  appreciably  dif¬ 
ferent  test  cases. 

The  experimental  data  and  the  theoretical  profiles  of  Figures 
70  and  71  at  first  observation  do  not  appear  to  be  in  good  agree- 
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ment.  The  experimental  points  shown  in  these  figures  are  in  their 
true  position  relative  to  the  abscissa  when  it  is  interpreted  as 
a  line  parallel  to  an  extension  of  the  splitter  plate.  The  theo¬ 
retical  computation,  on  the  other  hand,  assumes  coincidence  in 
the  directions  of  the  abscissa  and  the  flow  velocity  vector.  How¬ 
ever,  due  to  the  existence  of  a  small  mismatch  in  initial  stream 
pressures,  the  velocity  vector  does  not  lie  parallel  to  the  splitter 
plate,  but  rather  is  rotated  slightly  toward  the  hydrogen  stream. 

If  a  similar  rotation  is  applied  to  the  axis  of  theoretical  results, 
the  profiles  are  shifted  in  a  manner  that  places  the  minimum  fringe- 
shift  locations  much  closer  to  those  obtained  experimentally. 

Figures  72  and  73  show  the  comparisons  based  upon  this  assumption 
and  the  improved  agreement  that  results.  It  should  be  pointed  out 
that  the  first  experimental  profile  always  shows  the  greatest  in¬ 
fluence  of  the  splitter  plate  boundary  layers.  This  effect  has 
been  noted  in  earlier  work  and  can  be  accounted  for  in  the  analyt¬ 
ical  program  by  adjusting  the  initial  boundary  layer  profiles. 

A  frozen  flow  (a  =  37)  calculation  was  also  performed  using 
the  conditions  for  Run  72091301  in  order  to  compare  the  effects  of 
a  reacting  and  a  nonreacting  system.  The  results  of  this  calcula¬ 
tion  are  shown  in  the  top  half  of  Figure  70  and  it  is  clear  that 
the  chemically  reacting  results  are  very  much  different  than  those 
for  the  assumed  nonreacting  case.  The  basic  difference  between  the 
fringe-shift  profiles  computed  in  the  frozen-flow  case  (arbitrarily 
assuming  that  no  chemical  reaction  occurs)  and  the  nonequilibrium  case 
is  the  more  limited  extent  of  the  mixing  zone,  as  indicated  by  the 
widths  of  the  depressions  in  the  fringe-shift  profiles  at  the  down¬ 
stream  positions;  that  is,  the  growth  rate  for  a  chemically  reacting 
mixing  layer  is  greater  than  that  for  a  similar  nonreacting  mixing 
layer.  This  phenomenon  has  been  observed  in  previous  experiments 
in  this  laboratory  involving  mixing  layers  with  hydrogen/ oxygen 
combustion. 

Spectroscopic  data  from  the  parallel  mixing  experiments  clearly 
indicate  that  no  radiation  beyond  3.2  microns  was  present.  A  con¬ 
servative  estimate  of  the  number  of  excited  HF  molecules  obtained 
from  the  infrared  scans  indicates  that  approximately  8.5  x  1011* 
excited  HF  molecules  are  present  in  a  specified  volume  in  Run 
72111602.  Using  the  computed  HF  mass  fraction,  the  total  number 
of  HF  molecules  within  the  same  volume  was  found  to  be  2.4  x  1015. 
These  results  provide  additional  evidence  that  the  analytical  pro¬ 
gram  correctly  reflects  the  physical  T/henomena  in  the  mixing  region. 
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Flfur*  70.  hii  72091101:  fringe -Shift  Profile* 


1 


154 


SECTION  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


As  indicated  in  the  introductory  remarks,  initial  studies 
with  continuous  flow  chemical  laser  devices  at  several  labora¬ 
tories  have  demonstrated  that  a  complex  interrelationship  does 
indeed  exist  between  gas  dynamic  mixing  phenomena  and  chemical 
kinetic  reaction  processes.  An  important  objective  of  this 
study  was  to  demonstrate  the  feasibility  of  using  reasonably 
simple  experimental  techniques  to  simulate  these  complex  mixing 
phenomena  in  such  a  manner  that  a  variety  of  diagnostic  techniques 
could  be  employed  for  detailed  observation  of  the  processes  in¬ 
volved.  This  approach,  in  conjunction  with  supporting  analytical 
efforts,  should  lead  to  an  experimentally  verified  theoretical 
model  of  the  mixing  region  that  could  subsequently  be  used  to 
optimize  the  design  of  mixing  regions  for  operational  high-power 
chemical  laser  devices. 

Two  experimental  facilities  capable  of  closely  simulating 
the  operational  conditions  of  existing  continuous  flow  chemical 
laser  devices  have  been  described.  Each  of  these  facilities  uses 
the  reflected  shock  region  of  a  conventional  shock  tube  as  a 
supply  reservoir  for  supersonic  two-dimensional  nozzles.  One 
facility  is  designed  for  investigating  the  mixing  characteristics 
resulting  from  the  sonic  injection  of  cold  hydrogen  into  a  super¬ 
sonic  fluorine  carrying  stream.  The  second  facility  is  designed 
for  observing  the  shear  layer  mixing  between  two  parallel  supersonic 
streams  when  one  stream  contains  fluorine  and  the  other  contains 
hydrogen.  It  has  a  unique  split  driven  tube  construction  that  pro¬ 
duces  two  separate  reflected  shock  region  reservoirs  to  supply  two 
adjacent  two-dimensional  nozzles.  Both  facilities  include  viewing 
windows  with  fields  of  view  that  encompass  the  mixing  regions  and 
the  nozzle  exit  stations. 

In  our  opinion,  the  results  presented  in  this  report  clearly 
demonstrate  that  the  experimental  facilities  are  capable  of  simu¬ 
lating  the  mixing  layer  characteristics  of  either  a  hydrogen  in¬ 
jection  or  a  parallel  mixing  continuous  flow  laser  device.  Estab¬ 
lished  experimental  techniques  have  been  adapted  to  each  of  these 
facilities  to  permit  the  mapping  of  critical  flow-field  variables 
within  the  mixing  layer.  An  analytical  procedure  has  been  de¬ 
veloped  that  predicts  these  flow-field  properties  in  the  particular 
case  of  parallel  mixing  and  the  predictions  are  in  good  agreement 
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with  the  experimental  results  for  two  widely  different  cases. 

These  results  demonstrate  that  the  objective  of  an  experimentally 
verified  analytical  procedure  for  detailed  evaluation  of  the  re¬ 
active  mixing  region  within  the  optical  cavity  of  a  CW  chemical 
laser  device  is  realistically  obtainable. 

Recommendations  for  further  work  include  — 

1)  extension  of  the  analytical  model  to  the  injection  mixing 
configuration; 

2)  extension  of  the  experimental  program  to  include  wider 
selection  of  experimental  configurations,  including  close 
duplication  of  operational  device  geometry; 

3)  addition  of  optical  gain  measurement  to  the  diagnostics 
of  the  mixing  region. 
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